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1^20.  Continued. 

""In  this  report,  a  general  theory  is  developed  for  determining  the  complex 
reflection  and  transmission  coefficients  of  such  gratings  as  a  function  of 
incident  angle,  wavelength  and  polarization..  The  theory  does  not  require 
fhat  the  Bragg  condition  be  either  approximately  or  exactly  satisfied,  nor 
does  it  place  any  restrictions  on  the  degree  to  which  the  undiffracted  wave 
is  depleted. 

V 

--Theorems  are  stated  and  proved  about  polarization  effects  and  about  the 
relationship  between  the  spectral  and  angular  characteristics  of  such  grat¬ 
ings^,  Compact,  closed-form  expressions  for  the  complex  reflection  and 
-trlnsmission  coefficients  are  derived.  Strong  similarities  are  shown  to 
exist  between  these  gratings  and  conventional  multilayer  dielectric  stacks. 
For  a  particular  sinusoidal  grating,  our  results  are  compared  with  those 
derived  using  coupled  wave  theory. 

"'-'The  inverse  scattering  problem  is  considered.  Given  the  spectral  and  ang¬ 
ular  characteristics  of  a  horizontally  stratified,  quasi-sinusoidal ,  phase 
grating,  first  and  second  order  approximate  techniques  are  developed  for 
determining  the  corresponding  refractive  index  profile.  v  Examples  of  the 
procedures  are  gi yen.  -  -  _ . 

C^Tt  is  also  shown  that  horizontally  stratified  phase  gratings  can  be  used 
as  position-invariant,  frequency  domain  filters  in  coherent  optical  systems. 
Their  basic  filtering  properties  are  presented. 

Finally,  an  interferometric  technique  is  developed  for  fabricating  hori¬ 
zontally  stratified,  quasi-sinusoidal,  phase  gratings .^The  technique  uses 
a  spatially  incoherent  source  which  has  a  tailored  intensity  profile. 
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PREFACE 

This  final  report  was  prepared  by  the  Electro-Optics 
Department,  Radar  and  Optics  Division  of  the  Environmental 
Research  Institute  of  Michigan.  The  work  was  sponsored 
by  the  Air  Force  Office  of  Scientific  Research  under  Contract 
No.  F49620-78-C-  0061. 

This  report  covers  work  performed  between  1  May  1978  and 
30  September  1980.  The  contract  monitor  is  Dr.  John  A.  Neff, 
Directorate  of  Electronics  and  Solid  State  Sciences,  AFOSR/NE 
Building  410,  Bolling  Air  Force  Base,  Washington,  D.C.,  20332. 
Principal  Investigator  for  this  work  is  K.  Winick.  Major 
contributors  to  the  effort  are  Kim  Winick,  B.  Jin  Chang  and 
W.S.  Colburn. 
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1 

INTRODUCTION 

Holographic  optics  is  continuing  to  find  important  applications 
ranging  from  laser  scanning  systems  to  combiners  for  aircraft  head-up 
displays.  A  holographic  optical  element  is  produced  by  recording  an 
optical  interference  pattern  in  a  photosensitive  material.  When  the 
photosensitive  material  is  processed  following  exposure,  the 
recorded  interference  pattern  is  transformed  into  a  spatial 
variation  of  refractive  index  and/or  a  conductivity  within  the 
material  [1].  The  refractive  index  and  conductivity  will  vary 
significantly  over  distances  on  the  order  of  an  optical  wavelength. 
Therefore,  optical  propagation  through  a  holographic  element  will  be 
governed  by  the  laws  of  diffraction  rather  than  those  of  refraction 
or  reflection.  This  diffractive  optics  technology  is  an  emerging 
area  which  has  potential  for  important  new  optical  devices. 

When  the  spatial  variations  of  refractive  index  and  conductivity 
within  the  material  are  known,  it  is  then  necessary  to  determine  the 
field  diffracted  from  the  holographic  element  given  some  arbitrary, 
incident  optical  field.  For  completely  general  spatial  variations 
of  refractive  index  and  conductivity,  this  diffraction  problem  is 
extremely  difficult.  Many  researchers  have  put  forth  a  considerable 
amount  of  effort  toward  finding  a  solution  to  this  diffraction 
problem  for  limited  classes  of  refractive  index  and  conductivity 
variations.  As  part  of  this  report,  we  solve  the  diffraction 
problem  for  a  particular  class  of  refractive  index  and  conductivity 
variations. 

The  classes  previously  considered  all  have,  with  a  few  exceptions 
to  be  noted  later,  the  following  characteristics: 

1.  The  holographic  element  is  of  infinite  extent  in  the  x  and  y 
directions  (see  Figure  1-1). 
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Refractive  Index  =  n(x,  y,  z)  =  n(x  sin  *  + 


Figure  1-1.  Planar  Holographic  Optical  Element 


z  cos  i) 
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2.  The  element  is  flat,  i.e.,  its  front  and  rear  surfaces  are 
paral lei . 

3.  The  refractive  index  and  conductivity  variations  can  be 
written  as  a  function  of  (cos  0)z  +  (sin  0)x  where  0  is  an 
arbitrary  constant.  This  is  equivalent  to  saying  that  the 
refractive  index  and  conductivity  varies  only  as  a  function 
of  the  distance  measured  along  a  line  which  lies  in  the  x-z 
plane  of  the  element.  This  line  subtends  an  angle  of  0  with 
the  z  axis. 

4.  The  refractive  index  and  conductivity  variations  are 
pt.-  :cdic. 

Thus,  the  structures  which  have  been  stuoieo  are  planar,  periodic 
(though  not  necessarily  sinusoidal)  gratings.  The  methods  developed 
to  date  for  analyzing  such  structures  are  of  three  types  -  rigorous 
modal,  coupled-wave,  and  thin  grating  decomposition. 

The  rigorous  modal  approach  was  first  applied  by  C.B.  Burckhardt 
to  an  element  having  a  sinusoidal  refractive  index  variation  along 
the  x  direction  [2].  The  approach  -leads  to  an  eigenvalue  problem 
formulated  in  terms  of  a  matrix  equation.  The  matrix  is  infinite 
but  is  truncated  in  order  to  make  the  problem  solvable.  By  allowing 
tne  truncated  matrix  to  be  of  sufficiently  large  size,  good  approxi¬ 
mate  results  may  be  obtained.  F.G.  Kaspar  showed  that  the  modal  ap¬ 
proach  could  also  be  applied  to  periodic  gratings  having  spatial 
variations  of  refractive  index  and  conductivity  along  the  x  direction 
which  are  not  sinusoidal  [3].  Burckhardt's  and  Kasper's  numerical 
methods  are  rigorous.  They  are,  however,  beset  by  a  number  of  mathe¬ 
matical  problems. 

H.  Kogelnik  was  the  first  to  analyze  holographic  planar  gratings 
using  coupled-wave  methods.  In  a  classic  paper,  he  drived  closed 
form  expressions  for  the  field  diffracted  by  an  element  having  a 
spatial  variation  of  refractive  index  and  conductivity  which  is 
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sinudoidal,  out  arbitrary  in  direction,  i.e.,  t>  is  arbitrary  [4]. 
Kogelnik's  coupled  wave  theory  makes  a  number  of  approximations. 
These  approximations  are  valid  provided  that  the  Bragg  condition  is 
approximately  satisfied,  and  that  the  undiffracted  field  is  strongly 
depleted.  R.  Magnusson  and  T.K.  Gay  loro  have  extended  the 
couplea-wave  theory  approach  to  cases  where  the  gratings  are  periodic 
but  nonsinusoidal  [5].  Their  method  results  in  an  infinite  system 
of  coupled  wave  equations.  The  system  is  truncated  to  a  finite  size 
in  orcier  to  ootain  a  solution.  Kogelnik  has  also  studied  guided  wave 
devices  which  have  quasi-sinusoidal  perturbations  (i.e.,  the  ampli¬ 
tude  and  frequency  of  the  perturbations  slowly  change  from  one  cycle 
to  the  next)  along  the  guide  in  the  direction  of  wave  propagation 
[6]. 

The  thin  grating  decomposition  method  was  deveU pea  by  R.  Alfer- 
ness  [7].  In  this  technique,  the  holographic  element  is  divided  into 
thin  slabs,  parallel  to  tne  surface  of  the  element.  Each  slab  is 
treated  as  a  thin  transmission  grating.  A  transfer  matrix  of  infi¬ 
nite  size  relating  the  field  at  the  front  surface  of  the  slab  to  that 
at  the  rear  is  postulated.  When  strong  Bragg  effects  are  present,  a 
plausibility  argument  is  presented  for  truncating  the  matrix  to  a 
finite  size.  The  transfer  matrix  of  the  entire  element  is  then  ob¬ 
tained  by  multiplying  together  the  truncated  transfer  matrices  of 
the  individual  slabs.  The  theory  is  developed  only  for  a  TE  pola¬ 
rized  incident  field,  and  a  pure  dielectric  element  (conductivity  is 
zero).  Alferness  applies  this  technique  to  nonperiodic  gratings  as 
well  as  to  multiple  overlapping  gratings  having  grating  vectors  in 
different  directions  (though  still  in  the  x-z  plane).  The  theory  is 
nonrigorous.  W.  Wang  postulates  that  Alferness's  approach  can  be 
extended  to  reflection  gratings  [8].  Prior  to  Alferness's  work,  D. 
Kermisch  also  studied  a  thin  grating  uecomposition  technique  [9]. 
He  applied  his  matrix  method  to  structures  having  spatial  variations 
of  refractive  index  and  conductivity  along  either  the  x  or  z  direc¬ 
tions.  These  variations  were  assumed  to  have  the  form  of  a  slowly 


12 


RADAR  AND  OPTICS  DIVISION 


amplitude  modulated  sinusoid.  The  method  is  only  applicable  at  or 
near  the  Bragg  angle  and  under  the  conditions  that  strong  Bragg 
effects  dominate  the  diffraction  process.  All  of  the  above  theories 
require  that  the  incident  field  have  no  spatial  variations  in  the  y 
direction . 

D.  Young  nas  developed  a  coupled-wave  theory  tor  analyzing  the 
diffracted  field  from  a  holographic  element  having  arbitrary  spatial 
variations  of  refractive  index  and  conductivity  (10].  In  his  ap¬ 
proach,  infinite  sets  of  variable  coefficient  coupled  wave  equations 
are  obtained.  He  retains  only  a  single  pair  under  the  assumption 
that  strong  Bragg  effects  dominate  the  diffraction  process,  and  sol¬ 
ves  them  using  the  WKJB  approximations.  No  restriction  is  placed  on 
the  incident  field. 

In  summary,  we  make  the  following  observations.  Rigorous  modal 
methods  are  applicable  only  to  periodic  structures.  They  entail  the 
numerical  solution  of  a  complicateo  eigenvalue  problem.  Furthermore, 
the  eigenvalue  problem  is  of  infinite  dimension  and  must  be  truncated 
to  a  finite  size  in  order  to  obtain  a  solution. 

Coupled  wave  methods  have  been  applied  to  both  periodic  and  non¬ 
periodic  structures.  In  the  case  of  periodic  structures,  an  infinite 
set  of  constant  coefficient  coupled  wave  equations  is  obtained.  This 
set  is  truncated  to  a  finite  size  and  tne  resulting  equations  are 
solved  numerically.  For  non-periodic  structures,  an  infinite  set  of 
vanaule  coefficient  coupled  wave  equations  is  obtained.  Unaer  the 
assumption  that  strong  Bragg  effects  dominate  the  diffraction  pro¬ 
cess,  this  set  is  truncated  to  just  two  equations.  The  two  variable 
coefficient  coupled  wave  equations  are  solved  using  the  WKJB 
approximation . 

Thin  grating  decomposition  net hods  have  little  theoretical  jus¬ 
tification.  The  technique  requires  the  truncation  of  an  infinite 
transfer  matrix  to  a  finite  size.  The  method  does  not  consider  the 
case  of  TM  polari zation. 
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This  report  only  considers  structures  which  have  the  properties 
listed  below: 

1.  The  structure  is  a  pure  dielectric  (conductivity  is  zero). 

?.  The  structure  is  of  infinite  extent  in  the  x  arid  y 
dimensions . 

3.  The  refractive  index  of  the  structures  varies  spatially  only 
as  a  function  of  z.  (Structures  of  this  type  are  said  to  be 
horizontally  stratified.) 

4.  The  spatial  variations  of  the  refractive  index  are  quasi- 
sinusoidal . 

Holographically  produced  dielectric  structures  which  are  horizontally 
stratified  are  a  type  of  phase  reflection  gratings.  We  wish  to  de¬ 
termine  the  field  diffracted  from  such  a  structure  given  an  arbitrary 
incident  optical  field.  It  is  desirable  to  develop  an  analysis  tech¬ 
nique  which  is  rigorous,  computationally  simple,  and  allows  one  to 
deduce  general  properties  of  these  structures.  The  modal  technique 
cannot  be  applied  to  quasi-sinusoidal  structures.  Although  the 
coupled  wave  theory  approach  can  be  applied  to  quasi-sinusoidal 
structures  and  is  semi-rigorous,  it  is  only  justified  under  the  con¬ 
dition  that  strong  Bragg  effects  dominate  the  diffraction  process. 
It  is  not  immediately  apparent  whether  general  properties  of  the 
quasi-sinusoidal  structures  can  be  deduced  by  this  coupled-wave  ap¬ 
proach.  The  thin  layer  decomposition  method  is  non-rigorous.  Fur¬ 
thermore,  it  is  restricted  to  the  case  where  the  incident  field  has 
TE  polarization.  Thus,  none  of  the  existing  grating  theory  tech¬ 
niques  are  well  suited  for  determining  the  diffracted  field  from  the 
structure  being  considered  here. 

Electromagnetic  propagation  through  horizontally  stratified  media 
has  been  studied  extensively  [11].  A  general  theory  was  first  pub¬ 
lished  in  1950  by  Abeles  [12].  The  theory  has  been  successfully  ap¬ 
plied  to  a  number  of  areas  including  radio  wave  propagation  through 


14 


RADAR  AND  OPTICS  DIVISION 


the  ionosphere  ano  multilayer  dielectric  stacks.  Researchers  in  the 
area  of  holographic  gratings  have  neglected  Abele's  theory  for  the 
most  part.  W.W.  Rigrod  uid,  however,  draw  an  analogy  between  thick 
hologram  gratings  and  multilayer  dielectric  stacks  [13],  while  H. 
Koyelnik  has  mentioned  the  possibility  of  using  the  theory  of  strati¬ 
fied  media  to  study  thick,  holographic,  reflection  gratings  [6]. 

Research  has  not  been  pursued  in  this  area  for  fundamental  rea¬ 
sons.  First,  and  perhaps  foremost  is  the  fact  that  fabrication  of 
holographic  reflection  gratings  is  considerably  more  difficult  than 
fabricating  other  types  of  holographic  gratings.  This  is  because 
faoricating  reflection  gratings  requires  recording  n.atei ials  of  very 
high  resolution  (>5000  line  pairs/mm).  Secondly,  applications  of 
holograpnic  reflection  gratings  havt  developeu  more  si  .wly  tnan 
those  of  other  other  holographic  gratings.  Reflection  gratings  are, 
however,  becoming  increasingly  more  important.  The',  are  now  an  inte¬ 
gral  part  of  distributed  feedback  lasers  and  are  also  used  as  selec¬ 
tive  wavelength  filters  in  integrated  optics  devices  and  aircraft 
head-up  displays  [14-15].  Thirdly,  the  theory  of  stratified  media 
is  only  applicable  to  a  limited  class  of  gratings,  whereas  the  em¬ 
phasis  today  tends  to  be  on  developing  general  theories  which  have 
broader  ranges  of  applicability.  Finally,  no  serious  consideration 
has  been  given  to  the  problem  of  producing  holographic  reflection 
gratings  having  tailored  quasi-sinusoida  I  refractive  index  profiles. 
There  exists  the  possibility  of  designing  such  gratings  so  that  they 
have  specified  spectral  and  ungutar  characteristics.  The  ability  to 
design  these  gratings  would  require  extremely  accurate  analysis  tech¬ 
niques.  Since  the  grating  design  problem  has  not  been  seriously  con¬ 
sidered  to  date,  such  analysis  techniques  have  not  been  necessary. 

The  theory  of  stratified  media  provides  a  completely  rigorous 
method  of  determining,  to  any  degree  of  accuracy,  the  field  dif¬ 
fracted  by  the  structures  under  consideration.  Therefore,  it  pro¬ 
vides  a  means  for  checking  the  validity  and  accuracy  of  all  other 
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theories  as  they  apply  to  these  structures.  We  will  show  how  it  may 
also  be  used  to  gain  a  fundamental  understanding  of  the  spectral  and 
angular  characteristics  of  reflection  gratings.  Furthermore,  it  will 
aid  in  the  development  of  approximate  design  techniques.  A  design 
technique  is  a  procedure  which  allows  one  to  determine  the  refractive 
index  profile  of  a  reflection  grating  given  its  spectral  and  angular 
characteristics.  This  is  an  inverse  scattering  problem. 

This  report  is  organized  into  nine  chapters.  Chapter  two  dis¬ 
cusses  the  applicability  of  frequency  domain  concepts  to  the  analysis 
of  holographic  reflection  gratings.  The  use  of  these  gratings  as 
position-invariant  frequency  domain  filters  is  discussed. 

Chapter  three  reviews  the  general  theory  of  electromagnetic  pro¬ 
pagation  through  horizontally  stratified  media.  This  material  is 
not  new  and  is  presented  only  in  the  interest  of  completeness.  A 
computer  routine  is  written  to  calculate  the  field  diffracted  from  a 
horizontally  stratified  media.  Two  examples  are  given,  one  of  a 
multilayer  dielectric  stack  and  the  other  of  a  holographic  grating 
having  a  sinusoidal  refractive  index  profile.  The  field  diffracted 
from  the  holographic  reflection  grating  is  also  calculated  using 
coupled  wave  theory.  The  two  results  are  compared. 

Chapter  four  uses  techniques  from  the  theory  of  horizontally 
stratified  media  to  develop  an  efficient  method,  based  on  matrix 
multiplication,  for  calculating  the  field  diffracted  from  holographic 
reflection  gratings. 

Chapter  five  applies  the  theory  of  the  previous  chapter  to  derive 
analytic  expressions  for  the  field  diffracted  from  holographic  re¬ 
flection  gratings  having  quasi-sinusoidal  refractive  index  profiles. 
These  expressions  are  used  in  three  examples  to  determine  the  complex 
reflection  and  transmission  coefficients  of  gratings.  A  closed  form 
expression  is  also  derived  for  the  complex  reflection  and  trans¬ 
mission  coefficients  of  a  holographic  reflection  grating  having  a 
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sinusoidal  refractive  index  profile.  This  expression  is  universally 
valid.  It  does  not  require  that  the  Bragg  condition  be  satisfied 
either  exactly  or  approximately,  nor  does  it  require  strong  Bragg 
effects  to  dominate  the  diffraction  process. 

Chapter  six  discusses  some  general  properties  of  holographic  re¬ 
flection  gratings.  The  mutual  interaction  of  two  reflection  gratings 
placed  in  contact  is  investigated,  as  is  the  connection  between  the 
spectral  and  angular  character istics  of  a  grating.  A  theorem  on 
polarization  effects  is  stated  and  proved.  Using  the  matrix  multi¬ 
plication  technique  developed  in  chapters  four  and  five  combined  with 
a  mathematical  induction  argument,  general  expressions  for  the  com¬ 
plex  reflection  and  transmission  coefficients  are  derived  which  do 
not  require  matrix  multiplications  for  their  evaluation.  These  ex¬ 
pressions  are  used  in  the  following  chapter  to  develop  approximate 
techniques  for  solving  the  inverse  scattering  problem. 

Chapter  seven  considers  the  inverse  scattering  (design)  problem. 
Given  the  spectral  and  angular  characteristics  of  a  holographic  re¬ 
flection  grating,  first  and  second  order  approximate  techniques  are 
developed  for  determining  the  corresponding  refractive  index  varia¬ 
tion.  Examples  of  the  procedure  are  given. 

Chapter  eight  develops  and  analyzes  interferometric  techniques 
for  fabricating  holographic  reflection  gratings  having  quasi- 
sinusoidal  refractive  index  profiles. 

Chapter  nine  contains  our  conclusion,  and  a  discussion  of  areas 
for  future  research. 
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2 

FREQUENCY  DOMAIN  CONCEPTS 

In  this  chapter,  we  demonstrate  the  applicability  of  frequency 
domain  concepts  to  the  analysis  of  flat,  holographic,  reflection 
qratinqs.  We  restrict  our  discussion  to  gratings  which  have  refrac¬ 
tive  index  variations  only  in  the  z  direction  (see  Figure  2-1). 

'he  field  transmitted  by  such  a  grating,  when  it  is  illuminated 
by  a  plane  wave,  will  be  another  plane  wave  propagating  in  the  same 
direction  as  the  irrident  plane  wave  pi.  The  phase  and  amplitude 
of  this  transmitted  plane  wave  will,  however,  differ  from  the  inci¬ 
dent  one.  The  field  reflected  by  such  a  grating,  wher  it  is  illumi¬ 
nated  by  a  plane  wave,  will  be  another  plane  wave.  Its  x  and  y  di¬ 
rectional  cosines  will  be  identical  with  those  of  the  ircident  plane 
wave,  and  its  z  directional  cosine  will  be  of  opposite  sign.  As  in 
the  case  of  transmission,  the  nhase  and  amplitude  of  the  reflected 
plane  wave  will  differ  from  the  incident  one. 

Any  monochromatic  wavefront  can  be  described,  through  a  Fourier 
transform  decomposition,  as  the  superposition  of  plane  waves.  A 
flat,  holographic,  reflection  grating  is  a  linear  device,  therefore 
the  field  scattered  by  it  (i.e.,  the  transmitted  and  reflected 
fields)  can  written  as  the  sum  of  the  scattered  fields  resulting 
from  each  plane  wave  component  of  the  incident  field.  Thus,  linear 
system  theory  can  he  applied,  and  the  input-output  transfer  function 
of  a  flat,  holographic,  reflection  grating  derived. 

From  the  above  discussion,  we  note  that  the  character istics  of 

the  holographic  grating  can  be  completely  described  by  two  frequency 

domain  transfer  functions,  fy)  and  HR(fx,  *y)* 

H—( f  ,  f  )  is  the  transfer  function  for  the  transmission  case, 
I  x  y 

and  H„(f  ,  f  )  is  the  transfer  function  for  the  reflection 
R  x  y 

case.  Roth  HR(fx,  fy)  and  Hj^x’  are  comPlex  anrl 

may  be  written  as 
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Refractive  Index  =  n(z) 


Figure  2-1.  Holographic  Reflection  Grating 
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HT(fX-  V  ■  lHT(fx-  V1  eXP  (iVfx'  fy» 

WV  fy>  -IW  V1  «Vfx-  fy"  <«> 

HT(f  ,  f  )  is  the  ratio  of  the  complex  electric  field  ampli- 
I  x  y 

tude  of  the  transmitted  wave  to  the  complex  electric  field  amplitude 

of  the  incident  wave.  H0(f  ,  f  )  is  the  ratio  of  the  complex 

k  x  y 

electric  field  amplitude  of  the  reflected  wave  to  the  complex  elec¬ 
tric  field  amplitude  of  the  incident  wave,  xf  and  xf  are  the 

x  y 

x  and  y  directional  cosines,  respectively,  of  the  incident  vave  (x 
is  the  wavelength  of  the  incident  wave). 

Initially,  our  discussion  will  be  restricted  to  the  transmission 
case,  the  refection  case  is  analogous.  It  is  also  assumed  that  we 
are  working  in  spatially  and  temporally  coherent  light.  Consider 
the  situation  shown  in  Figure  2-2.  A  wavefront  is  propaaatinq  to 
the  right.  In  plane  one,  this  wavefront  has  a  spatial  distribution 
of  amplitude  and  phase  given  by 

at(x,  y)  =  |  a  j  ( x ,  y)|  exp  (i  ( x,  y))  (2-3) 


I n  plane  one  this  wavefront  can  be  equivalently  described  in  the 
frequency  domain  by  [2] 

w  y  =  lyv  yi exp  (ivy  y>  (?-4) 

where 


Aj(f 


x’ 


(x,  y)  exp[-i2ir(fxx  +  fyy)]dxdy (2-ba ) 


3j(x,  y) 


no  oo 


x’ 


fy)  exp[1^-(Xfxx  +Xfyy)dfxdfy  (2-5b) 


The  wavefront  in  propagating  between  plane  one  ana  the  front  surface 
of  the  holographic  grating  traverses  free  space.  The  transfer  fi no¬ 
tion  describing  free  space  propagation  between  two  planes  separated 
by  a  distance  d  is  [21 
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Plane  1  Holographic  Plane  2 

Grating 


W 


V 


Figure  2-2.  Transmission  through  a  Holographic  Reflection 
Grating 
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SH(f  ,  f  ;  d)  =  exp[-iuxd(f  2  -r  f  2)]  (2-6) 

A  y  a  y 

Thus,  the  wavefront  at  the  front  surface  holographic  grating  is  given 
in  the  frequency  domain  by 

AIffx»fy)<SW( fx»  fy;  d?)  (2-7) 

and  at  the  back  surface  by 

Vfx-fy>®<fx-  V  Wfx’  V  <M> 

The  wavefront  once  again  propagates  through  free  space  between  the 
back  surface  of  the  holographic  grating  and  plane  two.  The  wavefront 
at  plane  two  is  described  in  the  frequency  domain  by 

AI(fx’fy)<3U(fx’  fy;  dl)HT(fx’fy}^(fx’  V  d?}  U_9) 

We  note  that 

^fx’fy^^V  fy;  d2}  “®(fx»fy;  dl  +  d2}  (?"10) 

Thus  Fq.  (2-9)  can  be  rewritten  as 

W  fy)HT(fx’fyWfx’  V  dl  +  d2}  (?-H) 

Eq.  (2-11)  indicates  that  the  transfer  function  relating  plane  one 
to  plane  two  is 

W(f  ,  f  ;  d,  +  d,)H_(f  ,  f  )  (2-12) 

x  y  I  2  T  x  y 

The  following  conclusions  can  be  drawn  from  Eq.  (2-12) 

1.  The  holographic  grating  acts  as  a  frequency  domain  filter. 
It  changes  both  the  amplitude  and  phase  of  the  incident  fre¬ 
quency  components. 

2.  The  relationship  between  the  wavefront  at  plane  two  and  the 
wavefront  at  plane  -  ?  depends  only  upon  Hy  ( f  x  *  %) 
and  the  distance  d,+d0  between  plane  one  and  plane  two. 
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This  relationship  does  not  oepend  upon  tne  value  of  o,  or 
^2  alone.  For  example,  consider  tne  two  systems  shown  iu 
Figure  2-3.  If  the  wavefront  at  plane  one  is  identical  in 
both  systems,  then  Doth  systems  will  have  identical  wave¬ 
fronts  at  plane  two. 

We  investigate  the  use  of  the  holograpnic  grating  as  a  frequency 
domain  filter.  Suppose  the  holographic  grating  is  placed  immediately 
behind  plane  one  (i.e.,  a-|  =  0),  then  by  Eqs.  (2-6),  ('>- 1 )  and 
(2-8)  the  frequency  domain  description  of  the  wavefront  at  the  back 
surface  of  tne  element  will  oe 


¥V  VW  V  (2-13! 

Tne  hack  surface  of  the  holographic  grating  can  be  imaged  using  con¬ 
ventional  refractive  optics  onto  a  distant  plane  [3].  The  frequency 
domain  relationship  between  the  wavefront  at  plane  one  and  tne  wave- 
front  at  this  distant  image  is  given  by  Eq.  (2-14) 

v-v  -yy-v  -y  <2-i4> 


We  note  that  in  the  spatial  domain  the  wavefronts  given  by  Eqs. 
(2-13)  aid  (2-14)  are  conjugates.  This  conjugation  always  results 
from  conventional  refractive  optics  of  imaging  systems.  The  dis¬ 
tinction  is  actually  meaningless  in  tne  image  plane  where  i.he  inten- 

★ 

si  ty  of  conjugate  wavefronts  are  identical.  Thus,  H-p  ( -t  x , 
-f  )  is  the  transfer  function  between  plane  one  and  the  image 
plane.  A  system  for  performing  frequency  domain  filtering  is  shown 
in  Figure  2-4,  where  we  represent  the  conventional  refractive  optics 
imaging  system  by  a  single  lens.  It  should  be  noted  that  conclusion 
two  auove  implies  that,  the  holographic  grating  can  ue  located  any¬ 
where  between  plane  one  arid  the  refractive  optics  imaging  system 
without  Changing  the  system  transfer  function.  Therefore,  the  system 
in  Figure  2-5  has  the  same  transfer  function  as  the  system  in  Figure 
2-4. 
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Plane  1  Holographic 
Grating 

H,(f  ,  f) 
T  x  y 


Plane  2 


f 


Plane  1  Holographic  Plane  2 

Grating 


w  v 


Figure  2-3.  Position  Invariant  Property 
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Plane  1 
( Input) 


1 


Hoi ographic 
Grating 


M\.  K) 


I 

Image  Plane 
(Output) 


Imaging  Optics 
( Images  Plane  1 
onto  the  Image  Plane) 


Figure  2-4.  Single  Grating  Frequency  Domain  Filtering  System 
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Plane  1 
( Input) 


Holographic 

Grating 


V 


Image  Plane 
( Output ) 


Imaging  Optics 
( Images  Plane  1 
onto  the  Image  Plant) 


Figure  2-5.  Position  Invariant  Property  of  Frequency  Domain 
Filtering  System 
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We  have  described  a  system  which  is  capable  of  performing  fre¬ 
quency  domain  filtering.  The  principal  component  of  this  system  is 
a  flat,  holographic,  reflection  grating  which  has  a  refractive  iroex 

variation  only  in  a  direction  normal  to  its  surface.  The  transfer 

★ 

function  of  this  system  is  given  by  H^-f^,  -f  ).  It  is 

natural  then  to  investigate  the  characteristics  of  H,(f  ,  f  ). 

i  x  y 

From  symmetry  considerations  it  is  clear  that  HT(f  ,  f  ) 

I  x  y 

is  circularly  symmetric,  i.e., 

W  V  ■  Hi(Vfx  *  fy  )  <?-|5) 

Also,  it  should  be  noted  that  F  (f  ,  f  )  is  a  complex  quantity 

i  x  y 

as  is  indicated  by  Eq.  (2-1).  Given  the  refractive  index  profile, 

n(z),  it  is  possible  to  determine  HT(f  ,  f  ).  We  have  per- 

formed  such  a  calculation  (actually  we  computed  HR(fx,  f  )  for 

the  case  where  n(z)  is  sinusoidal).  A  holographic  grating  having 

such  a  refractive  index  profile  can  be  fabricated  by  recording  the 

interference  pattern  of  two  plane  waves  as  shown  in  Figure  2-6.  The 

HD(f  ,  f  )  corresponding  to  this  element  is  shown  in  Figures 
K  x  y 

2-7  and  2-8.  The  relationship  between  the  e  given  in  Figures  2-7 

and  ?-8  and  (f  ,  f  )  is 
x  y 

cos  e0  =-0  -  (Xfx)2  -  (Xfy)2  (2-16) 

Eq.  (2-15)  can  therefore  be  written  as 

HT(fx,  fy)  =  Ht(»1,  e0)  (2-17) 

Consider  the  setup  shown  in  Ficere  2-(;.  An  input  transparency 
is  illuminated  by  a  collimated  beam  at  normal  incidence.  The  light 
transmitted  by  the  transparency  passes  through  a  beamsplitter  and 
then  strikes  the  holographic  grating  shown  earlier  in  Figure  2-6. 
The  light  reflected  by  the  holographic  grating  is  redirected  6v  t> e 


UflVELENGTH  =  0.5145  MICRONS 
S  POLflRIZflRION 


n(z)  =  1.5  +  0.01  sin  2  t,  (5.8309)z 
z  is  in  microns 

Total  hologram  thickness  =  15  urn 
Hologram  imbedded  in  medium  having 
refractive  index  of  1.5 


Figure  2-8.  Phase  Change  Upon  Reflection  vs.  Angle  of 
Incidence  Characteristic 
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Figure  2-9.  Holographic  Mirror  System 
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beamsplitter  to  pass  through  some  refractive  optics  which  image  the 

input  transparency  onto  an  output  plane.  Following  our  earlier 

analysis,  we  designate  the  spatial  distribution  of  amplitude  and 

phase  of  the  wavefront  at  the  back  surface  of  the  transparency  by 

aj(x,  y).  The  frequency  domain  description  of  this  wavefront,  is 

A. ( f  ,  f  ).  In  the  output  plane  the  frequency  domain  uescrip- 

l  x  y 

tion  of  the  wavefront  is 

h<-v  -y  v-v  -v  (2-'8) 

where  HR(fx,  f  )  is  the  transfer  function  of  the  holographic 

grating.  Typically  a  transparency  will  have  spatial  frequency  com¬ 
ponents  up  to  100  n/rm.  Thus  we  ran  write 

Ai(fx,  fy)  =  0  for  fy  >  100  fc/mm  or  fy  >  100  ji/mm  (2-19) 

If  A  =  5145A,  then  by  Eqs.  (2-16)  through  (2-19)  it  is  clear  that 
the  only  segment  of  Hr(a,  ©o)  which  will  affect  the  wavefront  in 
the  output  plane  is  the  segment  described  by 

0°  <  ©  <  cos"1^!  -  (0.5145)Z(0.1)Z  -  (0.5145)2(0.1  )2)=  4.2°(2-20) 

From  Figure  2-7,  it  can  be  seen  that  |Hr(a,  ©  )|  is  reasonably 

flat  between  0°  and  4.2°.  Thus,  this  particular  holographic  grating 
does  not  change  the  relative  amplitudes  of  the  input  transparencies 
frequency  components.  Next  we  examine  the  exp  [ i <l>D ( x ,  9  )] 
term  of  the  frequency  response.  $R  vs  eQ  varies  by  about  30° 
(or  0.083  wavelengths)  in  the  range  between  0"  <_  ©Q  <  4.2°.  For  © 
small,  we  note  for  completeness  that  Eq.  (2-16)  can  be  written  ap¬ 
proximately  as 

©  =  xfx  +  xfy  (2-21) 

where  eQ  is  in  radians.  Small  variations  in  the  phase  change  upon 

reflection  as  a  function  of  (f  ,  f  )  will  introduce  aberrations 

x  y 
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into  the  image,  though  these  may  not  be  readily  discernible.  A  lin¬ 
ear  phase  shift  upon  reflection,  measured  as  a  function  of  (f  , 
f  )  has  no  effect  on  the  image  formed  in  the  output  plane  except 
that  of  translation . 

Most  individuals  familar  with  holographic  optical  elements  would 
expect  the  system  shown  in  Figure  2-9  to  behave  as  a  mirror.  it 
should  be  noted  however  this  result  is  actually  very  surprising.  It 
is  surprising  because  this  mirror  like  behavior  requires  the  holo¬ 
graphic  grating  to  display  as  a  function  of  and  f  a  very 
small  or  a  nearly  linear  phase  change  upon  reflection  over  the  prin- 
cipal  reflection  hand  of  the  filter.  We  note  that  a  similar  situa¬ 
tion  occurs  in  the  case  of  multilayer  dielectric  mirrors  [4].  These 
mirrors  are  observed  to  form  virtually  unaberrateri  imagec  upon  re¬ 
flection.  This  requires  that  the  phase  change  upon  reflection  of 
such  a  mirror  be  a  very  small  or  nearly  linear  funrtion  of  f  anc 
f  over  the  principal  ranae  of  frequencies  whirr  the  mirror  re¬ 
flects.  For  an  arbitrary  refractive  index  profile,  n(z), 

( f  ,  f  )  and  4 _ ( f  ,  f  )  are  neither  very  small  nor 
K  x  y  I  x  y 

linear  functions  of  f  anci  f  over  the  entire  ranae  of  t  arc1 

x  y  -  3  x 

f  .  Linder  what  conditions  and  over  what  ranges  are  tn<  y  np 
linear  This  question  remains  to  be  answered.  This  is  an  important 
consideration  when  filtering  is  to  he  per  former  using  this  type  f 
holographic  grating. 

We  explore  further  the  use  of  flat  holngraDhir  gratinos  whir* 
have  refractive  index  variations  only  in  a  direction  normal  to  then- 
surface.  Tn  the  previous  discussion  we  have  assumed  that  the  holo¬ 
graphic  grating  is  oriented  so  that  its  surface  is  parallel  tc  the 
input  plane.  This,  we  recall,  guaranteed  that  the  qratinq  response 

be  circularly  symmetric  in  (f  ,  f  ). 

x  y 

Now  suppose  that  the  grating  is  rotated  by  a  degrees  about  its 
y-axis  and  by  e  degrees  about  the  x-axis  (♦  degrees  ineicat.es 
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clockwise  rotation  while  -  degrees  indicates  counterclockwise  rota¬ 
tion).  Then  we  can  write 


H* (f x ,f y )  =  HyfVl  -sin2(sin~1(xfx)-a)-sin?(sin"1rxfy)+6))  (2-?2) 

where  H“,8(f  ,  f  )  is  the  transfer  function  of  the  grating 

1  x  y 

when  it  is  rotated  by  a  degrees  about  the  y  axis  and  by  e  degrees 

about  the  x  axis.  HT(f  ,  f  )  is  the  transfer  function  of  tnp 

grating  when  it  is  not  rotated  at  all.  It  should  be  noted  that  the 
qratinq,  when  rotated,  has  a  transfer  function  which  is  not  circu¬ 
larly  symmetric.  Furthermore  and  perhaps  more  significantly,  the 
response  of  the  qratinq  is  not  the  same  m  the  f  =  0  Diane  as  in 
the  f  =  0  plane 

Consider  now  the  situation  shewn  in  Fiqure  ?-10.  An  input  plane 
is  imaged  onto  an  output  plane  by  conventional  refractive  optics.  A 
""■’oi  atinp  havina  a  system  transfer  function  H  (f 

distance  d^  from  the  input  plane.  Another  holo- 
having  a  system  transfer  'unction  h  (f  , 
f  a  distance  d9  from  the  output  plane.  A  wav<?ro>t 

is  pi  from  left  to  right  in  this  system.  In  the  frequency 

■jumair,  tho  wave'ront  at  the  input  plane  will  be  described  bv 

A. (f  ,  f  ) 

I  '  x  ’  v 


Let,  us  imaqine  for  the  moment  that  Hy  ( f  x ,  f^ )  =  1.  this  is 

equivalent  to  removing  the  second  holographic  grating  from  the  sys¬ 
tem.  Thm  at  ’he  utput  plane 


★ 


(-* 


( ?-?3 ' 


gives  the  frequency  domain  description  of  the  wavefront.  At  a  dis¬ 
tance  *2  ^rom  the  output  plane,  ’ he  wavetron’  is  de'-rribeu  t  v 


A. ( -f  ,  -f  )HT f -f 
I  x  ’  y '  T  ) 


-f  ISfV  ,  f  ;dj 

y  x  y  2 
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Figure  2-10.  Two  Grating  Filtering  System 
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Thus,  with  Ht  (f  ,  f  )  not  identically  equal  to  1,  the  wave- 

1 9  x  y 

front  at  the  output  plane  will  be  given  by 


-%K/-v*v^V,^'vvv*<v  v0?1 


Hence,  we  rewrite  Eq.  (2-25)  as 


A,  ( -f 


-f  )H,  f-f  -f  )H,  (f v ,  f  ) 
*  y  1  ^  *  y  1 2  *  y 


Therefore, 


h!  f-f  ,  -f  1Ht  (f  ,  f  ) 

T^  x  y  T^  x*  y 


(2-26) 


(2-27) 


is  the  frequency  domain  transfer  function  of  the  system. 

Now  if 

H-  (f  ,  f  )  =  Ht  (-f  ,  -f  )  (2-28) 

Tp  x’  y'  T-|V  x’  y'  '  ' 

Then  Eq.  (26)  becomes 

A*(-f  ,  -f  )|H,  (-f  ,  -f  )|?  (2-29) 

V  x’  yM  T i '  x’  y'l  '  ’ 

We  note  that: 

1.  Eq.  (2-26)  describes  in  the  frequency  domain  the  wavefront 
at  the  output  plane. 

2.  This  '-'avpfront  is  independent  of  the  longitudinal  positions, 
i.e.,  d^  and  n^,  of  the  holographic  gratings 

3.  The  above  equations  remain  valid  in  the  case  where  the  holo- 
grapnical  optical  elements  are  rotated  about  their  x 

and  y  axes  if  H  (f  ,  f  )  is  replaced  by  H*,e(f  ,  f  )  and 
I  j  x  v  1 1  x  y 

and  replaced  bv  H“,e  (f  ,  *  1. 

'9  *  y 
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Jim 


4.  From  Eqs.  (2-22)  and  (2-2 7)  it  is  seen  that  the  frequency 
domain  transfer  function  of  the  system  will  be  real  if  the 
two  holographic  gratinqs  are  identical,  and  if  they  are  ro¬ 
tated  by  equal  amounts  but  in  opposite  directions.  By  rota¬ 
tion  of  equal  amounts  but  in  opposite  directions,  we  mean 
that  a i  =  ~  «2  and  e]  =  -e2' 

5.  Comment  (4)  indicates  that  a  real  transfer  function  can  al¬ 
ways  be  realized  by  using  two  identical  holographic  optical 
elements  which  have  an  arbitrary  and  perhaps  unknown  phase 
transfer  function  characteristic. 

6.  A  situation  analogous  to  what  has  been  described  above  exists 
when  the  holographic  gratings  are  used  in  the  reflection 
rather  than  the  transmission  mode. 
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Notes  to  Chapter  2 

[1 ]  For  simplicity,  we  assume  that  the  grating  is  imbedded  in  mediums 
which  match  it  in  refractive  index  at  its  front  and  rear  surfaces. 

[2]  J.W.  Goodman,  Introduction  to  Fourier  Optics,  McGraw-Hill,  New 
York,  1968. 

[3]  A  telescope  system  could  be  used  to  perform  the  imaging. 

[4]  W.G.  Driscoll,  Handbook  of  Optics,  Optical  Society  of  America, 
McGraw-Hill  Co.,  New  York,  1978. 
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3 

THE  THEORY  OF  DISCRETE  HORIZONTALLY  STRATIFIED  MEDIA 

This  chapter  presents  the  theory  of  electromagnetic  wave  propa¬ 
gation  through  discrete  horizontally  stratified  media.  The  theory 

was  first  published  by  Abeles  in  1950  [lj.  Our  development  is 

similar  to  one  given  by  Knittl  [2],  The  theory  is  presented  in  an 

effort  to  make  this  report  self-contained. 

Given  the  description  of  a  discrete,  horizontally  stratified, 
dielectric  media,  i.e.,  the  number  of  layers,  thickness  of  each 

layer  and  refractive  index  ot  each  layer,  we  will  determine  the 
reflectance,  phase  change  upon  reflection,  transmittance  and  phase 
chariye  upon  transmission  of  the  media  when  illuminated  by  an  s  (TE) 
or  p  (Til)  polarized  plane  wave  of  arbitrary  frequency  at  an 
arbitrary  angle  of  incidence.  In  the  notation  of  chapter  two,  tnis 
is  equivalent  to  determining  Hr[a,  oq)  and  Hj(\,  oq). 

Consider  the  structure  shown  in  Figure  3-1.  It  consists  of  N 
parallel  homogeneous,  isotropic  dielectric  layers  each  assumed  to  be 

f  h 

of  infinite  extent  in  the  x  and  y  directions.  The  j  '  layer  is  of 

thickness  d.  measured  in  microns  and  has  refractive  index  n,... 

J  \  \J  ) 

The  refractive  index,  n(jj  Is  given  by 

V  eo  yo 

where  e/-%  and  y,  >  are  the  perrni tti vity  and  permeability,  re- 
U  1  ^  t  n 

spectively,  of  the  j  layer,  and  eQ  and  yQ  are  tne  permitti¬ 
vity  and  permeability,  respectively,  of  tree  space.  It  is  assumed 
that  equals  yQ.  For  notational  convenience  the  region 

above  layer  1  will  b°  designated  as  layer  0,  while  the  region  be¬ 
low  N  will  be  jesignated  layer  N+l.  We  make  the  following 
det in i tions 
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©<  - ,  =  sin 

(j) 


.-1  "o 


n„  sin  e\ 

A  A  \ 


"(j)  J 


Y(j)  = 


n(j)  C°S  Q(j) 

Vu0U0 


,  for  s  polarization 


,  for  p  polarization 


(cos  9(j))v^ 


A  2 

=  r  n(j)d(j)  cos  e(j) 


A  feXP  (l0(j))  ° 


exp  (-i*(j) 


V)  =  U(N+1) 


A  1  0 


0  1 


M,  •  %  =  V,  .XU,  .XV, - \ 

(j)  (j)  (j)  (j) 


12  A 


n  "u, 


22J  j=0 


In  the  remaining  portions  of  this  chapter,  it  will  be  shown  that 


hr(x,  0q)  * 


^ —  ,  for  s  polarization 
bll 


,  for  p  polarization 
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—  ,  for  s  polarization 
HI 


HT(x,  e0) 


(N+l) 


,  for  p  polarization 


where  x  is  the  wavelength  of  the  incident  phase  wave  and  is  the 
angle  at  which  the  incident  plane  wave  is  propagating  with  respect 
to  the  z  axis. 

Maxwell's  phasor  wave  equation  for  the  electric  field  in  a  homo¬ 
geneous  isotropic  dielectric  is  given  by  Eq.  (3-1)  below 


32E 

32E 

32E 

X 

+ 

X 

+ 

X 

3X2 

ay2 

3Z2 

a2E 

32E 

32E 

y 

+ 

y 

+ 

y 

3X2 

ay2 

3Z2 

2r 

2r 

3  E 

3  E 

3  E 

z 

+ 

z 

+ 

z 

„  2 

,  2 

.  2 

3  X 

ay 

3Z 

+  K  Ex  =  0 


I  +  K6Ez  =  0 


(3-la) 


(3-lb) 


( 3—  1  c ) 


where  K  =  w/Ue.  u  is  the  frequency  in  radians  while  e  and  v  are  the 
permittivity  and  permeability  of  the  dielectric  respectively.  It  is 
noted  that 


F  -  F  +  f 
a  "  xR  lxL 


F  -  F  +  F 

y  "  yR  yL 


Ez  =  EzR  +  Ezl 


(3 -2a) 
(3-ztl 
( 3-2c ) 


where 
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exr 

aeiK((sin  e)y-(cos  e)2) 

( 3— 3a ) 

eyr 

=  0 

( 3— 3b ) 

ezr 

=  0 

( 3— 3c ) 

exl 

heik((sin  p).y+fcos  e)z) 

( 3-4a ) 

eyl 

=  0 

(3-4h) 

ezl 

=  0 

(3-4 c) 

with  a,  h,  and  e  arbitrary  constants  satisfies  the  wave  equations. 

A  solution  of  the  wave  equation  (3-1)  is  a  solution  of  Maxwell’s 

Equations  if  it  also  satisfies  v*(E  x  +  E  y  +  E  z)  =  0  which 

x  y  z 

is  one  of  Maxwell’s  four  equations.  Thus,  Eq.  (3-2)  is  a  solution 
of  Maxwell’s  Equations.  With  reference  to  Figure  3-1,  Eq.  (3-2)  re¬ 
presents  the  superposition  of  two  s-polarized  plane  waves,  one  tra¬ 
veling  in  the  direction  of  increasing  z  and  th«  others  traveling  in 
the  direction  of  decreasing  z.  Furthermore,  the  phasor  components, 

H  ,  H  and  H  of  the  magnetic  field  are  related  to  the  phasor 
x  y  z 

components,  E  ,  E  ,  and  E  of  the  electric  field  by  one  o* 
x  y  z 

Maxwell's  Equations.  The  relationship  is  given  below 


H  =  j.  Ph . 

x  wy  y  ay  az/ 

(3-5a) 

•  /af  aE  \ 

H  _  _L  ( _Ji  2  ) 

y  ~  wy  y  az  ax/ 

( 3-5h ) 

z  uu  \  ax  ay/ 

(3-Sc ) 

,i 

i  ■ 

Therefore,  the  magnetic  field  corresponding  to 
given  by  Eqs.  (3-2),  (3-3),  and  (3-4)  is 

the  elect. r ic  fiel  d 

t. 

1 

i 

H  -  H  +  u 
x  xR  xL 

( 3-6a ) 

|i 

i _ : _ j 


r 


where 
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H 

y 

=  H  _  + 

yR 

V 

( 3— 6b ) 

H 

z 

=  HzR  + 

HzL 

( 3—6  c ) 

’xR  = 

0 

(3-7  a) 

yR  = 

v?cos 

e  FxR 

(3- 7b) 

*zR  = 

■Ji s1n 

6  ExR 

( 3—7  c ) 

'xL  - 

0 

(3-8a) 

cos  e  E 


xL 


sin  e  E 


xL 


( 3-8b ) 
(3-Sc ) 


Similarly  Maxwell's  phasor  wave  equation  for  the  magnetic  Meld 
in  a  homogeneous,  isotropic  dielectric  is  given  Dy  Eg.  (3-9)  below 


3?H 

X 

+ 

2  7 

3  H  3'H 

*  + 

3 

2 

3  X 

ay 

3Z 

2,, 

2m 

2m 

3  H 

3  'H 

3  H 

o 

+ 

l*  - 

3  X 

ay 

3Z 

2m 

2m 

2m 

3  H 

3H, 

a  H 

£. 

+ 

_ _  4-  _ 

3X? 

ay2 

3Z 

4  +  K2H  =0 
2  x 


|  +  K?Hz  =  0 


( 3— 9  a ) 


( 3-Qb ) 


( 3- 9c ) 


where  K  =  as  before.  A  solution  of  the  wave  equation  is 


xR  +  HxL 

Mo 

1 

O 
C U 

yR  +  Hyl 

o 

1 

n 
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where 


Hz  =  HzR 


+  H 


zL 


(3- 10c) 


,  ,JK((sin  e)y  -  (cos  e)z) 

'xR  =  ae 

(3-1  la ) 

'yR  ’  0 

(3-1  lb) 

'zR  -° 

(3-llc) 

I  _  b  iK((sin  e)y  +  (cos  e)z) 
xL 

( 3— 12a) 

',i-° 

(3- 13b) 

'zL  '  0 

( 3— 12c) 

with  a,  b,  and  o  arbitrary  constant?.  With  reference  to  Figure  1, 
Eq.  (3-10)  represents  the  superposition  of  two  p-polarized  plane 
waves,  one  travelino  in  the  direction  of  increasing  z  and  the  other 
traveling  in  the  direction  of  decreasing  z.  A  solution  of  the  wave 
equation  (3-9)  is  a  solution  of  Maxwell's  Equations  if  in  addition 
it  also  satisfies 


V  •  (Hxx  +  H^y  +  Hzzl  =  0 


which  is  one  of  Maxwell's  four  eciuations.  Furthermore,  the  phasor 

components,  E  ,  E  ,  and  E  ,  of  the  electric  field  are  related 

x  y  z 

to  the  phasor  components,  H  ,  H  ,  and  H  of  the  maqnetic 

x  y  z 

field  by  one  of  Maxwell's  equations.  The  relationship  is  given  below 


E 


x 


E 


y 


( 3— 13a) 


(3- 13b) 
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Therefore,  the  electric  field  corresponding  to  the  magnetic  field 
given  by  Eqs.  (3-10),  (3-11),  and  (3-12)  is 


where 


F  -  F  +  F 

a  ~  xR  lxL 

(3- 14  a) 

Ey  =  EyR  +  EyL 

(3- 14b) 

Ez  -  EzR  +  ezl 

(3-1 4c) 

ExR  =  ° 

( 3—  1 5  a ) 

Ey«  'l/f  cose  HxR 

( 3— 15b) 

EzR  -  -yE  Sint,  HxR 

( 3- 15c 1 

ExL  ■  0 

(3- 16a) 

EyL  .  Vf  HxL 

( 3- 16h ) 

EzL  -  -Vr  Sinn  Hxl 

(3- 16c) 

Eqs.  (3-2),  (3-3),  (3-4),  (3-6),  (3-7),  and  (3-8)  constitute  an 

s-polarized  two  plane  wave  solution  of  Maxwell's  Equations  within 

each  layer  of  the  dielectric  stack  shove  in  Figure  3-1,  while  Eqs. 
(3-10),  (3-11),  (3-12)  (3-14),  (3-15),  and  (3-16)  constitute  a 

p-polarized  two  plane  wave  solution  of  Maxwell's  Fquat.ions  within 

each  layer  of  the  dielectric  stack.  In  order  for  the  corresponding 
sets  of  equations  to  he  valid  solutions  fnr  the  elect^omaoi  eti c 
fields  existing  within  the  entire  stack.  Maxwell's  boundary  condi¬ 
tions  must  be  satisfied  at  the  interfaces  between  the  layers.  We 
wish  to  show  that  Maxwell's  boundary  conditions  can  be  satisfied 
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by  tliese  equations  by  choosing  appropriate  values  for  a,  b,  and  e  in 
each  layer.  Consider  the  i^1'  boundary  and  translate  the  coordi¬ 
nate  system  so  that  z  =  0  at  this  boundary.  Maxwell's  six  boundary 
conditions,  i.e.,  continuity  of  the  tangential  components  of  E  aria  H 
and  continuity  of  the  normal  components  of  the  D  and  B  are  giver, 
be i ow : 


E(j)  =  E(J)' 
x  x 


c(j)  E(J)' 

y  y 

H(j)  =  H(jJ' 


H(j)  _  H( j ) ' 

y  y 


e  e  (j)  -  . 
c  /  •  \  l  _  -  e 


(j)  cz  =  e(j+i rz 

j(j)  u(j)’ 


r(j) 


( 3— 1 7  a ) 
( 3— 17b) 
( 3— 17c) 
( 3-1 7d ) 
(3-17e  1 
( 3—  1 7  f  ) 


*th 


where  superscript  (j)  denotes  at  the  j  boundary  on  tne  (j-1) 
layer  side 

superscript  (j)'  denotes  at  the  jth  boundary  on  the  jtn 
layer  side 

t 

subscript  (j)  denotes  in  the  j  layer 

Consider  first  the  s-polarized  case.  Eqs.  ( 3- 17b),  ( 3— 17c),  and 
(3-17e)  are  trivially  satisfied  for  arbitrary  a^_^, 

6(j  1)’  a(j)’  ^(j)  and  Boundary  condition  ( 3— 17a)  can 

be  explicitly  written  as 


iK,  .  n,sine,.  .  >y  ,  .  iK,  .  ,>sine,  ,,y  iK,.,< 

](j-l)6  ^  U"1)  +  (J"1)  (J-')  =  (J) 


sine  v 


’(J-U 


(j)' 

-  b(j)eiK  (j)sinw(j)y 


\  J 


1-  In 
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Equation  (3-18)  implies 

or  equivalently 

n(j-DS1nS(j-l)  =  n(j)sine(j)  ( 3-20) 

and 

au-i>  *  b(j-D  -  a(j) +  b»)  l3-21) 

Equation  (3-20)  will  be  recognized  as  Snell's  Law  of  Refraction. 
Boundary  condition  (3-17d)  can  be  explicitly  written  as 


7^°swuVK<3-'>sine<3-'>y 

eiK0)sineO)y 

C0M  (J)bcj>e<K(J)Sine(J)y  (3-2?) 

Equation  (3-22)  implies,  Eq.  (3-19)  above  and  also 
n(j-l)cose(j-l)(a(j_l)  “  b(j-l ))  =  n(j)C0S8(j)(a(j)  "  b(j)) 

(3-23) 

By  explicitly  writing  out  boundary  conditions  ( 3— 1 7 f ) ,  it  would  be 
seen  to  be  equivalent  to  boundary  condition  ( 3— 17a).  Solving  Eqs. 
(3-21)  and  (3-23)  simultaneously  yields 

[nmC0S9(j) +  "(j-if059  (j-D>u-D  *[nmcosem  - 

a(j)  ’  2n"(j)  co^i)  TT7 
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.[■ 


"u)cos#m  '  •  [\ 


J)C0**(J)*"!j-li£04*'j--JB 


i J )COi*( J ) 


Thus,  qiven  a(j_,).  h(j_i).  «(  j- 1 ) ,  Eqs  .  {  3-?V  , ,  '  3- :<  , 

and  (3-25)  determine  the  a^,  b^,  and  «^  )  necessary  ir 

der  to  satisfy  the  boundary  conditions. 

We  repeat  the  above  procedure  for  the  p-polari^ed  cases,  lqs. 
( 3— 17a),  (3-17d),  and  (3-17f)  are  trivially  satisfied  tor  jrnitra** 

a(j-l)’  h(j-D’  °( j-D’  "(j)’  b(j)’  and  HU).  Poundarv 

condition  (3- 17c)  can  be  explicitly  written  as 


i  iK ,  .  n^ins,;  ,  »y  .  .  i  K ,  ,  .sine,  ,.y 

d(j_l)e  (J-H  +  b(j-l)e  (J_1)  (J-1) 


.(j).iKU)Sin8[il’  *  bu)eiK(j)5lne(jl'' 


Eq.  (3-25)  implies 


K( j-1 )Sine( j-1 )  =  K(j)Sine(j) 


or  equivalently 


n(j-l)S1n*(j-l)  =  "(j)sine(j) 


3- H- 


and 


a(j-D  +  b(j-D  =  3(j)  +  b(j) 


3-  'u 


Boundary  condition  ( 3— 1 7 b >  can  be  explicitly  written  as 


C0S9(j-i)a(j-i)e 


iK(j-ns1n®(j-i  )y 


V  (j-i ) 


O  „  » K  ,  $  ’ • " 


-p- 

V  *  ( J ) 


cos® (j )e'K(d )s1n*(J )* 


(Equation  continued  on  next  pace) 
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„  /  uo  .  lK  ,  . .  s  me,  . .  y 

7Tm  ,J|  lj>e 

Eu .  '’-’t'1'-  implies,  Fq.  <: \-'‘7)  above  and  akn 

’•■'viihj-n  ■  "(Mi]  ■ 

Rv  explicitly  writmq  out  boundary  condition  (3-17e),  it  would  be 
seen  to  be  equivalent  to  boundary  condition  ( 3- 1 7c ) .  Solving  equa- 
1 1  ws  (  ?-■»«)  and  (3_31)  simultaneously  yields 


.--?S8  J 


1C0S8 


:ose, 


,„nC0S9- 

jsLL  u 


COSO  ,  ^ 


n .  .  .COSO.  .  ,3 
4-  _ L>_ 


Fnus,  q  l  v  en  i.j_]),  t>(._  ))t  and  Eos.  (3-?R), 

and  :  i.  a  detmnnm  ■>  the  a,.f,  t  and  neressa-'y 

der  to  satisfy  the  boundary  conditions. 


m  or- 


we  -live  a  i  v  en  a  solution  of  Maxwell's  Equations  and  toe  i.sn- 
riate<t  boundary  conditions  m  the  multilayer  dielectric  stac  and 
trie  idmin'iiij  free  spare  regions  shown  in  Figure  <- 1 .  Tne  melds  m 
'•icl'  layer  anl  the  adjoining  free  space  regions  can  tie  described  t  y 
t.ue  superposition  of  two  \  or  t.xO  p  polarized  plane  wa\>'<  .  These 
plane  waves  propaqate  m  ftic-  v-z  plane  and  have  directional  cosmos 
wb  irh  are  the  s  I'ne  in  the  f  *  irer  t  ton  and  "dual  in  maqn  ltud.  f  .*  p- 
p.rite  m  siqn  i-i  tin'  z  dire-tion.  T»ie  amplitudes  o<  these  p’,,a>... 
wave1  unit  tneir  iarer  t  ».--n,a  1  e  r  .  .  *  i  1  1  v  o  v  i  ■  .  ,m  1  .a a •  t  '  ’  -vr. 
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2e«im 


RADAR  AND  OPTICS  DIVISION 


Define  Y^  by  Eq.  (3-34) 

l  H 


Y  /  • 


YR 


(j) 


-XR 


(j) 


'XR 


(j) 


^  ™(j) 


>  for  s  polarization 


(3-3*1 


>  for  p  polarization 


Then  by  Eqs.  ( 3— 7b )  and  ( 3— 15b) 

'  n(j)C°S(j) 


Y,, 


V^O 


»  for  s  polarization 


(j)  =S 


'"(j) 


(3-35) 


(C0S®(j)^e0 


j=  .  for  p  polarization 


We  also  note  from  Eqs.  (3-8b)  and  ( 3- 1 6b )  that 
/ 


H 


YL 


(j) 


Y(j)  M 


'XL(j) 


XL 


(j) 


’YL(j) 


,  for  s  polarization 


,  for  p  polarization 


( 3-36 1 


Next  define 


vm 


1  1 


Ly(j)  ~y(j) 


(3-37) 
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by  Eqs.  ( 3—  3a )  and  ( 3—4 a ) ,  while  for  the  p-polar ization  case  we  can 
wr  i  to 


by  Eqs.  (3-1  la). 


or  equivalently 


'p(j),n 

rR 

e^(j)  0 

c(j)' 
LtYL  J 

_0  e~i0(j ) 

(3-12a),  (3-15b)  and  ( 3— 16b).  In  both  cases 

*(j)  =  K(J)  dj  C°S  9(j) 

*(j)  =  ^  "(3)  d(j)  C°S  °(j) 


(3-42) 


(3-43) 


(3-44) 


where  x  is  the  free  space  wavelength  corresponding  to  frequency  >j. 
Combining  Eqs.  (3-39)  and  (3-41)  yields  in  the  case  of  s-polarization 


where 


'c(j)' 

r(j  +  l  f 

tXR 

bXR 

- '(j-D  ,(i)u(j) 

r(j) 

F(J  +  1) 

.XL  . 

tXL 

— 

ei0m 

0 


0 


(3-*E) 


(3-46) 


Combinina  Eqs.  (3-40)  and  (3-42)  yields  in  the  case  of  p-polar ization 


c(j) 

r'YR 

tYR 

=  V(J-D  V(j)  U(j) 

r(j) 

r(j+l  ) 

L  'L  J 

Llyl  . 
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Terim 


Usinq  Eqs.  (3-45)  and  (3-47)  repeatedly,  we  can  write 


r  F(0)' 

|  EXR 

jF(°) 

L  xi.  J 


=  V(0)  (V(1  )U0  )V0  ))(V(2)U(2)V(2): 


*  *  •  (V(N-1)U(N-1)V(N-1))(V(N)U(N)) 


fF(N+ir 


v(N+l) 

L*xl  . 


for  s  polarization  (3-48) 


=  V(0)(V(1)U(1)V(0))(V(2)U(?)V(?)) 


*  *  (V(N-1)U(N-1)V(N-1))(V(N)U(N)) 


:(N+1)' 


.  ( N+l  ) 


for  p  polarization  (3-49) 

flow  with  the  aid  of  Eq.  (3-40),  Eqs.  (3-48)  and  (3-49)  become 

(VO)! 


/(C)(V(l)U(l)V(i))(V(2)U(2)V(?)) 


F ( N+1  ) ' 

(V  U  v_1  )V  ^ 

•  lV(N)U(N)V(N);V(N+l)  (N+1), 

.  EXL 


for  s  polarization  (3-50) 
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r(0)' 

rR 

F{0) 

L  VL  J 


v(o)  {vn)u(i)v(i))(v(?)u(2)v(2)) 


(V(N)U(N)VTn) 


(N+l) 


-  ( N+l ) '" 
YR 

-(N+l)1 


for  p  pol^riza' 


We  define 


Mm\  =  V/JXU,  .vVT1 


'(J) 


(j)U(j)V(j) 


N 


M  =  /7  M 

j=l 


(j) 


s  4  v7o)  M  V(N*1) 


tion  (3-51) 

(3-5 2) 

( 3-53 1 

(3-54) 


The  components  cc  the  2x?  matrices  S,  and  M  will  be  desig¬ 
nated  bv  S^,  ^ jK  and  MjK  respectively  (j,  K=l,  2). 

It  is  noted  that  Eqs.  (3-50)  and  (3-51)  relate  the  tangential 
electric  field  components  of  the  two  propaqating  plane  waves  in  the 
free  space  region  at  the  first  boundary,  to  the  tangential  electric 
field  components  of  the  two  propaoating  plane  waves  in  the  free  space 


written  as 


♦  1th 

boundary. 

Eqs.  (3-50) 

and  (3-51)  can 

c(0) 

^XR 

c  F(N+1)' 
=  ^irxR 

+  S  ) 

;>12tXL 

1 

for 

s  polarization 

O  _J 
^  X 

LlJ 

c  F(N+1)* 
-  znXR 

+  s  E  ^ 

22  XL 

1 

c(0) 

tYR 

-s  e(n+1)' 

"  HrYR 

+  s  F  ■  ' 

H2rYL 

for 

p  polarization 

F(0) 

LYL 

-s  e(n+1) 

-  a2ltyR 

+  <  F(N+1)' 
''?2tYL 
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If  a  plane  wave  is  incident  on  layer  No.  1,  then  it  is  easy  to  see  that 

4r+,)'  t\?'r 

Eq.  (3-55)  yields 


E^"'  and  E„i'-,/  must  be  identically  zero.  Thus,  for  s  polarization 


F(0) 

fcXL 

TToT 

ar 


J?1 

s 


11 


(3-57) 


and 


-(N+1 )' 
:XR 

TTor 

tXR 


’ii 


(3-58) 


Similarly,  for  p  polarization 


E 


(0) 

YL 

JW) 

CYR 


2ji 

>11 


(3-59) 


and 


We  define 


-(N+1 )' 
-YR 

TTOT 

tR 


’ll 


(3-60) 


:  (0 ) 
’XL 

7W' 

lxr 


:(0) 

’YL 


■Tor 

’YR 


for  ^-polarization 


for  p-polarization 


(3-61) 
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lF(N+1)' 

lXR 


7(or 

-XR 


,for  s  polarization 


1  H 


:  3-6? ) 


-  ( N+ 1  ) 1 
YR 


cos  e 


0 


■(0)  cos  e. 


-.for  p  polarization 


-YR  (N+l) 

It  will  be  noted  that  r,  which  we'll  call  the  reflection  coefficient, 
is  the  ratio  of  the  electric  field  of  the  plane  wave  reflected  from 
the  multilayer  stack  to  the  electric  field  of  the  plane  wave  incident 
upon  the  first  layer  of  the  stack,  t,  which  we'll  call  the  trans¬ 
mission  coefficient,  is  the  ratio  of  the  electric  field  of  the  plane 
wave  which  leaves  the  last  layer  of  the  stack  to  the  electric  field 
of  the  plane  wave  incident  upon  the  first  layer  of  the  stack.  Both 
r  and  t  are  complex  quantities,  each  having  a  magnitude  and  a  phase. 
This  is  noted  explicitly  bv  writing 


|r|e 


Jy. 


(3-63) 


t  =  |t(  ejYt 


Usinq  Povntinq's  Theorem,  the  power  density  across  the  x-y  plane 
due  to  the  incident,  reflected  and  transmitted  plane  waves  is  given 
bv  Eqs.  (3-64),  (3-65),  and  13-66)  respectively 

1 1 1 1  iiiii 

>  for  s  polarization 


ic  (0)  H(0) 

CXR  '  YR 


I 


(3-6«) 


F(0)  H(0) 

lYR  *  XR 


P  = 
p 


»  for  p  polarization 
f  I  “XL ^  ’  <1  ’  for  s  Polarization 

ft0’-  HXL  I  ’  ^or  P  polarization 


(3-65) 
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p ( N+1 ) 1  h(H+])'\ 


'XR 


'YR 


>  for  s  polarization 


(3-66) 


( |\|+ 1  ) 1  (N+l  ) 1 

(jrYR  ‘  HXR  ’  *r°r  P  polarization 


Eq.  f  3—64 )  is  rewritten  '»•  i  t  h  the  aid  of  Fqs.  (3-7b)  arm  (3-166)  as 

/ 


ho) 

— —  cose 

r(0) 

1  uo 

XR 

PT  -< 


(0) 


uo  coseo 


:  (0 ) 
YR 


for  s  polarization 


for  p  polarization 


(3-67) 


Eq.  (3-65)  is  rewritten  with  the  aid  of  Fqs.  (3-Sb)  and  (3-166)  as 

/  , - 


-^-^cose 

0 


,  for  s  polarization 


< 


(3-68) 


:  (0 ) 
YL 


cose 


for  p  polarization 


Eq.  (t_66)  is  rewritten  with  the  aid  of  ros.  (3-76)  and  ( 3-1 5b )  as 

( 


cose,.  ,  v 
u0  (N+l) 


.(N+l  )’ 
'XR 


,  for  s  polarization 


PT  =< 


-6«1 


ft? 

uo 


-N+l )’ 


1  )  1  YR 


cose , 


,  for  p  polarization 


'(N+l) 

Oefininq  the  reflectance,  ,  by 


APo 

P=^ 


-70) 
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and  the  transmittance,  T,  bv 


T  *!i 
PI 


Eqs.  (3-67),  (3-68),  and  (3-69)  yield 


F(0) 

lxl 

TToy 

lxr 

r  (0) 

-YL 

7m 

-YR 


,  for  s  polarization 


.  for  p  polarization 


(3-71) 


( 3-7? 1 


/n(N+l)  C0S9(N+1) 


'(o) 


cos  e 


0 


-  (N+l  ) 1 
-XR 


mr 

‘XR 


"(N+l)  C0S  60 


n  ( o )  C0S6(N+1) 


.(N+l) 

-YR 


mr 

YR 


»  for  s  polarization 

(3-73) 

,  for  p  polarization 


Combining  Eqs.  (3-61),  (3-62),  (3-7?),  and  (3-73),  we  have 


=  M 


(3-74) 


~ cos' ''IT 1  )>tl  ’  for  s  P°1arization 
n(ol  cos  90 


T  < 


(3-7*-' 


"(N+l )  C0S  9 (N+l)  |  t| 


^  "  (o)  cos  60 


,  for  p  polarization 


Thus,  we  have  derived  the  reflectance,  phase  chanqe  upon  reflection, 
transmittance  and  phase  change  upon  transmission  of  the  multilayer 
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dielectric  stack  shown  in  Fiaure  3-1  when  illuminated  by  an  s  or  p 
polarized  plane  wave  of  arbitrary  frequency  at  an  arbitrary  angle  of 
incidence.  Any  plane  wave  incident  upon  the  stack  can  be  written  as 
the  sum  of  an  s-polarized  and  a  p-polarized  plane  wave.  Due  to  the 
linearitv  of  Maxwell's  Equations,  each  of  these  plane  wave-  can  be 
treated  separately  using  the  above  procedure. 

Below,  we  give  an  algorithm  for  computing  the  reflectance,  phase 
change  upon  reflection,  transmittance  and  phase  change  upon  trans¬ 
mission  of  fhe  multilayer  dielectric  stack  shown  in  Figure  ?— 1  when 
illuminated  by  an  s  or  p  polarized  plane  wave  of  arbitrary  frequency 
at  an  arbitrary  angle  of  incidence  is  detailed. 


1.  An  s-polarized  or  p-polarized  plane  wave  having  a  wavelength 
x  in  layer  No.  0  is  incident  at  an  anqle  eQ  upon  the 
multilayer  dielectric  stack 


3. 


Fa.  (3-20),  which  is  a  statement  of  Snell's  Law  of  Refrac¬ 
tion,  is  used  to  compute  ©^  for  each  layer  of  the  stack. 


Fa.  (3-35)  is  used  to  compute  for  each  layer  of  trie 

stack 


A.  Fqs .  (3-371  and  (3-38)  a^e  used  to  compute  V,.,  ana 

_  i  v  J  J 

V(j)  for  each  laypr  of  ‘he  stack 

8.  Fqs.  (3-44)  and  (3-46)  are  used  to  compute  and  li^ 

for  each  layer  of  the  stack 


6. 


Eq  (3-5?)  is  used  to  compute  for  e*ch  layer  of  tun 

stack 


7.  Eq.  (3-53)  is  used  to  compute  M 

8.  Eq.  (3-54)  is  used  to  compute  S 

9.  r  and  t.  are  computed  b'-  Eqs.  (3-5Q)  through  H-fi?' 

10.  p  and  T  are  computed  by  Eqs.  (3-74)  and  (3-75) 


62 


RADAR  AND  OPTICS  DIVISION 


It  is  customary  to  ref  or  to  Y 
layer,  V 


(j) 


(j)  as 

as  the  admittance  matrix  of  the  j 


the  admittance  of  the  j 
th 


■  to 


layer,  (J 


.th 


(j) 


as 


the  phase  transformation  matrix  of  the  i  '  layer,  M  as  the  inter¬ 
ference  matrix  of  the  entire  stack  and  S  as  the  stack  system  transfer 
matrix. 


Later  on,  we  will  use  a  number  of  results  contained  in  t*'is  chap¬ 
ter.  It  would  be  more  convenient  at  that  time  if  these  results  were 
in  a  somewhat  different  form  that  which  they  appear  in  here.  Antici¬ 
pating  this,  we  will  present  the  alternate  forms  now. 


Consider  a  single  layer  stack.  Then  fqs.  (3-50)  and  (3-51) 
become 


>)■ 

~r(2)r 

tXR 

=  V(0)V(1)U(1)V(1)V(2) 

CXR 

p(0) 

r(2)' 

Llxl  J 

L  XL  j 

rF(2)’l 

lyr 

=  V(0)V(1)U(1)V(1)V(2) 

lyr 

c(0) 

c(  2)' 

L  yl  J 

_eyl  J 

,  for  s  polarization 

( 3-70) 


,  for  p  polarization 

(3-77) 


The  two  by  two  system  transfer  matrix  of  this  single  layer  stack  will 
be  designated  L.  Then 


L 


■  v(o)vn)u(i),n)v(n 


(3-78) 


Equations  (3-43)  and  (3-46)  indicate  that  approaches  the  identity 
matrix  in  the  limit  as  d^ ,  the  layer  thickness,  approaches  zero.  Thus 

G  =  lim  L  =  VTLv,-.  (3-79) 

d] *0  n 

The  theory  developed  in  this  chapter  was  used  to  determine  ‘he 
reflectance  versus  wavelength  characteristic  of  a  twelve  layer 
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"high-low,  high-low"  multilayer  dielectric  stack.  The  refractive 
index  profile  and  corresponding  reflectance  characteristic  are  shown 
in  Figures  3-2  and  3-3  respecti vely.  Similarly,  we  calculated  the 
reflectance  versus  wavelength  character istic  of  a  15  pm  thick  re¬ 
flection  hologram.  The  hologram  has  a  refractive  index  variation  in 
the  z-direction  (see  Figue  2-1)  given  by 

n(z)  =  1.5  +  0.01  sin  (2*  (5.285)  z) 

Furthermore,  the  hologram  is  assumed  to  be  surrounded  by  a  medium 
having  a  refractive  index  of  1.5.  The  solid  line  in  Figure  3-3  snows 
the  reflectance  versus  wavelength  characteristic  calculated  by  tne 
theory  of  stratifiea  media.  The  calculation  was  performed  by  divid¬ 
ing  the  hologram  into  7927  slabs  of  equal  thickness.  The  refractive 
index  within  each  slab  was  assumed  to  be  constant  and  equal  to  the 
refractive  index  of  that  slab  at  its  front  surface.  The  dotteu  line 
shows  the  reflectance  versus  wavelength  calculated  by  Kogelnik's 
coupleo  wave  theory.  Although  Kogelnik's  theory  is  not  exact, whereas 
the  theory  of  stratified  media  is,  Koqelnik 1 s  theory  is  seen  to  give 
reasonauly  good  results  in  this  particular  case. 
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Notes  to  Chapter  3 

[1]  F.  Abel es ,  Ann.  d.  Phys.,  5,  596  (1  950). 

[2]  Z.  Knittl,  Optics  of  Thin  Films,  John  Wiley  &  Sons,  New  York,  1976. 
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4 

THE  THEORY  OF  THICK  SLAB  DECOMPOSITION 


In  this  chapter  we  develop  a  technique  for  determining  the  frequency 
domain  transfer  functions,  HR(A,  0Q)  and  Hy(A,  0  ),  of  a  flat,  hologra¬ 
phic,  phase  reflection  grating.  As  before,  the  grating  is  assumed  to 
have  a  refractive  index  variation  only  in  the  direction  normal  to  its 
surface.  The  technique  will  be  referred  to  as  the  thick  slab  decompo¬ 
sition  (TSD)  method. 

Consider  the  situation  shown  in  Figure  4-1.  A  plane  wave  of  wavelength 
A  is  incident  upon  the  front  surface  of  a  flat,  holographic,  phase,  reflect¬ 
ing  grating  at  an  angle  e  .  Without  loss  of  generality  the  plane  of  inci¬ 
dence  is  taken  to  be  the  x-z  plane.  The  holographic  grating  is  assumed  to 
be  of  infinite  extent  in  the  x  and  y  directions,  and  it  has  a  refractive 
index  which  varies  only  as  a  function  of  z.  The  front  surface  of  the 
holographic  grating  is  adjacent  to  a  homogeneous  dielectric  of  refractive 
index  nf  while  the  back  surface  is  adjacent  to  a  homogeneous  dielectric 
of  refractive  index  n^.  The  complex  electric  field  amplitude  at  x  -  0, 
z  =  0"  of  the  incident  wave  is  designated  FR.  There  will  be  a  field  re¬ 
flected  from  the  holographic  grating  and  a  field  transmitted  through  the 
grating.  Both  of  these  fields  will  be  plane  waves.  The  complex  electric 
field  amplitude  at  x  =  0,  z  =  0  of  the  reflected  plane  wave  is  designated 
El  while  the  complex  electric  field  amplitude  at  7  =  0,  z  =  D+  of  the 
transmitted  plane  wave  is  designated  EJ  .  £  ,  E”  and  E'  are  complex  vec- 

K  K  L  R 

tor  quantities  having  x,  y  and  z  components.  The  right  going  reflection 
coefficient,  rR,  and  the  right  going  transmission  coefficient,  tR,  are 
defined  below  by  Eqs.  (4-la)  and  (4-lb),  respectively. 


AxL 

:xR 


rR  = 


EyR 


for  s  pol ari zation 


for  p  polarization 


(4-la) 
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Grating  Refractive  Index  =  n(z) 


Figure  4-1.  Holographic  Grating,  Definition  of  Reflection  and 
Transmission  Coefficients 
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-xl 

:xR 


for  s  polarization 


E'  cos  0 

■f—  •  or  .  for  p  polarization 
yL  cos  o 


(4-lb) 


It  immediately  follows  from  the  definitions  of  Chapter  2  that 


HR  <X-V  ■  rR 


(4-2a) 


and 


Ht  (X,0o^  =  tR 


(4-2b) 


If  the  grating  is  illuminated  from  the  reverse  side  analogous  left  going 
reflection  and  transmission  coefficient  (rL  and  tL)  can  be  defined. 
Furthermore,  knowledge  of  rR  and  rL  completely  determines  tR  and  tL 
This  result  follows  from  the  reciprocity  theorem  of  horizontally  strat¬ 
ified  media  [1  J.  Thus,  HR  ( x ,eo )  and  H-j.  (X,0Q)  can  be  determined  if 
rR  and  r^  are  known. 

The  right  going  reflection  coefficient  (and  consequently  the  left 
also)  can  be  rigorously  determined  using  Maxwell's  Equations.  To  be 
more  explicit,  let  us  suppose  that  the  section  between  z  =  0  and  z  =  zQ 
( zq  <  D)  is  removed  from  the  holographic  grating  shown  in  Figure  4-1. 

We  will  designate  the  right  going  reflection  coefficient  of  the  remaining 
grating  by  rR(zQ).  Then  Maxwell's  equations  yield  the  following  non¬ 
linear  differential  equation  and  associated  boundary  condition  for 
rR(z)  [  2  ]. 


drR(z) 

dz 


=  2iky(z)rR(z) 


dB(z)  _ 

dz  20  ( z 


[1  -  rR(z)] , 


0  <  z  <  D  (4-3a) 
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Rather  than  solving  Eq.  (4-3a)  directly,  an  alternate  solution 
is  to  divide  the  grating  shown  in  Figure  4-1  into  a  very  large  number 
of  thin  slabs  parallel  to  the  grating  surface.  The  slabs  are  chosen 
to  be  sufficiently  thin  so  that  the  refractive  index  across  each  thin 
slab  is  to  a  good  approximation  constant.  Thus,  each  thin  slab  may  be 
considered  to  be  a  homogeneous  dielectric.  The  holographic  reflection 
grating  is  being  modeled  as  a  multilayer  dielectric  stack.  In  a  straight¬ 
forward  and  rigorous  manner  the  reflection  coefficient  of  such  a  stack 
may  be  determined  (using  the  theory  of  stratified  media  given  in  Chap¬ 
ter  3).  The  approach  involves  determining  the  two  by  two  interference 
matrix  associated  with  each  slab  and  then  multiplying  together  the 
inteference  matrices  for  all  the  slabs.  It  can  be  shown  that  this  mul¬ 
tilayer  dielectric  stack  approach  is,  in  effect,  a  numerical  technique  for 
solving  Eq.  ( 4-3a )  [2]. 

The  approach  discussed  above  has  three  principal  drawbacks.  First, 
it  does  not  result  in  a  compact  solution;  second,  it  does  not  yield  in¬ 
sight  into  the  problem;  and  third,  it  requires  extensive  "number  crunch¬ 
ing."  Consider  a  15  Mm  thick,  holographic,  reflection  grating  having 
a  sinusoidal  refractive  index  variation  with  a  period  of  0.1715  Mm. 

Such  a  hologram  can  be  recorded  in  dichromated  gelatin  from  the  inter¬ 
ference  pattern  of  two  normally  incident  plane  waves  of  wavelength  0.5145 
urn.  The  recording  geometry  is  shown  in  Figure  4-2.  The  refractive  index 
of  dichromated  gelatin  is  nominally  1,5.  We  wish  to  determine  the  reflec¬ 
tion  coefficient  of  this  element  by  using  the  multilayer  decomposition  ap¬ 
proach.  Into  how  many  thin  slabs  should  each  period  of  the  refractive  index 
profile  be  divided  in  order  to  get  a  good  approximation  of  the  reflection 
coefficient?  The  answer  to  this  question  will  of  course  depend  upon 
the  magnitude  of  the  refractive  index  modulation.  Typically,  this  might 
be  on  the  order  of  0.02.  It  would  probably  be  unwise  to  choose  less 
than  ten  slabs/period  and  fifty  slabs/period  seems  more  reasonable. 

Ten  slabs/period  corresponds  to  875  matrix  multiplications  for  each 
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X  -  0.514um 


Figure  4-2. 


15um  l~«- 


-  0.514um 


Holographic  Reflection  Grating  Recording  Geometry 
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value  of  0Q  and  \  while  forty  slabs/period  corresponds  to  4375  matrix 
multiplications  for  each  value  of  dQ  and  X.  This  is  a  rather  large 
number  of  calculations. 

For  some  refractive  index  profiles  one  may  have  to  resort  to  the 
multilayer  decomposition  approach  in  order  to  determine  the  reflection 
coefficient.  We,  however,  are  interested  in  refractive  index  profiles 
which  are  holographically  produced.  Let  n(z)  denote  the  refractive 
index  profile.  n(z)  is  assumed  to  be  a  continuous  function  of  z  and 
its  Fourier  transform  will  be  designted  N(fz).  We  restrict  our  analy¬ 
sis  to  refractive  index  profiles  which  satisfy  the  conditions  below. 

|n(z)  -  nQ|  -  0.03  where  nQ  is  nominally  1.5  (4-8a) 


and 


N(f  )  =  0  for  4  pm"1  S  <  6  pm'1  (4-8b) 

These  conditions  are  not  artificial,  they  are  reasonably  consistent  with 
currently  available  holographic  recording  materials  and  recording  tech¬ 
niques.  We  will  elaborate  on  this  point  in  Chapter  8. 

The  technique  we  have  developed  for  determining  the  reflection  co¬ 
efficient  of  a  flat,  holographic,  reflection  grating  will  be  referred 
to  as  the  thick  slab  decomposition  (TSD)  method.  In  principle  its  quite 
simple.  The  holographic  grating  is  divided  into  a  relatively  small 
number  of  slabs  each  parallel  -to  the  surface  of  the  holographic  grating. 
The  slabs,  however,  are  chosen  to  be  sufficiently  thin  so  that  the  re¬ 
flection  coefficient  of  each  slab  by  itself,  when  adjacent  to  dielectric 
mediums  which  index  match  the  slab  at  its  front  and  rear  surfaces,  is 
small.  Next,  it  is  shown  that  for  our  case  of  interest  y (z)  in  Eq.  (4-3a) 
can  to  a  good  approximation  be  considered  constant.  The  reflection  co¬ 
efficient  of  each  slab  is  now  calculated  using  Eq.  (4-3a)  and  the 
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p 

approximation  that  1  -  rR(z)  r  1 .  It  should  be  noted  that  this  approxi¬ 
mation  together  with  the  approximation  that  y[z)  is  constant  reduces 
Eq.  (4-3a)  to  a  first  order,  linear  differential  equation  which  can 
be  solved  in  closed  form  by  integration. 

With  the  reflection  coefficient  of  each  slab  in  hand,  the  recipro¬ 
city  theorem  of  stratified  media  is  applied  in  order  to  determine  the 
two  by  two  transfer  matrix  of  each  slab.  The  matrices  for  all  the  slabs 
are  then  multiplied  together  to  yield  the  transfer  matrix  of  the  entire 
grating.  The  reflection  coefficient  of  the  grating,  as  well  as  HR  (A, 

0Q)  and  H-j.  (A,  0Q),  are  easily  determined  from  this  matrix.  The  tech¬ 
nique  is  described  in  detail  below. 

We  begin  by  considering  a  slab  between  z  =  a  and  z  =  b  of  the  grat¬ 
ing  shown  in  Figure  4-1.  Our  immediate  goal  is  to  derive  the  right  going 
reflection  of  this  slab.  The  slab  thickness  A 

A  -  b  -  a 

is  assumed  to  be  sufficiently  small  so  that  the  magnitude  of  the  right  go¬ 
ing  reflection  coefficient  of  this  slab  is  much  less  than  one  when  it  is 
adjacent  to  dielectric  mediums  which  match  the  slab  in  refractive  index 
at  its  front  and  back  surfaces.  The  discussion  which  follows  will  be  re¬ 
stricted  to  this  slab  under  such  conditions.  The  coordinate  system  for 
the  slab  will  be  as  shown  in  Figure  4-3.  The  slabs  refractive  index  pro¬ 
file  will  be  designated  n(z). 

With  the  approximation 

1  -  rR(z)  3  1  0  -  z  -  A  (4-9) 

Equations  (4- 3a)  and  ( 4 - 3 b )  become 
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drRlz) 
d  z 


-  2iky(z)rR(z)  +  ^  23(z)  ,  0 


<  z  <  A 


rR(A)  =  0 


This  equation  can  also  be  written  as 


(4-1 Oa ) 
(4-1  Ob) 


_d_ 

dz 


rR(z)  exp  (■ 


z  -|  z 

-21k  /  Y(C)  dC)  =  exp  (-21k  /  Y(C 

•'  +  •'  + 

nT  1 


C)  dO 

(4-11 ) 


Integrating  both  sides  of  Eq.  (4-11)  yields 


A 

/ 

d 

dz  1 

0+ 

rD(z )  exp  (-2ik  J  y(C)  d£) 
0 


z 

( 


dz  = 


A 

f  iFwyexp  {-2ik 

„+  L 


z 

/  y(C)  dfjj 


dz 


(4-12) 


A 

rR( a)  exp  ( -2 i k  J  y (0  40  -  rR(0+)  exP  ('2ik  J  y (0  d<;)  = 


Z 


exp 


(-21  k 


y(C)  d.O  dz 


(4-13) 
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The  boundary  condition  rR(A)  =  0,  Eq.  (4-10b),  is  now  applied  and 
Eq.  (4-13)  becomes 


rR<0 


*)  ■-  j  [^26TTT«P  <-2ik  j 


yU)  <k) 


dz  (4-14) 


0  0 
From  Eq.  (4-7),  i.e.,  rR  =  rR(0+),  one  gets 


rR  = 


j  ^2BWexp  (-21k  f 


yU)  *0 


dz 


or  equivalently 


rR  = 


ZJ 

[ [Mil  |  exD  r 
J  [  dz  MU  exp  1 


•2,k/ 


yU)  dc) 


dz 


0 


0 


(4-15) 


If  we  define  e(z)  by  Eq.  (4-16)  below 

n(0-)sin  e 


in  e(z)  = 


sin 


n(z) 


0 


(4-16) 


Then 


cos  e(z)  = 


2(0-)sin2  Q, 


n2(z) 


(4-17; 
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Using  Eq.  (4-17)  we  rewrite  Eqs.  (4-4),  (4-5a),  and  ( 4- 5 b )  as  given 
below. 


y(z)  =  n(z)  cos  0(z) 


li(z)  = 


0 

n(z)  cos  0(z) 
n(z) 

zn  cos  0(z)  ’ 


for  s  polarized  illumination 


for  p  polarized  illumination 


(4-18) 

(4-1 9a) 

(4-1 9b) 


We  wish  to  evaluate  /  3 ( z ) .  Consider  first  the  s  polarized  case. 


Then  from  Eq.  (4-1 9a) 


d8(z) 

dz 


[ dn(z) 


'  P(z)  =  ^Tr^TT(I)  [“dT^  cos  0(z)  -  n(z)  sin  0(z) 


dr(z) 


dz 


(4- 


20) 


Differentiating  Eq.  (4-16)  yields 

sin  0(z)  +  n(z)  cos  0(z)  d^Z^  =  0  (4-21) 


or  equivalently 


du(z)  _-dn (z)  tan  0( z) 
dz  '  dz  n ( z ) 


Combining  Eqs.  (4-20) and  (4-21) 


dB(z) 

dz 


/  EJ(z)  = 


d»(z). 

dz 


/  n(z) 


[1  +  tanZ  y ( z ) ] 


dn( z)  1 _ 

dz  n(z)  cos^  0(z) 


dn(z) 

dz 


n^  cos 


0 


(4-22) 


(4-23) 
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The  last  line  follows  from  the  fact  that  n(z)  -  nQ  (see  Eq.  ( 4-8a ) )  and 
so  e(z)  -  6q  (Eq.  (4-16)).  For  the  p  polarized  case  Eq.  (4-1 9b )  yields 


ygf-t  8(2)  - 


cos  0{ z)  d^z)  +  n( 2)  sin  e( z ) 


n(z)  cos  0(z) 
Next,  Eqs.  (4-22)  and  (4-24)  are  combined  to  yield 


dg(z) 

dz 


/  S(z)  = 


^-/  n(z) 


[1  -  tan^  e(z)] 


(4-24) 


dn(z)  (1  -  tan  90) 
dz  n„ 


(4-25) 


The  last  line  follows  from  the  approximations  n(z)  -  nQ  and  e(z)  -  eQ. 
Eq.  (4-15)  now  becomes 


'"‘I 


"/ 


exp  [-2ik  /  y(0  d£]  dz 


0 


(4-26) 


where 


a  = 


n0  cos  60 


1  -  tan  9, 


,  for  s  polarization 


,  for  p  polarization 


( 4-2  7a ) 


( 4-2  7b) 


We  wish  to  show  that 
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Z 


-2  i  k 


/ 


YU)  dS  -  -2iknQ(cos 


0  4  z  S 


(4-28) 


Equation  (4-4)  is  rewritten  as  shown  below 


V(z)  =  f/[n0  +  An(z)]2  -  (nQ  +  An(0))2  sin2  0Q  (4-29) 

i  ( z )  =  ^n2  -  n2  sin2  0Q  +  2nQ(An(z)  -  An(0)  Sin2  6Q)  +  [An(z)]2  -  [ '  n(0)]2  sin2  •  ( 


(4-30) 


U) 


--  ^  -  4 


sin 


2nQ(‘.n(z) 


An(0)  sin^  eQ) 


(4-31) 


where  n(z)  has  been  expressed  as  n^  +  .\n(z).  Note  that  |_\n(z)j  -  0.03 
from  1q.(4-8a).  In  Eq.  (4-31),  the  second  order  terms  in  An(z)  and 
,\n(0)  have  been  neglected.  Equation  (4-31)  is  now  expanded  in  a  power 
series 


2.2  1/2  !  2n0(AnU|^n(0)  sin  Q) 

Y(z)  <n0  '  n0  sin  V  2  ',2  2  .2  J/2 

(nQ  -  nQ  sin  "Q) 


1 

8 


2nQ(An(z)  -  An(0)  sin2  Oq)  2 

~2  2  .  2  „  TlfT 

(n0  -  nQ  sin  uQ) 


+  .  .  . 


(4-32) 
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From  Eq.  (4-32)  it  can  be  seen  that 


Iy(z)  -  nn  cos  0 


<  1 2An(  z ) 


0  cos 


<  < 

0  -  z  -  A 


(4-33) 


and  so 


z 

/ 


2K  |  y(5)  dS  -  2K  I  nQ  cos  0Q  d£ 


/"< 


2K 


cos  ar 


/2in 


U)  d£ 


.  0  <  z 


(4-34) 


We  want  to  find  a  strong  upper  bound  on 


2K 


COS  0 


z 

/ 


2 An (C)  dt 


0  -  z  -  A 


We  have  been  unable  to  do  this,  but  hasten  to  point  out  that  it  seems 
likely  that  a  rather  weak  bound  is  given  by  Eq.  (4-35)  below. 


2K 


cos  0 


L. 

I 


2AnU)  d£ 


<  0.18 


_  <  <  A 

0  -  z  -  a 


(4-35) 


The  following  plausibility  argument  is  presented  to  support  this 
conjecture.  Let  us  write 

An(z)  =  A(z)  sin  [2uf(z)z  +  <f>(z)] 

Equation  ( 4- 8b )  i n di cates  that  An(z)  is  to  an  extent  narrow  band.  Arbi¬ 
trarily  pick  p  such  that 
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Since  An(z)  is  to  an  extent  narrowband,  we  would  expect  that  there  exist 
some  :z  such  that  A(z),  f(z)  and  <P(z)  are  reasonably  constant  in  the 
interval  p  <  z  <  sz  and  such  that  f(p)6z  is  some  integer.  Now  partition 
the  interval  from  [0,  z]  into  sub-intervals  [p  ,  p  +  (<$z)  ],  m  =  1,  ... 

L  where  p^  =  0.  A(z),  f(z)  and  <t>(z)  are  reasonably  constant  throughout 

each  sub-interval  and  f(pm)(i5z)  is  some  integer  for  all  m  except  m  =  L. 
For  m  =  L,  f(p^)(6z).  is  less  than  one.  It  is  noted  that 


J A(£)  sin  [2Ttf(^)C  +  4>(C)]  dC 
0 


L-l  Pm+(6z) 

■£/ 


m=1  p 
rm 


m 


A(£)  sin  [2vrf(C)C  +  <£(<=)]  dc 


PL+(d2)L 


A(U  sin  [27TfU)C  +  4>U)]  dC 


L-l  p  +(<5z) 
'm 


/•I 

m=l  p 

r  n 


A (Pm)  sin  [2Trf(pm)C  +  d>( Pm) ]dC 


/ 


PL+(  6z)l 


+  I  A(pl)  sin  [2iTf(pL)L  +  4>( P^)  ]  d£ 

pL 


-  A(p, ) 


/ 


Pl+(6z)l 


sin  [2-nf(pL)L  +  PL ) 3  d^ 


(note  that  J 

Pm 


P  +(‘‘>z) 
Km  v  m 


A(pm)  sin  [2tt f ( Pm )r  +  ;(pm)]dr  =  0) 
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|a(pl) 1 

-  TTf(pLJ 
Consequently 

z 

cos  e0  J  2An^  j  -  Tif ( pL)  cos  eQ 
0 

"  a" f"(pL)'  cos  eQ  lA(pi_)l 
0  -  Z  -  A 

From  Eqs .  (8a)  and  (8b) 

| A(pL) [  -  0.03  ,  f(pL)  -  4  Mm_] 

Restricting  A  to  the  visible  spectrum 

0.45  pm  -  A  -  0.65  ym 

and  noting  that  eQ  <  42°  [3],  we  have 
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cot\ /  “<«>  d« 

0 


<  8 

"  (0.45) (T)  cos  42° 


(0.03)  -  0.18 


In  any  case,  the  validity  of  Eq.  (4-28)  can  always  be  checked  by 
evaluating  the  integral  on  the  right  hand  side  of  Eq.  (4-34).  We  will 
assume  that  Eq.  (4-28)  is  valid.  Equation  (4-26)  then  becomes 

A 

rR  -  -ot  J  d-|]yZ-  exp  [( -2i knQ  cos  0Q) z]  dz  (4-36) 

0 

The  above  equation  is  the  desired  result. 

We  briefly  comment  on  Eq.  (4-36).  By  rewriting  Eq.  (4-36)  as  given 

below 


r 


R 


rect  (  )  exp  (~i2nfzz)  dz 


(4-37) 


where 


rect  (  ^  )  = 


0  -  z  -  A 


0,  elsewhere 


fz  =  n0  cos  G0 


It  is  seen  that  rR  is  the  Fourier  transform  of 


dn  (z) 
dz 


rect  (  j  ) 


(4-38) 


(4-39) 
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One  can  also  give  a  physical  interpretation  of  Eq.  (4-36).  Such  an  in¬ 
terpretation  goes  as  follows.  The  slab  can  be  thought  of  as  being  com¬ 
posed  of  a  very  large  number  of  parallel  sub-slabs  each  of  thickness  dz 
and  constant  refractive  index.  Now,  if  n(zg)  is  the  refractive  index 
of  the  sub-slab  whose  front  surface  lies  in  the  plane  z  =  Zg  then  to  a 
good  approximation 

n(zg  +  dz)  -  n(Zg)  dz 


will  be  the  refractive  index  of  the  sub-slab  whose  front  surface  lies 
in  the  plane  Zg  +  dz.  As  a  ray  having  an  amplitude  of  one  propagates 
into  the  slab,  some  of  it  is  reflectecd  at  each  of  the  interfaces  sep¬ 
arating  the  sub-slabs.  The  amplitude  of  the  reflected  ray  at  the  plane 
z  =  Zg  is  given  by  Fresnel's  formula  [4]. 

n(zn)  cos  [0 ( z« ) ]  -  n(zn  +  dz)  cos  [0(z  +  dz)] 

,  for  s  polarization 

( 4- 40a ) 


n(z0)  cos  0(zo] 


+  n(zn  +'  dz)  cos  [8(zn  +  dz)' 


n(Zg  +  dz)  cos  [0(Zg) ]  -  n(Zg)  cos  [0(Zg  +  dz)] 


for  p  polarization 

( 4-40b) 


It  can  be  shown  with  the  use  of  Snell's  Law  and  a  bit  of  algebra  and 
trigonometry  that  Eqs .  (40a  and  40b)  become 


in  the  limit  as  dz  -*•  0.  It  is  assumed  that  the  amplitude  of  the  ray 
does  not  attenuate  significantly  as  it  propagates  into  the  slab,  and 
that  multiple  reflections  at  the  sub-slab  interfaces  can  be  ignored. 
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Then  the  right  going  reflection  coefficient  of  the  slab  will  be  the  sum 
of  the  reflections  from  each  of  the  sub-slab  interfaces.  The  phase  asso¬ 
ciated  with  the  reflected  ray  from  the  sub-slab  interface  at  z  =  z^  is 
given  by  the  total  optical  path  length  that  ray  travels  in  going  to  the 
plane  z  =  zQ  and  then  returning  to  the  front  surface  of  the  slab.  This 
optical  path  length  is  approximately 

2KnQ(cos  6q)Zq 


Thus , 


j 

— 1 

\ 

v-Z 

/  dn ( z ) 

01 

j  dzjexp  [(-2iKnQ  cos  eQ ) 2 j ] 

d2j 

\ 

z  =  z  J 

J  -J 

which  approaches  Eq.  (4-36)  in  the  limit  as  dz.  -*•  0. 

We  have  succeeded  in  deriving  a  closed  form  expression  for  the 
right  going  reflection  coefficient  of  a  slab  section  of  the  grating. 
Below,  the  two  by  two  transfer  matrix  of  this  slab  is  derived  from  know¬ 
ledge  of  its  right  going  reflection  coefficient.  Figure  4-4  shows  the 
slab  illuminated  from  the  front  side  with  a  plane  wave.  Similarly, 
Figure  4-5  shows  the  slab  illuminated  from  the  back  side  with  a  plane 
wave.  Er,  El,  Er  and  E^  are  the  complex  electric  field  amplitudes. 

They  are  vector  quantities  having  x,  y  and  z  components.  In  Figures  4-4 
and  4-5,  0Q  is  related  to  0g  by  Snell's  Law 


eA  =  sin"  ( 


1  n ( 0 )  sin  Gg 


n  (a) 


) 


(4-41) 


But  n(0)  -  n(A)  by  Eq.  (4-8a)  and  so 


u0  '  °0 


(4-42) 
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Suppose  now  that  the  holographic  grating  is  illuminated  simultaneously 
from  both  sides.  The  incident  angles,  0q  and  0q,  on  the  front  and  back 
sides  respectively,  are  assumed  to  be  related  by  Snell's  Law  (Eq.  (4-41)). 
The  principle  of  superposition  applies.  The  situation  is  depicted  as 
shown  in  Figure  4-6.  ER  and  IL  will  be  uniquely  determined  by  ER  and  . 
The  relationship  can  be  written  in  matrix  form  as  given  below 


and 


rE  i 

xR 

E 

xL 

qll  q12 


q21  q22 


'xR 


‘xL 


,  for  s  polarization  (4-43a) 


-yR 


■yL 


qll  ql 2 

q21  q22  i 


yR 

E' 


,  for  p  polarization  (4-43b) 


The  Q-jj's  in  Eq.  (4-43b)  are  not  the  same  as  those  in  Eq.  (4-43a).  The 
two  by  two  matrix  Q 


Q  = 


qll  q12 


q21  q22 


(4-44) 


is  called  the  transfer  matrix. 

The  right  going  reflection  and  transmission  coefficients  rR  and 
tR,  respectively,  are  defined  as  before 
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,  for  s  polarization 


,  for  p  polarization 


t 


R 


EyR 


cos 


0 


=  0 


for  s  polarization 


for  p  polarization 


(4-45) 


(4-46) 


Similarly,  the  left  going  reflection  and  transmission  coefficients,  r L 
and  t^ ,  respectively,  are  defined  by 


r 


L 


xR 

xL 


,  for  s  polarization 


,  for  p  polari zation 


cos  0 


0 


COS  0Q 

er  =  0 


for  s  pol  ari  zation 


for  p  pol ari zation 


(4-47) 


(4-48) 
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We  once  again  invoke  the  fact  that  0q  -  9„ .  Thus, 

cos  0q 


cos  6 


f  -  1 
0 


With  the  use  of  Eq.  (4-43a)  and  (4-43b)  it  is  easy  to  show  that 


rR 


qjM 

qn 


,  for  s  or  p  polarization 


(4-49) 


PR  ’ 


Ml 


Ml 


cos  0 


0  -  1 


cos  0q  qn 


,  for  s  polarization 


,  for  p  polarization 


(4-50) 


r,  = 


ql  2 


L  q 


11 


,  for  s  or  p  polarization 


(4-51) 


( 


-q21qi2  _ 

~  *22 

"q21q12 


,  for  s  pol  ari  zation 
(4-52) 


Ml 


cos  0 q'  _q21q12 

+  q22  >  *  cos"e^  -q^-  '  q22  ’  for  P  polanZatl°n 


Solving  for  qn  ,  q12>  q21  and  q22  in  terms  of  rR,  tR,  rL  and  tL  one 
gets 


Ml 


i/tr 


,  for  s  polari zation 


(4-53) 


COS  0 


( 1 /to)  •  - -f  -  l/tD  ,  for  P  polarization 

K  COS  i' 
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'rL/tR 


,  for  s  pol ari zation 


(-rL/tR! 


COS 

cos  Ql 


(4-54) 


- 5/ 1 r  ,  for  p  polarization 


V‘r 


,  for  s  polarization 


cos  0r 


(4-55) 


(rW  •  coFsr  :  rRytR  •  for  p  PoTari zation 


tLtR  '  rRrL 


tLtR  '  rRrL  x  Cos  90  -  tLtR  '  rRrL 


,  for  s  polarization 


(4-56) 


r.r~nr  -  - f - -  ,  for  p  polarization 

cos  Og  cR 


It  will  be  shown  below  that  tR,  t^  and  rL  can  be  derived  from  know- 
leged  of  rR.  We  start  with  rR  given  by  Eq.  (4-36),  which  is  repeated  belt 


exp  [{-2iknQcos  Oq)z]  dz  (4-57) 


The  left  going  reflection  coefficient  can  be  written  as 


r r  -aj  d-(-Ad~  exp  [(-2iknQ  cos  0g)z]  dz 


where 
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(4-59a ) 


( 4- 5  9  b ) 

With  a  simple  change  of  variable,  Eq.  (4-58)  becomes 

rL  -  a'  exp  [ ( - 2 i k n Q  cos  Og)A]  /  exP  [(2ikng  cos  0g)zl  dz 

{  (4-60) 

Combining  Eqs.  (4-57)  and(4-60)  and  noting  that  -  '0  yields 

rL  -  -r*  exp  [(-2iknQ  cos  &g)A]  (4-61) 

where  *  denotes  the  conjugate. 

Consider  Figure  4-4.  Energy  conservation  requires  that  the  z  com¬ 
ponent  of  the  Poynting  vector  of  the  transmitted  wave  equal  the  z 
component  of  the  Poynting  vector  of  the  incident  wave  plus  the  z  compo¬ 
nent  of  the  Poynting  vector  of  the  reflected  wave.  This  condition  is 
expressed  mathematically  by  Eq.  (4-62). 


1 


nQ  cos2  0' 


,  s  polarization 


1  -  tan2  Oq 


n 


0 


,  p  polarization 


=  1 


(4-6; 


Similarly  from  Figure  4-5,  energy  conservation  dictates 

.irj2  +  )  t  L  j  2  /  g  =  1 


t 
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where 


g  =  o(a)  cos  0g/n(O)  cos  0q,  for  s  or  p  polarization 


But 


n( 0)  -  n(A)  and  0q  -  0Q 

[see  Eqs.  (4-8a)  and  (4-42)].  Thus  Eqs.  (4-62)  and  (4-63)  reduce  to 


1‘rI  -V'  -  lrRl2 

(4-64) 

|tLl  -Vl  -  |rL|2 

(4-65) 

From  Eq.  (4-61)  one  gets 

1 rL 1  -  IV 

(4-66) 

Combining  Eqs.  (4-65)  and  (4-66) 

■  I'-rI2 

(4-67) 

It  remains  only  to  determine  the  phases  associated  with  tR  and  t^. 

We  do  this  with  the  aid  of  the  reversibility  theorem  of  stratified  media  [1] 
which  is  given  without  proof  in  Eqs.  (4-68a)  and  (4-68b). 

rR*R  *  Vl  =  0 

(4-68a) 

rL*i  +  lLrR  °  0 

(4-68b) 
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Combining  Eqs.  (4-61),  (4-68a)  and  (4-68b)  immediately  yields 

tR  -  |tR|  exp  [( -iknQ  cos  0q)a]  (4-69) 

tL  -  ) tL |  exp  [(-iknQ  cos  0q)a]  (4-70) 

From  Eqs.  (4-61),  (4-64),  (4-67),  (4-69)  and  (4-70)  it  can  easily  be 
shown  that 

tLtR  "  rRrL  r  exp  K-2ikn0  C0S  90^  (4-71) 

We  will  use  this  result  later. 

From  Eqs.  (4-53)  through  (4-56)  and  Eq.  (4-44)  we  have 


Q 


-r, 


rR  tLtR  "  rRrL 


(4-72) 


Equations  (4-61),  (4-71),  and  (4-72)  can  now  be  combined  to  yield 


where 


1  r*  exp  [(-2iknQ  cos  0Q)A] 
rR  exp  [(-2iknQ  cos  0q)A] 


rR- 


a 


dn(z) 

dz 


exp  [(-2iknQ  cos 


e0)Z] 


dz 


(4-73) 


(4-57) 
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tR  r  V'  -  lrR|2  exP  C(-i knQ  cos  e0)]  (4-74) 

Equation  (4-74)  is  an  immediate  consequence  of  Eqs.  (4-64)  and  (4-69). 

In  the  previous  section  the  two  by  two  transfer  matrix  of  a  slab 
of  thickness  A  adjacent  to  dielectrics  which  index  match  the  slab  at 
its  front  and  rear  surfaces  was  derived.  The  derivation  required  that 
the  magnitude  of  the  right  going  reflection  coefficient  r^  of  that  slab 
be  small.  Let  us  now  divide  the  grating  shown  in  Figure  4-1  into  N  slabs, 
each  of  which  has  a  small  right  going  reflection  coefficient. 

The  slabs  will  be  numbered  1  through  N,  with  the  slab  at  the  front 
surface  of  the  grating  being  the  first  and  the  slab  at  the  rear  surface 

it  i.L 

being  the  N  .  The  transfer  matrix,  ,  of  the  jL  slab  is  given  by 
Eqs.  (4-57),  (4-73)  and  (4-74)  which  are  repeated  below. 


Q,  =  — 


1  rRj  exp  [(-2iknQ  cos  60j)Aj] 

rR  exp  [(-2iknQ  cos  90j)Aj] 

J 


(4-73) 


Rj 


/  dn  .  ( z ) 

-ir- 


exp  C(-2i knQ  cos  eQ ^ ) z]  dz  (4-57) 


tR  -  J^l  -|rR  | 2  exp  C(-iknQ  cos  90j)Aj] 
j  J 


(4-74) 


The  refractive  Index  profile  of  the  grating  will  be  designated 
by  n(z)  as  before.  Thus 


,i 
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n^z)  =  n  ( z  + 


n-j  (z )  =  n(z) 


xU) 


(4-75) 


Furthermore,  from  Snell's  Law  and  Eq.  (4- 8a) 


-l  nf  sin  6n 

eoj  =  1  <  -Sr — a  5  J*'. 


,  N 


(4-76) 


Recall  that  nf  is  the  refractive  index  of  the  dielectric  in  contact 
with  the  front  surface  of  the  grating  while  n^  is  the  refractive  index 
of  the  dielectric  in  contact  with  the  back  surface  of  the  grating.  It 
should  be  noted  that  the  transfer  matrix,  ,  of  the  first  slab  is 
given  by  Eq.  (4-73)  only  when  nf  *  n(0).  This  is  because  Eq.  (4-73) 
was  derived  under  the  condition  that  the  jl  slab  is  adjacent  to  a  di¬ 
electric  which  index  matches  the  slab  at  its  front  and  rear  surfaces. 
Similarly,  the  transfer  matrix,  Q^,  of  the  Nth  slab  is  given  by  Equa¬ 
tion  (4-73)  only  when  nfa  =  n(D).  We  circumvent  this  problem  by  con¬ 
ceptually  adding  two  additional  "slabs"  to  the  grating,  one  at  its  front 
surface  and  one  at  its  back  surface.  These  slabs  are  of  thickness  6 
and  have  refractive  indices  n(0)  and  n(D),  respectively.  The  situation 
is  illustrated  in  Figure  4-7.  In  this  case  Q.j  and  QN  as  given  by  Equa¬ 
tion  (4-73)  are  valid.  As  6  shrinks  to  zero  the  new  grating  shown  in 
Figure  4-7  becomes  indistinguishable  from  the  original  one  shown  in 
Figure  4-1.  All  that  is  necessary  then  is  to  determine  the  transfer 
matrices  of  the  two  additional  layers  shown  in  Figure  4-7  in  the  limit 
as  <5  -*■  0.  This  has  already  been  done  in  Chapter  3  (see  Eq.  (3-79)). 

The  transfer  matrix  of  the  front  layer  will  be  designated  while  the 
transfer  matrix  of  the  back  layer  will  be  designated  G^. 

The  transfer  matrix,  Q,  of  the  entire  grating  will  then  be 
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(4-77) 


The  right  going  reflection  coefficient,  r^,  of  the  grating  is  given 
by  Eq.  (4-49) 


9?! 

rD  =  -=4.  ,  for  s  or  p  polarization  (4-49) 
R  qll 


and  the  right  going  transmission  coefficient  by  Eq.  (4-50). 


- —  ,  for  s  polarization 

qll 


1  cos  eo 

- —  *  ---  -a-,-  ,  for  p  polarization 
q11  cos  °0 


(4-50) 


We  have  seen  that  the  transfer  matrix  relates  the  right  and  left 
going  electric  field  immediately  outside  of  the  front  surface  of  a  grat¬ 
ing  to  the  right  and  left  going  electric  field  immediately  outside  of 
the  back  surface  of  the  same  grating. 

The  transfer  matrix  will  depend  on  the  refractive  index  of  the 
material  which  is  adjacent  to  the  front  surface  of  the  grating  and 
the  refractive  index  of  the  material  which  is  adjacent  to  the  back 
surface  of  the  grating.  It  is  more  convenient  to  redefine  the  transfer 
matrix  so  that  it  relates  the  right  and  left  going  electric  field  imme¬ 
diately  inside  of  the  front  surface  of  the  grating  to  the  right  and  left 
going  electric  field  immediately  inside  of  the  back  surface  of  the  grating. 
Thus,  the  redefined  transfer  matrix  no  longer  depends  on  the  refractive 
indices  of  mediums  which  are  adjacent  to  the  grating.  This  redefined 
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transfer  matrix  is,  of  course,  identical  with  the  original  transfer 
matrix  if  the  two  mediums  which  are  adjacent  to  the  grating  match  it 
in  refractive  index  at  the  interfaces. 

Now  suppose  that  two  gratings,  called  one  and  two,  are  placed  to¬ 
gether  with  their  surfaces  in  contact.  Gratings  one  and  two  have  trans¬ 
fer  matrices  Qj  and  respectively.  What  is  the  transfer  matrix  Q1  and  2 
of  the  combined  grating?  If  the  refractive  index  is  continuous  across  the 
grating  to  grating  interface,  then  the  right  and  left  going  fields  must 
be  continuous  across  this  interface.  Thus,  it  immediately  follows  that 

Q,  and  2  -  Q,  -  Q2  (4-78) 

If  the  refractive  index  is  not  continuous  across  the  grating  to 
grating  interface,  then  the  right  and  left  going  fields  are  not  contin¬ 
uous  across  the  interface.  Therefore,  Eq.  (4-78)  is  no  longer  valid. 

Let  |G^  denote  a  matrix  which  relates  the  right  and  left  going  fields 
immediately  inside  the  back  surface  of  optical  element  one  to  the  right 
and  left  going  fields  immediately  inside  the  front  surface  of  optical 
element  two.  Then 


and  2  ■  «1  '  1G2  '  l4-79> 

Clearly,  ^  relates  the  right  and  left  going  fields  on  one  side  of  a 
dielectric  to  dielectric  interface  to  those  on  the  other  side.  The  form 
of  ^  was  derived  in  Chapter  3  (see  Eq.  (3-79)). 

Note  that  In  Eq.  (4-77)  the  Q j ' s  can  be  multiplied  directly  to¬ 
gether  together  without  the  use  of  ^G^'s.  This  follows  from  the  fact 
that  there  is  no  discontinuity  in  the  refractive  index  at  the  inteface 
between  the  jth  and  j+lth  slabs. 
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Notes  to  Chapter  4 

[1]  Knittl,  Optics  of  Thin  Films,  John  Wiley  &  Sons,  New  York,  1976. 

[2]  E.  Wolf,  Progress  in  Optics,  Vol .  5  (Chapter  5),  North-Hoi  land 
Publishing  Co.,  Amsterdam,  1966. 

[3]  The  propagation  angle,  eo  cannot  excede  the  critical  angle  of  the 
holographic  material.  This  critical  angle,  6  ,  is  given  by 


[4]  A.  Vasicek,  Optics  of  Thin  Films,  North-Holland  Publishing  Co., 
Amsterdam,  1960. 
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QUASI -SINUSOIDAL  REFRACTIVE  INDEX  PROFILES 
AND  THE  THICK  SLAB  DECOMPOSITION 

In  Chapter  4  a  technique  for  determining  the  two  by  two  transfer 
matrix  [and  consequently  Hr(a,  6q)  and  Ht(a,  6g)]ofa  flat,  holographic 
grating  having  a  refractive  index  variation  only  in  a  direction  normal 
to  the  grating  surface  was  presented.  Here,  we  will  apply  this  method, 
known  as  Thick  Slab  Decomposition  (TSD),  to  several  holographic  gratings 
having  quasi-sinusoidal  refractive  index  profiles.  The  results  will  be 
compared  with  those  obtained  by  the  multilayer  dielectric  stack  approach. 

We  will  also  derive  simple  formulas  for  HR(\,  eQ)  and  Ht(a,  eQ)  when 
the  refractive  index  profile  of  the  grating  is  sinusoidal. 

The  following  observations  can  be  made  concerning  the  TSD 

1.  The  holographic  grating  is  divided  into  N  slabs.  Each 

of  these  slabs  must  be  sufficiently  thin  such  that  the  mod¬ 
ulus  squared  of  the  right  going  reflection  coefficient  is  much 
less  than  one.  This  condition  places  a  lower  limit  on  the  value 
of  N. 

2.  The  number  of  matrix  multiplications  required  by  the  TSD  method 
is  proportional  to  N. 

3.  The  TSD  method  offers  computational  advantages  over  the  multi¬ 
layer  dielectric  stack  approach  only  when  the  right  going  reflection 
coefficient  of  each  thick  layer  can  be  evaluated  without  sub¬ 
stantial  calculation. 

4.  The  evaluation  of  the  right  going  reflection  coefficient  of 
each  slab  generally  becomes  simpler  as  the  number  of  slabs  N, 
increases . 
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as 


We  assume  that  the  refractive  index  profile,  n(z)  can  be  written 

m 


n(z)  =  nQ  +  A(z)  sin  [2rf(z)z  +  r(z)]  (5-1) 

It  was  noted  in  Chapter  4  that  refractive  index  profiles  of  reflection 
gratings  which  can  be  holographically  fabricated  will  not  have  frequency 
components  outside  of  the  range  4  pm-1  to  6  ym"1 .  Thus,  n(z)  is  to  an 
extent  narrowband  with  f(z)  nominally  in  the  range  4  y  m”^  to  6um’'. 
Therefore,  A(z),  f(z)  and  r(z)  will  be  approximately  constant  over  an 
interval  of  length,  ,  where 

5  =  fiW  (5-2) 

If  the  grating  is  divided  into  N  slabs,  the  jth  one  being  a^  ym  thick 
(Aj  <  6)  then  the  refractive  index  profile  within  the  jth  slab  can  to  a 
reasonable  approximation  be  written  as 

Rj(z)  s  nQ  +  Aj  sin  (2nfjZ  +  r^)  (5-3) 

where  A.,  f.  and  r.  are  constant. 

J  3  J 

Below,  we  derive  the  right  going  reflection  coefficient  rR  ,  of  a 

J 

slab  A.  ym  thick  which  has  a  refractive  index  profile  given  by  Eq.  (5-3). 
J  2 

We  assume  for  the  moment  that  |rR  |  «  1.  From  Eqs.  (4-27a),  (4-27b)  and 
(4-57)  of  Chapter  4.  J 


AJdn  .(z) 
.  J 
dz 


exp  [(-2iknQ  cos  eQ)z]dz 


(5-4) 
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where 


no  cos  e0 


1  -  tan  e, 


for  s  polarization 


for  p  polarization 


(5- 5a) 


( 5- 5b) 


Substituting  Eq.  (5-3)  into  Eq.  (5-4)  yields 


Aj 

R  -  I  2irfjfl.  cos  +  rj)  exp  C  ( “2i  knQ  cos  @0)z]dz  (5-6) 


-2tt  f  .A. a  f  Aj 


iVT  or 

exP  ( i  r  j )  exp  [  i  (  2tt  f  ^  -  2knQ  cos  0Q )  z ) 

n 


+  exp  ( -ir j )  exp  [-i(2irfj  +  2knQ  cos  eQ)z] 


dz 


-2Trf.A.a  exp  [i(2irf.  -  2kn_  cos  en)A.]  -  1 

- J-J~  exD  (ir  \ _ _  J  -  0 _ 0  J 

2  P  'irj' i  ( 2tt f j  -  2knQ  cos  6g) 


2w f j^ja 

- 2 -  exp  Mrj) 


exp  [-i  ( 2tt f  j  +  2knQ 
-  i  (  2tt  f  j  +  2knQ 


cos  90)a.]  -  1 
cos  e"0) 


(Equation  continued  on  next  page) 
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=  ■1rfjAjoAj  exT  (i tr j  +  (*fj  -  knQ  cos  90)Aj])  sine  [("f.  -  knQ  cos  e0)Aj] 

_1Tf3AjaAj  exp  ("itrj  +  (irfj  +  kn0  cos  e0)Aj^)  sinc  +  kno  cos  0o)a 

The  above  equation  is  repeated  below  for  reference 

rR  =  -TrfjA^aAj  exp  (i [r ^  +  (tt f j  -  knQ  cos  eQ)Aj.])  sinc  [(irf  -  knQ  cos  9q)a.] 
J 

-TTfjAj.aA  exp  (-i[r\  +  iirfj  +  knQ  cos  eQ)A^])sinc  [ (irf ^  +  knQ  cos  9q)a j ] 

(5-7) 

From  Eq.  (5-7)  we  have 


R.1 

J 


.  2 
sinc 


[(*fj  -  knQ  cos  e0)Aj] 


+  sinc 


2 


[(^fj  +  knQ  cos  eQ)Aj] 


<  2(nfjAjaAj)2 


(5-8) 


As  in  Chapter  4,  the  index  modulation.  A.,  is  assumed  to  be  less  than 

J 

0.03.  Therefore,  if  A-  is  chosen  such  that  A.  =  1/f.  then 

J  J  J 

| rD  |2  <  2*2(f.)2(0.03)2a2(l/f,)2  =  0.018a2  (5-9) 

Kj  ~  J  J 

Since  9^,  the  propagation  angle  within  the  grating,  can  never  exceed  the 
critical  angle  (42°)  and  since  nQ  is  nominally  1.5,  Eqs.  (5-5a)  and  ( 5- 5b) 
indicate  that 
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a  <  1.46 


(5-10) 


Combining  Eqs.  (5-9)  and  (5-10)  yield 


I rR  |  <  0.027 

Hence,  our  original  assumption 


(5-11) 


I  «  1 

J 

was  valid. 

We  combine  Eq.  (5-7)  with  Eqs.  (4-73)  and  (4-74)  of  Chapter  4.  The 
results  is 


^  '  tc 


exp  C(-2iknQ  cos  eQ)A..] 


(5-12) 


where 


C  =  -irf jO&jAj  exp  (i Cr j  +  (irf j  -  knQ  cos  Sg^-O' 


sine  [(itfj  -  knQ  cos  9q)a ^3  - 


irfjAjaAj  exp  (-iCTj  +  (irf j  +  knQ  cos  0 


sine  [(irfj  +  knQ  cos  6g)&j] 


( 5- 1 3a  ) 


109 


B  =  -"fjAjaAj  exp  (-i C r j  +  ( nf j  +  knQ  cos  0Q) Aj]) * 

sine  [Ufj  -  knQ  cos  eQ)Aj]  - 

TtfjA^aAj  exp  (i  [Tj  +  (ufj  -  knQ  cos  eQ' Aj])  • 

sine  [(ufj  +  knQ  cos  eQ)Aj]  (5-1 3b) 

tR  =  V1  -  lc|2  exp  [(-iknQ  cos  0q) Aj 3  (5-13c) 

Equations  (5-12),  (5-13a),  ( 5-1 3b)  and  (5-13c)  give  the  transfer  matrix, 
Q.,  of  a  slab  A.  microns  thick  which  has  a  refractive  index  profile  n.(z) 

J  J  J 

described  by 

n(z)  =  nQ  +  A.  sin  (2irf  .z  +  r  .) 

J  J  J 

| A, |  <  0.03  and  4  ym’^  <  f.  <  6  ym  ^  (5-14) 

J  J  ” 

The  results  obtained  so  far  are  summarized  below.  We  consider  a 
flat,  holographic  grating  of  thickness  D  microns  which  has  a  refractive 
index  that  varies  only  in  a  direction  normal  to  the  grating  surface. 

It  is  assumed  that  the  refractive  index  can,  to  a  good  approximation,  be 


Jim 


RADAR  AND  OPTICS  DIVISION 


n(z)  =  n0  ♦  Aj  sin  (2.fjZ  *  Tj)  for  £  4,,  <  z  <  £  4.  (5-15) 


(noteY*  a  -  0) 

L—t  m 


where  nQ  is  nominally  1.5  m=l 

| A  - 1  <  0.03 

.  -1  ,  ,  - 1 
4  ym  <  f .  <  6  nm 

Then  Eq.  (5-12)  is,  to  a  good  approximation,  the  transfer  matrix  of  the 
section  of  the  HOE  between 


j-1  j 

T.  \  ;  2  <  I  4m 

m=  m=l 


The  two  by  two  transfer  matrix,  Q,  of  the  grating  as  well  as  rR  .  HR(x, 
eQ),  tR  and  H^.(x ,  eQ)  of  the  HOE  will  be  given  by  Eqs .  (4-77),  (4-49), 
(4-2a),  (4-50)  and  (4-2b)  of  Chapter  4.  These  equations  are  repeated 
below. 


qll  q12 


-  s  (n  ^  Gt 


(5-16) 


h21  m22 


q2i 

Hr(X,3q)  =  rR  =  - —  ,  for  s  or  p  polarization 


(5-17) 


111 


HT(x ,  0 0 ) 


1_  _  cos  eQ 
11  cos  0A 


RADAR  AND  OPTICS  DIVISION 


for  s  polarization 


— ,  for  p  polarization 


(5-18) 


It  should  be  noted  that  Eq.  (5-12)  is  valid  regardless  of  how  large  Ai 

2 

is  provided  that  |r0  I  <<  1 .  Recall  that  Eq.  (5-8)  gives  an  upper  bound 
2  j 

for  | rD  ]  which  can  be  used  to  determine  the  maximum  allowable  a.. 

Rj  3 

It  can  be  shown  that  calculation  of  Hr(a,  6g)  and  Ht(a,  eQ)  for  ten  or, 

more  different  values  of  6g  or  \  via  Eqs.  (5-12),  (5-16),  (5-17)  and 

(5-18)  will  be  at  least  an  order  of  magnitude  faster  than  performing  the 

same  caluclations  using  the  multilayer  dielectric  stack  approach.  Equation 

(5-12)  can  be  further  simplified  under  tl?e  following  conditions: 

1.  a.  is  chosen  such  that  f-A.  is  an  integer 
J  J  J 

Note  that  for  f,A.  an  integer,  ri  must  be  0  or  n  since  n(z)  is  a 

J  J  J 

continuous  function  of  z.  When  the  above  condition  is  satisifed  Eq.  (5-12) 
becomes 


i  -3j 


exp  [( i knQ  cos  eQ)A j ] 


exp  [(-ikn0  cos  Wq)a^] 


a j  exp  [{-iknQ  cos  eQ)Aj]  exp  [(-2iknQ  cos  eQ)L. 


(5-19) 


where 
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IT  f  j  A  *CtA  • 
J  J  J 


(*' 


nc 


knA  cos 


+  sine  [ ( tt f j  +  knQ  cos 


(5-20) 


The  plus  sign  corresponds  to  fjAj  an  odd  integer  and  =  0,  or  fjA^ 
an  even  integer  and  r.  =  it.  The  minus  sign  corresponds  to  f.A.  an  odd 

J  J  J 

integer  and  r.  =  u,  or  f.A.  and  even  integer  and  r.  =  0.  In  the  remain- 
J  J  J  J 

der  of  this  report  we'll  assume  that 

A.  -  1/f .  (5-21) 

J  J 

Then  Eqs.  (5-19)  and  (5-20)  can  be  rewritten  as 


Q 


0 


exp  (i<i>j ) 

Trr^i 


exp  (-i$j) 


a.  exp  ( - i 4> j )” 
exp  ( -2i<j> j )_ 


(5-22) 


where 

•(•j  =  (knQ  cos  eQ)/fj  (5-23) 

aj  =  7rAj01  (sin  "  kno  cos  e0)/fj]  + 

sine  [ ( tt f j  +  knQ  cos  eQ)/fj])  (5-24) 

The  above  three  equations  are  the  principal  result  of  this  chapter. 
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On  the  following  pages  plots  are  given  of  the  reflectance,  ( |HR( 


»Q)  I  ,  phase  change  upon  reflection 


y>-  v 


2 

transmittance  (|H^.(x,  e q )  |  ),  and  phase  change  upon  transmission 

ht(a,  e0) 

|ht(x,  e0)| 

vs.  angle  of  incidence  for  three  flat  holographic  optical  elements 
under  s  polarized  plane  wave  illumination.  All  of  the  gratings  con¬ 
sidered  have  a  refractive  index  which  varies  only  in  a  direction  normal 
to  the  surface  of  the  grating,  and  for  simplicity,  they  are  assumed  to 
be  surrounded  by  a  homogeneous  dielectric  of  refractive  index  1.5.  The 
plots  were  generated  using  both  the  multilayer  dielectric  stack  and  thick 
layer  decomposition  techniques. 

Each  of  the  elements  has  a  refractive  index  profile  which  to  a 
good  approximation  can  be  given  by  Fq.  (5-15)  with  the  r/s  =  0  and  nQ  =  1.5. 
Figures  5-1,  5-6  and  5-11  are  plots  of  the  refractive  index  modulation 
A,  (solid  line),  and  the  frequency  f.  (dashed  line)  as  a  function  of  the 

w  J 

cycle  number  j.  Figures  5-2  through  5-5  are  plots  of  the  reflectance, 
phase  change  upon  reflection,  transmission  and  phase  change  upon  trans¬ 
mission  vs.  angle  of  incidence  for  the  element  whose  refractive  index 
profile  is  given  in  Figure  5-1.  The  refractive  index  profile  of  this 
element  can  be  analytically  described  by  Eq.  (5-25)  below 

n(z)  =  1.5  +  0.02  sin  ( 2 tt - 5 . 1 4237 - z )  (5-25) 
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Figure  5-1.  Refractive  Index  Modulation  and  Frequency  vs. 
Cycle  Number  (Example  1) 
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ANGLE  OF  INCIDENCE  IN  OEGREES 

Figure  5-2.  Comparison  of  Reflectance  vs.  Angle  of  Incidence 
Characteristics  Computed  by  the  TSD  Method  and  the 
Multilayer  Dielectric  Stack  Method  (Example  1) 
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Figure  5-4.  Comparison  of  Transmittance  vs.  Angle  of  Incidence 
Characteristics  Computed  by  the  TSD  Method  and  the 
Multilayer  Dielectric  Stack  Method  (Example  1) 
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Figure  5-6.  Refractive  Index  Modulation  and  Frequency 
Cycle  Number  (Example  2) 
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Figure  5-7. 


Comparison  of  Reflectance  vs.  Angle  by  Incidence 
Characteristics  Computed  by  the  T5D  Method  and 
the  Multilayer  Dielectric  Stack  Method  (Example  2) 
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where  z  is  measured  in  microns  and  the  element  thickness  is  14.58  pm. 

The  solid  lines  in  Figures  5-2  through  5-5  correspond  to  results  obtained 
by  the  thick  layer  decomposition  technique,  while  the  dashed  lines  cor¬ 
respond  to  results  obtained  by  the  multilayer  dielectric  stack  technique. 
The  thickness  of  each  layer  in  the  multilayer  dielectric  stack  technique, 
was  taken  to  be  0.0039  pm.  It  is  noted  that  the  two  techniques  yield 
results  which  agree  quite  closely.  The  set  of  very  short  dashed  lines, 
which  appear  as  dots  in  Figure  5-3,  indicate  a  phase  change  upon  re¬ 
flection  of  approximately  356°  (or  equivalently  -4°).  The  discontinuities, 
in  the  phase  change  upon  reflection  shown  in  Figure  5-3,  are  physically 
acceptable  since  they  occur  only  at  angles  where  the  reflectance  is  zero. 

By  examining  the  numerical  values  of  the  points  plotted  in  Figure  5-3, 
one  would  find  that  the  multilayer  dielectric  stack  theory  predicts  that 
the  discontinuities  occur  at  phase  changes  of  approximately  -4°  and  176°. 
The  thick  layer  decomposition  theory  predicts  that  the  discontinuities 
occur  at  phase  changes  of  approximately  0°  and  180°. 

Figures  5-7  through  5-10  are  plots  of  the  reflectance,  phase  change 
upon  reflection,  transmittance,  and  phase  change  upon  transmission  vs. 
angle  of  incidence  for  the  element  whose  refractive  index  profile  is 
given  in  Figure  5-6.  The  refractive  index  profile  of  this  element  can 
be  analytically  described  by  Eq.  (5-26)  below 


n(z)  =  1.5  +  0.01  ^1  +  sin  ( 2tt . f .  1 42 37 - z ) 


(5-26) 


where  z  is  in  microns  and  the  element  thickness  is  6.99  ym.  The  solid  lines 
in  Figures  5-7  through  5-10  correspond  to  results  obtained  by  the  thick 
layer  decompositon  technique,  while  the  dashed  lines  correspond  to  re¬ 
sults  obtained  by  the  multilayer  dielectric  stack  technique.  The  thickness 
of  each  layer  in  the  multilayer  dielectric  stack  technique  was  taken  to  be 
0.005  ym.  Once  again,  the  two  techniques  yield  results  which  agree  quite 
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Finally,  Figures  5-12  through  5-15  are  plots  of  the  reflectance, 
phase  change  upon  reflection,  transmittance,  and  phase  change  upon  trans¬ 
mission  vs.  angle  of  incidence  for  the  element  whose  refractive  index 
profile  is  given  in  Figure  5-11.  The  refractive  index  profile  of  this 
element  does  not  have  a  simple  analytic  description.  The  solid  lines 
in  Figures  5-12  through  5-15  correspond  to  results  obtained  by  the  thick 
layer  decomposition,  while  the  dashed  lines  correspond  to  results  obtained 
by  the  multilayer  dielectric  stack  method.  The  thickness  of  each  layer  in  the 
multilayer  dielectric  stack  technique  was  taken  to  be  0.005  ym.  Both  tech¬ 
niques  yield  results  which  are  in  good  agreement.  Recall  that  a  number 
of  restrictions  were  placed  on  the  refractive  index  profile  in  order  to 
justify  the  approximations  used  in  the  development  of  the  thick  layer 
decomposition  theory.  One  of  these  assumptions  placed  a  limit  of  0.03  on 
the  maximum  refractive  index  modulation.  This  limit  is  exceeded  for  the 
grating  having  the  refractive  index  profile  shown  in  Figure  5-11.  De¬ 
spite  this  fact,  plots  in  Figures  5-12  through  5-15  indicate  that  the 
theory  still  yields  fairly  accurate  results  in  this  particular  case. 

In  this  section,  we  will  derive  analytic  expressions  for  HR(>,  eQ) 
and  Ht(x,  0q)  of  a  flat,  holographic,  reflection  grating  having  a  re¬ 
fractive  index  profile,  n(z),  given  by 

n(z)  =  nQ  +  A  sin  ( 2tt f z )  (5-27) 

For  convenience,  the  grating  is  assumed  to  be  surrounded  by  a  dielectric 
of  refractive  index  and  the  thickness,  D,  of  the  grating  is  taken  to 
be  an  integer  multiple  of  1/f  (i.e.,  D  =  N/f  where  N  is  some  integer). 

We  start  by  stating  the  following  theorem  from  linear  algebra  [2] 
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Theorem 

Let  B 


B  = 


J11 


J21 


1 2 


22 


(5-28) 


be  a  two  by  two  unimodular  matrix  (i.e.,  ^1^22  -  2 ^21  =  !)•  Then  the 


,th 


N  power  of  B  for  an  arbitrary  integer  N  is  given  by 


/bl  1  Fn-1  ‘  Fn-2  b12FN-1 


,b21  FN-1 


where 


b22FN-l  "  FN-2 J 

cos'1  (  2-(bn  +  b22)),  ’  (bn  +  b22) 
cosh  (|(b11  +  b22))  p  (bi ;  +  b?2 ) 


1 


(5-29) 


(5-30) 


FN  ‘ 


(sin  (N  +  1  )  >  )/sin  , ,  ^bll  +  b22)  '  -  1 


(sinh  (N  +  1 )) ,/sinh  t,  '^(b^  +  b??)  1 


tquation  (5-27)  indicates  that  the  f^’s  and  A^’s  of  Eqs .  (5-22),  (5-23), 
and  (5-24)  are  independent  of  j.  Therefore,  the  Q.'s  of  Eq.  (5-22)  are 
also  independent  of  j.  Consequently,  Eqs.  (5-22),  (5-23),  and  (5-24)  can 
be  rewritten  as 
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where 

=  (knQ  cos  0Q)/f  (5-33) 

dj  =  a.|  =  Aa^sinc  [ ( ..- f  -  knQ  cos  *)Q)/f]  + 

sine  [(irf  +  kn  cos  OgJ/f])  (5-34) 


We  note  that  Qj  is  unimodular.  It  follows  from  the  theorem  above  that 


1  33 


Now  by  Eqs .  (5-16),  (5-17)  and  (5-18) 


Q  = 


hr(x,  e0) 


(5-36) 


(5-37) 


HT{A’  V  =  exp  (i^ 


(5-38) 


VT-  a, 


_L_  p  -  F 
—  rN-l  hN-2 


Equations  (5-37)  and  (5-38)  above  are  completely  general  expressions  for 
Hr(-,  <'q)  and  Hy( x ,  Og).  On  the  other  hand,  analytic  expression  for 
HR(  ,  ,  *)g)  and  Ht(a,  p*q)  derived  from  a  coupled  wave  theory  approach  re¬ 
quire  "g  to  either  approximately  or  exactly  satisfy  the  Bragg  condition 
and  also  require  that  the  transmitted  wave  be  strongly  depleted. 
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Notes  to  Chapter  5 

[1]  See  Chapter  8  for  a  discussion  of  this  assumption. 

[2]  P.  Yeh,  A.  Yariv,  C.  Hong,  "Electromagnetic  Propagation  in 
Periodic  Stratified  Media.  I.  General  Theory,  Journal  of 
the  Optical  Society  of  America,  Vol .  67,  No.  4,  April  1977 


135 


RADAR  AND  OPTICS  DIVISION 


6 

GENERAL  PROPERTIES  OF  HOLOGRAPHIC 
REFLECTION  GRATINGS 

In  the  previous  chapters  we  developed  techniques  for  determining 
the  diffracted  field  from  a  flat,  holographic,  phase,  reflection  grat¬ 
ing  which  has  a  retractive  index  variation  only  in  the  z  direction. 

The  refractive  index  was  assumed  to  be  a  continuous  function  of  z.  In 
this  chapter  we  will  simplify  the  form  of  the  transfer  matrix,,  prove 
a  theorem  concerning  the  mutual  interaction  between  two  gratings  when 
placed  in  contact,  discuss  the  effects  of  polarization,  indicate  the 
interrelationship  between  wavelength  and  incident  angle,  and  determine 
a  formula  for  the  transfer  matrix  which  does  not  involve  matrix 
multipl i cations . 

6.1  SIMPLIFICATION  OF  THE  TRANSFER  MATRIX 

We  start  by  considering  a  flat,  holographic,  phase  reflection 
grating  which  has  a  refractive  index  variation  only  in  the  z  direction. 
This  refractive  index  variation  is  arbitrary  and  not  necessarily  contin¬ 
uous.  The  grating  is  assumed  to  be  of  infinite  extent  in  its  x  and  y 
dimensions.  The  front  and  rear  surfaces  are  adjacent  to  dielectrics 
having  refractive  indices  n^  and  respectively  (see  Figure  6-1).  In 
all  future  dicussions,  the  word  grating  will  refer  to  a  structure  having 
the  above  characteristics. 

Many  of  the  results  derived  in  Chapter  4  do  not  explicitly  depend 
upon  the  form  of  the  refractive  index  variation.  In  particular,  the 
grating  or  any  slab  thereof  can  be  characteri zed  by  a  two  by  two  trans¬ 
fer  matri x  Q. 
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Refractive  Index  =  n(z) 


Figure  6-1.  Holographic  Reflection  Grating 
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Hrim 


Q  = 


-r. 


rR  tLtR  ‘  rRrL 


(6-1) 


where  r  ,  rL,  tp,  and  tL  are  the  right  going  reflection  coefficient, 
left  going  reflection  coefficient,  right  going  transmission  coefficient, 
and  the  left  going  transmission  coefficient,  respecti vely.  Recall  that 
r^  and  tp  are  functions  of  the  free  space  wavelength,  a,  and  the  inci¬ 
dent  angle,  0q,  at  the  front  surface  of  the  grating,  r^  and  t^  are  func 
tions  of  the  free  space  wavelength  and  the  incident  angle  at  the  back 
surface  of  the  grating.  We  wish  to  simplify  the  form  of  the  transfer 
matrix  Q.  The  incident  angles  at  the  front  and  back  surfaces  are  as¬ 
sumed  to  be  related  by  Snell's  Law. 


The  reversibility  theorem  of  stratified  media  is  given,  without 
proof,  by  Eqs.  (6-2)  and  (6-3)  below.  [1] 


rRtR  +  V*  =  0 

(6-2) 

rLt*  +  tLr*  =  0 

(6-3) 

rR’  rL * 
through 

tD  and  t.  can  be  written  in  polar  form. 

K  L 

(6-7)  below. 

We  do  so  in  Eqs.  (6- 

rR  “  lrRl  e>1>  "V 

(6-4) 

rL  =  |rL|  exp  (i,L) 

(6-5) 

tR  -  |tR|  exp  (1rR) 

(6-6) 
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tL  =  |tL|  exp  (irL) 


(6-7) 


It  is  worthwhile  at  this  point  to  make  explicit  and  to  justify 
the  underlying  assumption  that  has  been  made  throughout  this  report. 
Suppose  that  the  front  surface  of  the  optical  element  in  Figure  6-1 
is  illuminated  by  a  plane  wave  at  an  incident  angle  6q.  Then  the  as¬ 
sumption  is  that  the  field  at  any  plane  z  =  Zg  within  the  element  can 
be  written  as  the  field  due  to  the  coherent  superposition  of  two  plane 
waves.  Both  plane  waves  lie  in  the  plane  of  incidence  of  the  grating, 
but  one  travels  in  the  +z  direction  while  the  other  in  the  -z  direction. 
The  plane  waves  make  angles  of  6 ( Zq )  and  180  -  0(Zg)  degrees  with  the 
z-axis,  where  eU0)  is  given  by  Snell's  Law 


,  nf  sin  en 


(6-8) 


This  assumption  can  be  rigorously  justified  by  viewing  the  grating  as 
a  dielectric  multilayer  stack  composed  of  an  infinity  of  infinitesimally 
thin  layers.  The  theory  developed  in  Chapter  3  can  then  be  applied  to 
the  grating.  The  validity  of  the  above  assumption  follows  from  that 
theory.  Equations  (6-2)  and  (6-3)  become  with  the  aid  of  Eqs.  (6-4) 
through  (6-7). 

rL  =  |rj  exp  (i yl)  =  -|rR|  exp  [i(2TR  -  Yr)J  (6-9) 

rL  =  j rL |  exp  (iyL)  =  -|rR|  e*P  [i(2rL  -  yR)]  (6-10) 


It  immediately  follows  that 


irLl  ■  i<-ri 
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Furthermore,  Eqs.  (6-4)  and  (6-9)  yield 

rL  =  ~rR  exP  (i2IV 
Energy  conservation  dictates  that 

!rR I2  *  I‘rI29  =  ' 

|rL|2  +  |tL|2/9  -  t 

where 


n,  cos  Oq 

-5- - .  for  s  or  p  polarization 

nf  cos  0O’  r  K 


Combi ni ng  Eqs . 


Combining  Eqs. 


(6-6)  and  (6-15)  yields 

(1  -  |rR|2)’/Z 
tR  -  t/2  6XP  ^ 1 

(6-7),  (6-11)  and  (6-16  yields 

tL  =  (1  -  |rRj2)1/2  g1/2  exp  (irL) 


DIVISION 

(6-12) 

(6-13) 

(6-14) 

(6-15) 

(6-16) 

(6-17) 

(6-18) 

(6-19) 
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Therefore , 


Vl  =  (1  ‘  I  rRl2>  exP  Ci(rL  +  VI  (6-20) 

which  with  the  aid  of  Eq.  (6-13)  becomes 

tRtL  =  ^  "  lrR^  exp  ^i2rR^  (6-21) 

Now  by  Eq.  (6-14) 

rLrR  =  “ 1 rRi 2  exp  (i2lV  (6-22) 

Combining  Eqs.  (6-21)  and  (6-22)  one  has 

tLtR  ‘  rLrR  =  exp  (i2IV  (6-23) 

With  Eqs.  (6-14)  and  (6-23),  the  transfer  matrix,  Q,  given  by  Eq.  (6-1) 
finally  becomes 


1 


r*  exp  ( i 2 T R) 
exp  (12Tr) 


(6-24) 


6.2  THE  ELEMENT  ADDITION  THEOREM 

We  are  now  in  a  position  to  prove  a  useful  theorem. 

Theorem  6-1 : 

Let  two  gratings,  called  one  and  two,  be  placed  in  contact 
surface  to  surface.  Grating  one  has  transfer  matrix  Q1  and  corre¬ 
sponding  rD  and  tD  .  Grating  two  has  transfer  matrix  Q„ 

R1  R1  2 


f 
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Proof: 

Equation  (6-25)  was  developed  earlier  as  Eq.  (4-78)  of  Chapter  4. 
Equations  (6-26)  and  (6-27)  follow  immediately  from  Eqs.  (6-24)  and 
(6-25).  In  the  remainder  of  this  report  we  will  assume,  for  simplicity, 
that  the  refractive  indices  of  the  holographic  element  at  its  front  and 
rear  surfaces  match  those  of  the  surrounding  medium. 


6.3  POLARIZATION  EFFECTS 

The  transfer  matrix,  Q,  and  consequently  the  right  going  reflection 
and  transmission  coefficients  of  a  grating,  are  functions  of  three 
parameters  --  polarization,  wavelength,  and  angle.  It  was  shown  in 
Chapter  4  that  the  transfer  matrix,  Q,  can  be  decomposed  into  a  product 
of  transfer  matrices,  Qj,  as  given  below. 

N 

q  =  n  Qj  <6-28) 
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The  Qj's  are  given  by  Eqs.  (5-5a),  (5-5b),  (5-19),  and  (5-20)  of  Chap¬ 
ter  5.  These  equations  are  repeated  below  for  convenience. 


exp  C ( i knQ  cos  9Q )Aj) 


1  aj  exp  [(-ikn0  cos  eQ)A .] 

aj  exp  [ (-1  knQ  cos  0q)Aj]  exp  [( -2i knQ  coseQ)A^ 

(6-29) 

aj  =  +  TTfjAjOiAj^sinc  [(irfj  -  knQ  cos  0Q)Aj]  +  sinc  [(^fj  +  knQ  cos  0o^Aj-l) 

(6-30) 


nQ  cos  9q 


1  -  tan^  Qr 


•5— —  ,  for  s  polarization 

C  A  4*. 


(6-31) 


0 


for  p  polarization 


A  number  of  conclusions  can  be  drawn  from  the  equations  above.  The 
first  of  these  is  stated  as  a  theorem. 


Theorem  2 : 

Consider  two  gratings,  called  one  and  two.  Grating  one  is  illumi¬ 
nated  at  an  angle  of  incidence  9g  with  s  polarized  light  and  grating 
two  is  illuminated  at  the  same  angle  of  incidence  with  p  polarized 
light.  Grating  one  has  a  refractive  index  profile,  n(z),  given  by 

n(z)  =  nQ  +  An(z) 

j 

and  grating  two  has  a  refractive  index  profile,  n(z),  given  by 

n(z)  =  nQ  +  e  •  An(z) 
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The  constant  e  is  defined  below 


e 


tan2  eQ) 


cos20q(1 


Then  the  transfer  matrix,  Q,  of  grating  one  as  a  function  of  Ag  is  iden 
tically  equal  to  the  transfer  matrix,  Q,  of  grating  two  as  a  function 

of  Ag. 


Proof: 

Equations  (6-28),  (6-29),  and  (6-30)  indicate  that  the  transfer 
matrix  Q  is  independent  of  polarization  for  a  fixed  0g  provided 

A- a  for  s  polarization  =  A, a  for  p  polarization  (6-32) 

J  J 

Equation  (6-32)  can  be  combined  with  Eq.  (6-31)  to  yield 


A  - -  for  s  polarization  =  A. 

n0  cos  0O 

It  then  follows  immediately  that  Eq.  (6 
A.  for  p  polarization  =  A.  - - 

J  J  COS^  6g(l 


1  -  tan2  6g 

- - -  for  p  polarization 

0  (6-33) 

-32)  is  equivalent  to 

- - p -  for  s  polarization 

-  tan*  0O) 


or 


An(z)  for  p  polarization  =  e  •  An(z)  for  p  polarization 

Theorem  6-2  states  that  for  a  fixed  angle  of  incidence  changing  the 
polarization  is  completely  equivalent  to  multiplying  the  refractive 
index  modulation  by  a  constant.  Although  the  theorem  is  strictly 
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2Frim 


valid  only  for  dg  fixed,  it  gives  reasonably  accurate  results  over  any 
range  of  incident  angles  for  which  e  does  not  change  appreciably. 

6.4  THE  INTERRELATIONSHIP  BETWEEN  WAVELENGTH  AND  INCIDENT  ANGLE 

To  emphasize  that  Q.,  and  consequently  Q,  is  a  function  of  three 

J 

parameters  --  polarization,  wavelength,  and  angle  of  incidence  --  let 
us  wri te 

Q.  =  Q.  (polarization,  A  ,  e  )  (6-34) 

J  J  u 


and 

Q  =  Q  (polarization,  X  ,  0Q)  (6-35) 

We  now  make  the  following  assertion. 

Assertion  1 : 

Suppose 

cos  0n  cos  0n 

- 9-  = - 2.  (6-36) 

A  A 

and 

a(at  6q)  -  u(at  0Q)  (6-37) 


Then 


qj 


(polarization,  A 


9o> 


qj 


(polarization,  A 


e0)  (6-38) 


and 


Q  (polarization,  a  ,  0Q)  -  Q  (polarization,  '  ,  0Q)  (6-39) 
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Plausibility  Argument 

It  can  be  seen  from  Eqs.  (6-29),  (6-30),  and  (6-31)  that  the  excep¬ 
tion  of  the  single  factor,  a,  0g  and  A  appear  in  the  matrix  (and 
consequently  Q)  only  in  the  combination  cos  6g/A.  On  the  other  hand,  we 
see  by  Eq.  (6-31)  that  a  is  a  function  of  0g  but  not  A ,  Experience  indi¬ 
cates  that  small  changes  in  a  result  in  small  changes  in  Q.  (I  am  unable, 
however,  to  prove  this  analytically.)  Thus,  the  above  assertion  follows 
immedi ately . 

Now  let 


A  C0S  60 


(6-40) 


Then 


-A(sin  6n)  d6n  +  (cos  6n)  dA 
dX  =  - - - % - - 2 -  (6-41) 

h 


and  so 


dA  =  A  tan  0g  dfig 


(6-42) 


Therefore,  by  Assertion  1 

Q  (polarization,  A  +  (A  tan  eQ)  deQ,  eQ)  -  Q  (polarization.  A,  f>0  +  d-Q) 

(6-43) 

Equation  (6-43)  shows  the  connection  between  angular  selectivity  and 
wavelength  selectivity  in  a  holographic  grating.  In  particular,  it 
indicates  that  changing  the  wavelength  byA(tan  0Q)  d* jq  is  roughly 
equivalent  to  changing  the  angle  by  dOg. 
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In  the  remainder  of  this  chapter  we  prove  four  theorems  which  will 
be  used  in  the  following  chapter  to  obtain  an  approximate  solution  to 
the  inverse  scattering  problem.  The  first  two  of  these  theorems  have 
strikingly  similar  counterparts  in  the  theory  of  multilayer  dielectric 
stacks. 


6.6  DESCRIPTION  OF  THE  TRANSFER  MATRIX,  Q,  VIA  THE  PARAMETERS 

f,  AND  a . 

J  J 

Equations  (6-28)  through  (6-31)  express  tht  transfer  matrix,  Q, 


as  a  product  of  transfer  matrices  Q 


Each  matrix  0.  in  turn  is  deter- 
J 


mined  by  the  refractive  index  profile.  Without  loss  of  generality,  one 

may  assume  that  A  =  1/f..  We  wish  to  derive  an  expression  for  g  which 
J  3 

involves  only  a.'s  and  f.'s  and  no  matrix  multiplications.  First  define 
J  J 

♦j  by 


knQ  C0S  "0 


(6-44) 


Then  Eq.  (6-29)  becomes 


(.'o  = 


Qo ' 


exp  (i$j) 


a  . 
J 


i  a  exp  -  i 

J  J 

a  .  exp  (-iv  .)  exp  (-i2;  •  ) 

J  J  J 


( 6-4S 


We  make  the  definitions  given  below 
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M 

J 


1  a  exp  (-i '  )  j 

J  J  j 

dj  exp  ( - i  ;  j )  exp  ( - i )  | 

-j 


(6-46) 


The  class  of  even  products  of  order  m,  designated  Ep(m)  ,  are  ex- 
iu-pss i  ons  of  the  form 


a 


,-1 


a  exp 
in 


m- 1  j  ♦  1 

HO  •••  E 


j=0  n-Lj+l 


'0  '  0 


\i*r' 


E  «<- 


!)J  •  ')c„  ■  0  for  (j  .  1  •  >t  -1 


n  .  +  1 
J 


(6-47) 


whore  in  is  a  positive  even  integer  and  .  is  a  mono  ton i- 

I  c  m 

cally  increasing  sequence  of  positive  integers.  An  element  of  this 
class  will  be  denoted  by  ep(m).  The  class  of  odd  products  of  order  in, 
designated  op(m),  are  expressions  of  the  form 


m-1  >‘j  +  l 


-1 


a  a  ...  a 
'  1  'l  "m 


exp  +  ij  hEE 


(M)J  +  1  -l). 


j=0  n=t 
J  J 


■o  =  ° 


-1 


((-1) 


J*1 


I  )  •  '-  0  for  >  + 

"  J 


-  1 


n  = ,  ♦  1 

J 
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where  m  is  a  positive  odd  integer  and  >  •••  »  is  a  monotoni- 

cally  increasing  sequence  of  positive  integers.  An  element  of  this 
class  will  be  denoted  by  op(m). 

Using  the  definitions  above,  we  state  the  following  theorem. 


Theorem  6-3: 


where 


nv 


;=i 


'n 


21 


-(N) 

X(N) 

'll 

i  2 

,(N) 

y(  N  ) 

'21 

*2  2 

ep(m) 

’  "m  " 

1  ^  ’  ’  m 

-  N) 

,(N) 

1 2 


X<JJ*  exp 


[_i2(Z)^)] 

j=l 


X^  -  X^*  exp  [-i2(£  Oj)3 


j  =  l 


(6-50) 


(6-51) 

(6-52) 


(6-53) 


(6-54) 
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Proof: 

Theorem  6-3  will  be  proved  by  mathematical  induction.  It  is  trivial 
to  show  that  Theorem  6-3  is  true  for  N  =  1.  Now  suppose  the  theorem  is 
true  for  N  =  N  ,  then  it  is  sufficient  to  show  that  it  is  also  true  for 
N  =  N0  +  1  .  We  wri  te 


e  b 

c  d 


N  +1 
o 


n  «,  ■  <nv 

j=l  J=1 


<N0> 

1 

o 

xn 

x12 

(N1 

<No> 

CM 

X 

x22 

L, 

aN  +1  exp  <-i:N  +  1} 

O  0 


aN  +1exp('ipN  +1}  exp(-i2;Vl} 

O  0  u 


Therefore , 


( 6-56 ] 


(N  )  (N  ) 

X11  +  x12  aN  +1  6XP  +1^ 

o  o 


(6-56) 


(NJ  (N  ) 

b  =  X11  aN  +1  6XP  (_i^  +1}  +  x12  6XP  ("i2;N  +1}  (6’57) 

0  0  o 


(»„)  !"o> 

c  x21  +  x22  a^  +1  exp  ('"» 4>N  +-\ 

oo 


(6-58) 


(N  )  (N  ) 

d  =  x 2 1  °  aN  +1  exp  +1}  +  x22  exp  (_2:N  +1]  (6_59) 

0  0  0 
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It  remains  to  show  that 


e  =  x 


(Vl) 

11 


(6-60) 


c  =  x 


<v1} 

21 


(6-61) 


N  +1 
o 


b  =  c*  exp  <Pj)3 


j  =  l 


(6-62) 


N  +1 
0 

d  -  e  exp  [-12(7: 


(6-63) 


We  prove  Eq.  (6-60)  first.  Equation  (6-56)  can  be  rewritten,  by  using 
Eqs .  (6-51)  and  (6-53) 


(Nq)  (N0)' 

X11  X21 


N 


exp  [-12(7:  ^j)aN  +1  exP  +1* 


i  =  l 


=  1  +  (sum  of  all  ep(m),  l  -  NQ)  +  (sum  of  all  op(m),  £n  -  Nq)* 


exp  [-  i  2  (7. 4^)]  aM  +1  exp  (-i<*>N  +1 ) 
j-1 


(6-64) 


Consider  now  the  arbitrary  odd  product  of  order  m  given  below  (m  -  N 


0; 
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op(m)  =  a4  a£  . ..  exp 
12  m 


iU„  +  +  ...  +  4>0  )  +  i 


m-1  ‘‘j+l’1 

EE, 

J  =0  n=£  j+1 


((-Dj+1-D^ 


(6-65) 


Then 


op(tn)*  exp  [-i2(V'  40]  a  .  exp  (-1$N  ,)  = 

J  o  o 


a  a  ...  a  a  exp 
1  2  m  m+1 


-i  (4>j  +  4>5  +  . . .  +  4>p  +  4>f,  ) 

1  k2  m  m+1 


m  Vr1 


+  i 


E  E 


j=0  n=£.+l 

J 


where  l  =  N  +1 
m+ 1  o 


(6-66) 


Note  that  the  right  hand  side  of  Eq.  (6-66)  is  just  an  element  of  EP(m), 

-  N  +  1 .  Furthermore,  it  is  easy  to  see  that  all  elements  of  EP(tn) , 

-  N  +  1  and  only  those  elements  can  be  generated  by  multiplying 

some  op(m)*  £  -  N  by 

exp 


aN  +1  exP 
o 


1-1 +1^ 
o 


Therefore 


e  =  1  +  sum  of  al 1  ep(m) ,c 


+  1 


=  X 


<VD 

n 
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Next  we  will  prove  Eq.  (6-61).  Equation  (6-58)  can  be  rewritten  using 
Eqs .  (6-51)  ,  (6-52)  and  (6-54)  as 


c  =  x 


(No) 

21 


+  x 


(V* 

11 


exp  [-12  (y;  4).)]  aN  +]  exp  (-i<j>N  +] ) 

J  ^ 


(s 


=  (sum  of  al  1  op(m) ,  £  - 


N 

o 

+  (1  +  the  sum  of  all  ep(m),  £^  -  N  )*  •  exp  [-i2(^  ^  4 ■ . )] 

m  0  J  + 

J=1  < 


*  aN  +1  exp  (“1<J)N  +1}  (6"67) 

O  0 

Consider  now  the  arbitrary  even  product  of  order  m  given  below  (^  -  Nq) 

ep(m)  =  a£  a^  . . .  a£  exp  [-i(<Pt  +  +  . . .  +  <t>£  ) 

12m  12  m 

m-1  Vl'1 

+  "1)iSn]  (6'68) 

j=0  n=£.+l 
J 


Then 


ep(m)*  exp 


exp  ^'i(J,N  +1^  = 
0 


(Equation  continued  on  next  page.) 
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a 


a  a  exp  [-i($  +  <*>-  + 

m  xm+l  *1  £ 


+ 


) 


m  £j+l_1 

+ 1  2  {(_1)J+1  -  ■'K1  where  Vi  =  No + 1 

j=0  n=£ j  +1  (6-69) 


Note  that  the  right  hand  side  of  Eq.  (6-69)  is  just  an  element  of  0P(m) , 

2.  -  N  +1.  Furthermore,  it  is  easy  to  see  all  elements  of  0P(m) , 
mo 

£  -  Nq  +  1  ,  and  only  those  elements  can  be  generated  by  multiplying 

1  or  some  op(m) ,  JL  -  N  by 
'  m  o 


exp 


aN  +1  exp  (-S  +1} 
o  o 


Therefore 


c  =  sum  of  al  1 


op(ro) >  £m  -  N0  +  1 


=  x 


(No+l) 

21 


Next  we  will  prove  Eq.  (6-62).  Equation  (6-57)  can  be  rewritten  using 
Eqs.  (6-53)  and  (6-54)  as 

N„ 

(N  )* 

b  =  x22°  exp  [-i2(>  ^j)]  aN  +1  exp  (-i<J>N  +1 ) 

j  =  l  0  0 


(Equation  continued  on  next  page.) 
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{N  )’ 


+  x21°  exp  [-12(^  *.)]  exp  {-12^  +]) 

j=l  0 


(NJ1 


N  +1 
0 


v  o' 

aN  +1  exP  +1  ^  x22  exP  E~i2(  /  *  $>•;)] 

0  0  T“T  J 

J=1 


N  +1 
o 


(N0) *  ^ 

+  x21  exp  [-i2(  }  j  4>j)] 

j=l 


(6-70) 


Now  from  Eqs.  (6-58)  and  (6-70), 


N  +1 
0 

b  =  c*  exp  [-12(7;  <j»j)] 

j=l 

Finally,  Eq.  (6-6.1)  will  be  proven.  Equation  (6-59)  can  be  rewritten, 
using  Eqs.  (6-53)  and  (6-54)  as 

N 

(V*  . 

d  =  x]2  exp  [-*2(>  <f>J]  aN  +]  exp  (-i$N  +,  ) 

N 

0U*  v™' 

+  xn  exp  [-i2(>  #.)]  exp  (-i2<j>N  +]) 

J-l 


(Equation  continued  on  next  page.) 
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here  N  is  the  number  of  slabs  into  which  the  grating  has  been  divided. 
Then  by  Eqs.  (6-1)  and  (6-72) 
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TI^-‘  j 

tR  ”  1  +  sum  of  all  ep(m},  H  -  R  ex*3  ^  ^ 

m  3  =  1 

sum  of  al  1  op(m) ,  £  -  N 

rR  ~  T  +  sum  of  all  ep(m)  ,£  ^  N 

m 

Eqs.  (6-73)  and  (6-74)  can  be  easily  expanded  when  N  is  small.  As  an 
example,  for  N  =  2 ,  we  have 

f/l  -  af  •  |/l  -  a2z 

*R  =  1  +3^2  exp  (-i^  +  d>2))  6xp  ('l(q>l  +  *2^ 

a}  exp  (-i <t>i )  +  a2  exp  (-i  (2(p]  +  <f2) 
rR  1  '+  a^  exp  (-"i"(01  +  42)) 

and  for  N  =  3  we  have 


(6-73) 


(6-74) 


-  a,  ^  1  -aj  j/l  -  a3  exp  (-♦(*,  ♦  *2 


*3> 


1  ♦ 


a,a?  exp  (-1(41,  ♦  «2))  ♦  a2a3  exp  (-1U2  ♦  $.,))  *  a,a3  exp  < - 1  ( i 3  *  2c 2  ‘  <3)) 


a,  exp  (-($,)  ♦  a2  exp  (-i(2»,  ♦  »2))  *  a3  exp  ( -1  ( 2c,  *  ♦  i3 ) )  *  ‘ji2a3  *  *2  * 

1  +  a,a2  exp  (-1(4,  *  c2 ) )  ♦  a?a3  exp  (-iU2  *  ;3))  ♦  a,a,  exp  ( - f ( c ,  ♦  2*2  ♦  i 3) ) 


It  should  be  pointed  out  that  Eqs.  (6-73)  and  (6-74)  are  analogous  to 
the  complete  Airy  summation  formulas  for  multilayer  dielectric  stacks 
[1]. 
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Note  that  Eqs.  (6-73)  and  (6-74)  can  be  combined  to  yield  the  fol¬ 
lowing  result 


( /7y^l  -  a j )  *  exp  (-i  )>  ^  <pj)  =  sum  of  all  op(m) ,  -  N 

j=1  (6-75) 


j=l 


2 

But  f)R,  the  right  going  reflectance,  equals  |rRj  ,  and  the  right  going 
transmittance,  TR,  must  equal  1  -  pR  because  of  energy  conservation. 
Thus,  Eq.  (6-75)  becomes 


N  __ 

(  R  —  )1/?2  (FT  |/ 1  -  a2  )  =  |sum  of  all  op(m) ,  (  <  N 

I  -  pp  *  "  J  ^ 

J=1  (6-76) 


We  will  return  to  Eqs.  (6-73),  (6-74)  and  (6-75)  when  we  investigate 
the  theory  of  holographic  grating  filter  design  in  the  following 
chapter. 

We  note  from  Eq.  (6-24)  that,  in  general. 


and 


r*  exp  (i 21 R  ) 

j  J 

exp  (i2TR  ) 
j 


(6-77) 


N 

Q  =  /  /  Q j  (6-78) 

j  =  l 


< 


f 
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Theorem  6-3  is  generalized  below.  We  start  by  making  the  following 
defini  tions . 


1 


Z. 

J 


A 


r*  exp  ( i 2Fr  ) 
exp  { i 2Tr  ) 

J  J 


(6-79) 


The  class  of  generalized  even  products  of  order  m,  designated 
GEP(m)  are  expressions  of  the  form 


b 


i 


m-1  4  j+r1 

b£  exp  [i  E  Z  (1'  ("1)J)‘R.] 

J  J 


for  j  an  odd  integer 
for  j  an  even  integer 


(6-80) 


where  m  is  a  positive  even  integer  and  •••  >  m  ’s  d  monotoni- 

cally  increasing  sequence  of  positive  integers.  An  element  of  this 
class  will  be  denoted  gep(m).  The  class  of  generalized  odd  products 
of  order  m,  designated  GOP(m) ,  are  expressions  of  the  form 
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m-1  £j+l_1 


bo  bo  ...  bj,  exp  [ 


EE  " - <- 


j-O  n*lj 


Uj+1)r  ] 

j 


r^  for  j  an  odd  integer 


!rs  Tor  j 
* 

r£  for  j 


j  an  even  integer 


£,  =  0 


r  =  o 

Ro 


(6-81) 


where  m  is  a  positive  odd  integer  and  ^ »  •••  *  a  nnonotoni- 

cally  increasing  sequence  of  positive  integers.  An  element  of  this 
class  will  be  denoted  by  gop(m). 

Using  the  definitions  above,  we  state  the  following  theorem 
which  is  a  generalization  of  Theorem  6-3. 

Theorem  6-4; 


11  h 


,(N)  V(N) 

V11  v12 


V(N)  V(N) 

*21  v22 


(6-82) 


where 


=  1  +  sum  of  allgep(m),  -  N 


(6-83) 


V^i  =  sum  of  allqop(m),  -  N 


(6-84) 
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_ 

12 


V<?}*  exp 


j=l 


(6-85) 


,(N) 

*22 


W(N)* 

V-j-j  exp 


t’2<ErR; 


j=i 


(6-86) 


Proof: 

Theorem  6-4  will  be  proved  by  mathematical  induction.  It  is  trivial 
to  show  that  Theorem  6-4  is  true  for  N  =  1.  Now  suppose  the  theorem  is 

true  for  N  =  N  ,  then  it  is  sufficient  to  show  that  it  is  also  true  for 
o 

N  =  N  +  1 .  We  wri te 
o 


N  +1 
o 


■  n  v  <n  v  2V, 

j=i  j=i  0 


(N  )  (N  ) 

/  0  v  0 
11  12 


(N)  (N  ) 

V  0  V  0 
21  22 


R 


exp  ( i  21' 

lN  +1  N  +1 

o  o 


n  e*P  (121'n  ) 

kN  +1  +1 

0  0 


There  fore , 


(N  )  (N  ) 
e  =  Vn  +  V1 2°  rR, 


(6-88) 


N  +1 
o 
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c  =  V 


(VD 

21 


N  +1 
o 

b  =  c*  exp  [i 2 (  rR  )  ] 

j-1  'j 


(6-93) 


(6-94) 


Ue  prove  Eq.  (6-92)  first.  Equation  (6-88)  can  be  rewritten  by  using 
Eqs .  (6-83),  (6-84),  and  (6-85)  as 
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(N>  (N0)* 

e  =  V11  +  V21  exp  ti2  2_j  "jl 

j  =  l 


N  +1 
o 


(1  +  sum  of  all  gep(m),  -  N) 


+  (sum  of  all  gop(m)  ,  -  N)* 


exp 


-r.J' 


J  N  +1 

J-I  0 


(6-96) 


Consider  now  the  arbitrary  generalized  odd  product  of  order  m  given 
below  Um  -  Nq) 


gop(m)  -  b.  b{  ...  bf  exp 
1  2  m 


m-1  "j+1'1 


j=0  n=ij 


(6-97) 


Then 


gop(m)*  exp  [i2(^ 


j'J  rR, 


N  +1 
o 


m  \j+l 

v*,  ••• b* b^*,  £  Z  ('-(-i>j,rR]  (6-9s) 

'  2  m 
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+  (1  +  sum  of  allgep(m),  £  -  N  )* 
r  m  o 


v.* 

•  exp  [12  rR  )]  rR 

Ki  RN  +1 


TT  j  V 


(6-99) 


Consider  now  the  arbitrary  even  product  of  order  m  given  below  (£  -  N  ) 

m  o 

m-1  £j+l"1 

gep(m)  =  b£  b£  ...  b£  exp  [V'  V'  (1  -  (-l)j)rR  ]  (6-100) 

12  m  Tt. 

J 


gep(m)*  exp  [i2(^  rR  )]  rR 


b£  b*0 


m  ^j+r1 

b?  bt . exp  1 Z)  y,  ° 

m  m+1  *r~r 


j=0  n-tj 


(6-101 ) 


where 


£  ■,  =  N  +  1  and  b„  =  r„ 

m+1  0  Vl  RN  +1 

o 


Note  that  the  right  hand  side  of  Eq.  (6-101)  is  just  an  element  of 
GOP(m),  £m  =  +  1 •  Furthermore,  it  is  easy  to  see  that  all  elements 

of  gop(m),  £m  =  Nq  +  1  and  only  those  elements  can  be  generated  by  mul¬ 


tiplying  1  or  some  gop(m),  ^  -  Nq  by 
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Finally,  Eq.  (6-95)  will  be  proven.  Eq.  (6-91)  can  be  rewritten  using 
Eqs.  (6-85)  and  (6-86)  as 


(N  )*  u  y 

-  V0  exp  C12<2J  rR  >]  r*  exp  <12r  > 


n  +1 

j=l  J  o 


(M  )*  ^ 

V,!0  exp  [12  (^J  rR_)]  exp  (12r, 


(N  )*  (NJ*  ^ 

=  ^V12  rR„  +  V11  ^  exp  rR,)-1 


(6-103) 


Now  from  Eqs.  (6-88)  and  (6-103) 


' J 

d  =  e*  exp  [12(^  rR  )] 


The  proof  is  complete. 


Using  the  above  theorem  it  immediately  follows  that 


1  1  +  sum  of  all  gep(m),  -  N  (sum  of  all  gop(m),  -  N)*  exp  ( i 2  r  j 

N  1  J  =  1  J 

o  =  (  n  t- '  )  | 

M  R,  .  N 

J  ;  sum  of  all  gop(m),  n  -  N  ( 1  +  sum  of  all  gep(m)  ,  -  N)  exp  ( i  2  X)  ) 

I  m  j  =  l  J 


(6-104) 
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and  so,  by  Eqs.  (6-1)  and  (6-104) 


= 


N 

JT  tR  . 

J  =  1 


R  1  +  sum  of  all  gep(m),  <  N 


(6-105) 


sum  of  all  gop(m) ,  2m  -  N 
R  T  +  sum  of  all  gep(m),  £  <  N 


Combining  Eqs.  (6-105)  and  (6-106)  yields 


r  N 

1  n 

t~(J7  tR  )  =  sum  of  a11  g°p(Ri) .  £m  -  n 
r  j-i  j 


But  ! tR  | 2  =  1 

j 


(6-106) 


(6-107) 


,  2 

| rR  [  because  of  energy  conservation.  Furthermore,  .  R, 

j  J  2 

the  right  going  reflectance  equals  j  rRj  and  the  right  going  transmittance, 
Tp,  must  equal  1  -  pR  because  of  energy  conservation.  Thus,  Eq.  (6-107) 
becomes 


N  r—— 

(  1  _-R—  ) 1/2  J7  '  lrR.i^  =  ^sum  of  a11  V  ■ 

R  i-l  j 

J  1  (6-108) 

At  this  point  we  wish  to  prove  two  additional  theorems  which  will  be 
used  in  the  following  chapter. 
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Jim 


Theorem  6-5 : 


u 

sum  of  all  op(l)  =  a(eo) J  n'(z)  exp  (-i(^r  nQ  cos  0o)z)  dz 


=  a(0  )  grtn'(z)  rect(f)}  _  2  no  cos  9o 
0  U  TZ  A 


(6-109) 


where  { • } *  denotes  Fourier  transform  evaluated  at  f 

2  N  , 

D  =  the  grating  thickness  =  2_,  T~ 

j  =  l  Tj 


rect  ( Z)  M  1 •  0  -  z  -  D 
"  /  0,  otherwise 


denotes  differentiation. 


Proof: 


u 

f  n'(z)  exp  (-i(xn0  cos  e0)z)  dz 


l/f 

r 


^n'(z)  exp  nQ  cos  °0)z)  dz  + 


m+1 

E  l/f i 

N-l  -j*l  J 


.,4" 


^  I  n'(z)  exp  (-i(xn0  cos  »*0)z)  dz 

m=l  Ja 

E  i/fj 

j=i  J 
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But 


m+1 

E  i/f, 

j=i  J 


/ 


4ti 

n'(z)  exp(-i(—  n  cos  6  )z)  dz 

AO  0 


m 

E  i/f. 


becomes  with  a  change  of  variable 


exp 


-i(Tno cos  eo>  Et;1 


1/f, 


/ 


m+1 


J  Q 


n'(z  +  Ef1 

j  =  l  J 


•  exp  (-i(-^nQ  cos  6q)z)  dz 


Furthermore , 


n(z  +^77)  = 


nm+l(z)  =  Am+1  Sin  (2*fm+lz) 


j  =  l 


Integrals  of  the  form 


1/f, 


/ 


m+1 


4  7!" 

n'm+l(z)  exp  ('1(T  no  C0S  Go)z)  dZ 
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have  already  been  evaluated  in  Chapter  5  (Eqs.  ( 5-4) - (5-7) ) .  It  was 
shown  in  Chapter  5  that 


n'm+l(z)  exp  (-1  (—  nQ  cos  eQ)z)  dz 


m+1  ,  . 

(”Vi> 


Also  note  that 


exp  (-i(ifnQ  cos  0Q)  (^77))  = 


j=l  J 


exp  (-i2(y;  ) ) 

j=1 


Therefore 


v/" 


(z)  exp  (-i(xno  cos  eo^z^  dz  = 


N  m-1 

a1  exp  (-i 0-j )  exp  (-i(2(^  +  <j>m)) 

m=2  j=l 

By  Equation  (6-49),  op(l)  is  an  expression  of  the  form 

exp  (-i  (2(^-1  +  <j>2  +  ...  +  4*^^  )  +  4^)) 


And  so 
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VeRIM. 


N  m-1 

sum  of  all  op(  1 )  =  a]  exp  (-i^)  exp  (-i(2(^  ♦  ;m)) 

m=2  j=l 

D 

»  a(dQ)  j  n 1  (z)  exp  (-i(4^  nQ  cos  ejz)  dz 
0 

GO 

=  a(0Q)  j  n 1  (z)  rect  (-^)  exp  (-i(^f  nQ  cos  0O)Z)  dz 

-CO 


=  o(0Q)  ^{n1  (z)  rect  (^)l  f 


2no 

—  cos 


1 1 

o 


The  proof  is  complete. 

Theorem  6-6: 

D 

gopd)  -  a(e0)  J  n'(z)  exp  (-i(^  %  cos  6Q)z)  dz  (6-110) 
0 

for  small  | rR .  | ,  j=l ,  . . . ,  N 

*J 

where  Ai  =  the  thickness  of  the  jth  slab 

«J 

N 

D  =  the  grating  thickness  =  £  a. 

j=l  J 

'  denotes  di  fferentation . 
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Proof: 

We  showed  in  Chapter  5  (Eqs.  (5-3)  and  (5-4)  with  =  0)  that 

f3  j"1 

rR.  '  'a{eo^  /  n'(z  +  XI  Am^  exp  n0  cos  Vz)  dz 

J  n  m=l  (6-111) 


for  small  i rRj ! •  Similarly,  for  the  left  going  reflection  coefficient 
of  the  jth  slab,  r|_ . 

AJ  j 

rL-  =  a(0o)  /  Am)  "  z^  exp  n0  cos  ‘Vz>  dz 

J  i  m=l 


for  small  |r(.|.  Now  with  the  substitution  t  =  A,  -  z,  the  above  equa- 
J  J 

tion  becomes 


rL.  "  I  -a (ej 


P  j-' 

>o>J 

n  m=l 


Am)  exp  (i  (^  nQ  cos  »0)t)  dt 


exp  ( -i  (~  nQ  cos  fcQ)Aj)  (6-113) 


for  small  j  rj_ . j  .  Therefore,  by  Eqs.  (6-111)  and  (6-113) 

J 

r  ■  -rR*  exp  (-i (~  no  cos  e0);.j) 

J  J 
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for  small  j  r  |_ .  j  and  ,  r  r  -  i .  But  by  Eq.  (6-14) 

J  J 

r,  =  -r*  exp  ( i 2 r R  ) 

j  j  j 


and  so 


2rR.  *  -  x  no(cos  6o)Aj 


(6-114) 


for  small  jrR.|  and  |  r  .  j.  Note  that 
J  J 


U 

I90)  j 


a(6  )  /  n 1  (z)  exp  (-i(y-  nQ  cos  GQ)z)  dz 


oi(e)  /  n'(z)  exp  (-Ktt  n0  cos  0o)z)  dz 


J 

y  a. 

m=l  J 


+  u(e) 


S/ 

J"2  j-i 

£ 


n'(z)  exp  (-i(  —  nQ  cos  ' Q) z )  dz 


=  a(0) 


/- 


(z)  exp  ( -i (—  nQ  cos  hq)z)  dz 


(Equation  continued  on  next  page.) 
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N 

+  ct(e)  ^ 
j=i 


•  exp  (-i (j-  nQ  cos  0Q)t)  dt 


j-1 

(We  made  the  substitution  t  =  z  -  ^  ^  ). 

m=i 


Now  for  | rR j I  small  j  =  1,  ...  ,  N,  Eq.  (6-115)  becomes 


D 

ct(6°)  J  n'  (z)  exp  (-i(y-  nQ 
0 


cos  t3Q)z)  dz  - 


r 


R 


1 


i  ■  i  4"n 
(-1  ( —  n 

A  o 


cos 


But  we  have  already  shown  that 


2TD  =  "V  n  (cos  0  )A . 
R •  X  ov  o'  j 


for  | rRj  i  and  |  rLj 


smal 1  ,  and  so 


DIVISION 


(6-115) 


(6-116) 


(6-114) 
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U 

v/"' 


,  4tt 


a(0j  I  n '  (z)  exp  (-i(—  nQ  cos  bQ)z)  dz 


N 


j-1 


-  rRl  +  _  .. 

1  j=2  J 


Z!rR  Cxp  1  12  2  rR> 


m=0 


(6-117) 


for  i  r^,  [  and  (  rL  j  |  smal  1  ,  j  -  1  ,  . . .  ,  N . 

By  Eq.  (6-10),  |  rLj.  j  =  j  rR^  1  and  so  j  rR  ■  ]  small  implies  that  rL. 
is  also  small.  But  by  Eq.  (6-81)  gop(l)  can  be  written  as 


j-1 


g°P{1)  =  r  +  rR-  6XP  (i2S  PR  } 


1  “  J 


m=l 


and  so  the  proof  is  complete.  We  make  three  remarks.  First, 


1  i  m  r  D  =  0 . 

y°  rj 


Secondly,  gop(l)  asymptotically  approaches 


D 

>*(0o)  j  n '  (z ) 

0 


exp 


COS 


dz 


as  the  A  . 

J 


s  approach  zero. 
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,  /  • 


/  n’(z)  exp  (-i(— nQ  cos  Qq)z)  dz 


rect  (^)}f  =  _Q 


2n  cos  5 


and  so 


gop(l)  -  ot(6  )  9T{ n'(z)  rect  (|)}f  =  2no_C°S  9_o 

Z  A 
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Notes  to  Chapter  6 

[1]  Z.  Knittl,  Optics  of  Thin  Films,  John  Wiley  &  Sons,  New 
York,  1976 

[2]  *  denotes  conjugation. 

[3]  See  Equation  (5-3)  with  r.  set  to  zero. 
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7 

THEORY  OF  HOLOGRAPHIC  GRATING  DESIGN 

The  ultimate  goal  of  any  applied  science  is  the  successful  applica¬ 
tion  of  fundamental  principles  to  the  design  of  useful  devices.  In  al¬ 
most  all  cases,  the  development  of  techniques  for  analyzing  observed 
phenomena  predates  the  development  of  techniques  for  designing  devices 
based  on  those  phenomena.  Furthermore,  design  capability  is  generally 
far  more  difficult  to  attain  than  analysis  capability.  Design  is  an 
engineering  problem,  and  as  such,  is  often  as  much  an  art  as  a  science. 

The  design  of  holographic  gratings  is  in  many  respects  similar  to  the 
design  of  multilayer  dielectric  devices.  A  multilayer  dielectric  device 
is  a  stack  of  thin  homogeneous  dielectric  layers  (see  Figure  3-1  ).  The 
study  of  such  structures  is  often  referred  to  as  the  optics  of  thin  films. 

We  digress  briefly  to  discuss  the  history  of  the  optics  of  thin  films. 

The  interference  colors  produced  by  soap  bubbles  and  by  oils  on 
water  were  described  as  early  as  the  17th  century.  Isaac  Newton  was 
unable  to  explain  this  phenomena  with  his  corpuscular  theory  of  light. 

It  was  not  until  the  development  of  physical  optics  by  Young  and  Fresnel 
that  the  colors  of  thin  films  was  discovered  to  be  caused  by  the  inter¬ 
ference  of  light.  The  Fabry-Perot  interferometer  was  the  first  prac¬ 
tical  device  based  on  the  principle  of  interference. 

In  1917,  the  spectral  colors  of  the  covers  of  some  beetles  observed 
in  reflection  was  explained  by  Lord  Rayleigh  [1].  The  covers  of  the  beetles 
was  essentially  a  laminar  structure  analogous  to  the  "high-low,  high-low" 
stacks  of  multilayer  dielectric  interference  filters.  At  that  time,  how¬ 
ever,  no  techniques  for  fabricating  multilayer  dielectric  structures  were 
available  and  Lord  Rayleigh's  paper  received  little  attention.  The  first 
real  appearance  of  thin  film  optics  can  be  traced  to  1904  when  a  patent 
was  granted  for  an  antireflection  technique  [2]. 
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Practical  developments  of  thin  film  optics  didn't  occur  until  the 
1930' s  when  they  sprouted  from  advances  in  vacuum  technology.  It  be¬ 
came  possible  then  to  vacuum  deposit  thin  transparent  materials  of  spec¬ 
ified  thicknesses  onto  substrates.  Thus  was  the  atmosphere  in  the  period 
of  1937  to  1947  during  which  the  first  theories  for  anti refl ecting  coatings 
and  monochromatic  interference  filters  were  developed.  It  was  not  until 
1950  however,  that  a  general  treatment  of  electromagnetic  wave  propagation 
through  stratified  media  was  published  [3].  This  treatment  was  later  sim¬ 
plified  by  the  introduction  of  matrix  techniques  (see  Chapter  3). 

In  this  early  period  the  primary  focus  of  research  was  mainly  in  the 
area  of  analysis,  closed  form  expressions  for  reflectivities  of  multi¬ 
layer  dielectric  stacks  were  conceptually  simple  but  had  large  numbers  of 
parameters  (the  refractive  index  of  each  layer,  the  thickness  of  each 
layer,  and  the  number  of  layers).  Determining  these  parameters  to  yield 
desired  reflectivity  characteristics  is  an  extremely  difficult  design  prob¬ 
lem.  Great  progress  has  been  made  in  this  area  [4,  5,  6].  Initially,  de¬ 
sign  theory  was  limited  to  very  special  problems  such  as  anti refl ection 
coatings.  Later,  however,  theoretical  work  was  undertaken  to  develop 
systematic  design  procedures  applicable  to  wider  classes  of  problems. 
Although  a  tremendous  amount  of  effort  has  been  expended  in  an  effort  to 
develop  these  procedures  and  significant  progress  has  been  made,  the  goal 
has  not  been  completely  reached.  The  design  of  mul ti 1 ayer  dielectric 
stacks  is  to  an  extent  still  an  art  as  is  the  design  of  optical  lens 
systems.  Approximate  techniques  are  often  used  to  obtain  a  first  cut 
solution  to  design  a  problem.  The  design  is  then  analyzed  using  simple, 
exact  techniques,  and  refinements  are  made  as  needed  to  the  original  de¬ 
sign.  The  use  of  a  digital  computer  greatly  facilitates  this  iterative 
procedure . 

The  theoretical  problem  of  design  consists  of  specifying  the  param¬ 
eters  in  the  expressions  for  reflectivity  so  as  to  yield  desired  device 
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characteristics.  In  order  to  physically  realize  the  device  however, 
these  parameters  must  satisfy  certain  technological  constraints.  For 
multilayer  dielectric  stacks,  these  constrains  limit  the  maximum  number 
of  layers  in  the  stack,  and  also  the  materials  which  can  be  used  to  form 
dielectric  layers. 

We  have  developed  a  theory  for  analyzing  holographic  grating  filters. 
An  alternate,  currently  available  analysis  technique  based  on  multilayer 
dielectric  structures  was  already  known,  though  not  used,  prior  to  our 
work  (see  Chapter  4).  The  original  motivation  for  our  research  was  to 
develop  an  analysis  technique  which  would  be  computationally  more  effi¬ 
cient.  We  were  successful  in  this  endeavor.  Furthermore,  we  will  show 
that  our  analysis  technique  can  be  used  to  gain  considerable  insight  into 
the  design  problem.  Equations  (6-74),  and  (6-107)  of  Chapter  6  are  key 
results  which  we'll  exploit  in  order  to  develop  design  techniques. 

The  design  parameters  in  holographic  grating  filters  are  the  A,'s 

J 

and  f . ‘ s .  As  in  the  design  of  multilayer  dielectric  stacks,  these  param- 

J 

eters  must  satisfy  certain  technological  constraints.  For  the  holographic 
case,  these  constraints  place  restrictions  on  the  types  of  refractive  index 
profiles  which  can  be  holographically  recorded  in  available  materials.  It 
has  already  been  mentioned  in  earlier  chapters  that  these  refractive  index 
profiles  are  quasi-sinusoidal  (i.e.,  4  urn'1  <  f  <  6  urn-1)  of  small  index 
modulation  (i.e.,  |A  .  |  ^  0.03).  Furthermore,  the  element  thickness  (i.e., 

N 

/  T~)  is  limitedto  athickness  of  less  than  approximately  30  pm.  Refrac- 

j-1  j 

tive  index  profiles  which  do  not  lie  in  this  class  are  not  valid  solutions 
to  the  design  problem.  Discussions  of  actual  procedures  for  fabricating 
a  holographic  filter  will  be  deferred  to  Chapter  8.  The  design  problem  can 
now  be  succinctly  stated  as  follows: 

Given  tp  and  r^  as  a  function  of  6  ,  \  and  polarization,  deter¬ 
mine  the  A.’s  and  f.'s  which  yield  those  characteristics.  These  A  's 

J  J  j 

and  f j ' s  must  also  satisfy  the  constraints 
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|  A  .  |  <  0.03 


4  pm 


<  f.  < 

-  J  - 


6  pm 


emulsion  thickness  = 


30  pm 


A  solution  for  a  design  problem  may  not  always  exist. 

Equations  (6-73),  (6-74) , (6-1 05)  and  (6-106)  of  Chapter  6  are  anal¬ 
ogous  to  the  Airy  summation  formulas  for  multilayer  dielectric  stacks  [4], 
There  is  good  reason  to  believe  that  a  number  of  techniques  used  in  the 
design  of  multilayer  dielectric  stacks  may,  with  modification,  be  appli¬ 
cable  to  the  holographic  grating  design  problem.  This  area  has  not  been 
explored  as  of  yet,  but  it  appears  promising.  Our  discussion  will  be 
limited  to  a  simple,  approximate  holographic  design  technique. 

Before  we  begin,  let  us  discuss  the  motivation  behind  this  work. 

There  is  a  market  for  multilayer  dielectric  devices  as  attested  to  by 
the  wide  spread  use  as  antirefl ection  coatings,  dielectric  mirrors, 
polarizing  beam  splitters,  and  interference  filters.  All  of  these  de¬ 
vices  are  currently  fabricated  using  evaporated  thin  film  technology. 

The  possibility  exists  for  fabricating  holographically ,  devices  similar 
to  these.  Also,  one  would  hope  to  produce,  holographically,  devices 
which  cannot  be  realized  using  multilayer  dielectric  structures.  An  ex¬ 
ample  of  such  a  device  is  a  head-up  display  unit  for  aircraft  [7].  The 
fabrication  of  a  holographic  element  is  conceptionally  simpler  than  the 
fabrication  of  multilayer  dielectric  structures  since  it  is  basically  a 
photographic  process.  At  this  point,  it  is  only  fair  to  say  that  the 
practical  realization  of  sophisticated  holograhic  grating  devices  is  un¬ 
certain.  Only  time  will  tell. 

In  almost  instances  where  a  holographic  optical  element  is  used,  the 
element  is  embedded  in  a  medium  having  a  nominal  refractive  index  of  one. 
Such  a  medium  is  of  course  air.  We  will,  however,  restrict  our  discussion 
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to  the  case  where  the  element  is  imbedded  in  an  index  matching  medium 
assumed  to  have  a  nominal  refractive  index  of  1.5.  The  reason  behind 
this  decision  is  mainly  pedagogical,  though  it  does  simplify  the  develop¬ 
ment  of  design  techniques.  In  any  case,  the  analysis  techniques  developed 
so  far  place  no  restrictions  on  the  imbedding  medium.  As  a  last  point, 
we  note  that  Eqs .  (6-73),  (6-74),  (6-105)  and  (6-106)  of  Chapter  6  are 
strictly  applicable  when  the  imbedding  medium  is  of  refractive  index  1.5. 
They  can,  however,  be  extended  to  the  more  general  case. 

Equations  (6-74),  (6-107)  and  (6-108)  of  Chapter  6  are  repeated  below 
for  reference. 

sum  of  al  1  ou(m)  ,  <  N 

r  =  - - -  (7.-1) 

1  +  sum  of  al  1  ep(m)  ,2.  <  N 

m  - 


tR  is  the  right  going  transmission  coefficient  of  the  grating,  rR  is 
the  right  going  reflection  coefficient  of  the  grating,  and  Pd  is  the 
right  goint  reflectance  of  the  grating.  rR  and  tR  ar  the  right  going 

J  J  +  h 

reflection  and  transmission  coefficients  repectively  of  the  j  slab  of 

the  grating.  We  note  that  in  Eqs.  (7-2)  and  (7-3)  the  thickenss  of  this 
t  h 

j  slab  is  arbitrary  and  therefore  not  necessarily  an  integral  number 
of  periods  of  the  refractive  index  profile.  In  Eq.  (7-1),  however,  the 
j  slab  has  a  thickness  equal  to  the  period  of  the  jtn  cycle  of  the  re¬ 
fractive  index  profile. 
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Let  us  examine  Eq.  (7-1)  first.  We  note  from  Eqs.  (6-47)  and  (6-49) 
of  Chapter  6  that  both  op(m)  and  ep(m),  m  >  2  are  sums  of  terms  which  are 
of  order  greater  than  or  equal  to  two  in  the  aVs  (i.e.,  each  term  con¬ 
tains  a  product  of  at  least  two  a^.'s).  As  can  be  seen  from  Eq.  (5-24)  of 
Chapter  5,  the  a^'s  are  rather  small  quantities.  Thus,  as  a  first  approxi¬ 
mation  all  ep(m)  and  op(m),  m  >  2  will  be  neglected  in  Eq.  (7-1).  Then 
this  equation  becomes  to  a  first  order  approximation 

r„  =  sum  of  al 1  op(l )  (7-4) 

(note  ep(l )  s  0) 


From  Theorem  6-5,  Eq.  (7-4)  is  rewritten  as 


rR  * 


n' (z)  rect 


cos  e  g 


(7-5) 


Equation  (7-5)  is  a  well  known  result.  Recall  that  it  was  rigorously  de¬ 
rived  from  Maxwell's  equations  in  Chapter  4  for  the  case  of  low  reflectivity 
(i.e.,  |rR|  small).  A  much  better  approximation  (i.e.,  one  which  is  less 
restrictive  on  the  magnitude  of  r  )  can  be  obtained  from  Eqs.  (7-2)  and 
(7-3)  by  neglecting  all  terms  which  involve  products  of  three  or  more 

rR  's  when  the  rR  's  u.c  small.  We  note  from  Eqs.  (6-80)  and  (6-81  of 
i  i 

Chapter  6  that  gep(m),  m  .  3  and  gop(m),  m  1  3  are  sums  of  terms  which  ea>  h 

term  contains  a  product  of  at  least  three  rD  ’s.  Thus,  as  a  second  order 

Rj 

approximation  all  gop(m),  m  :  3  and  gep(m),  m  4  3  will  be  neglected  in 
Eqs.  (7-2)  and  (7-3).  Then  these  equations  become  to  a  second  order  ap¬ 
proximation. 


(see  comments  following  Eq.  (7-^. 


Note  thatgop(2)  -0  and  also  that 
N  N 


But  from  the  Eq.  (6-114)  of  Chapter  6,  we  have 
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From  Theorem  6-6,  Eqs .  (7-6)  and  (7-7)  can  be  rewritten  as 


Equations  (7-4)- (7-9)  above  are  analogous  to  results  derived  in  multi¬ 
layer  dielectric  stack  theory  [4].  In  that  theory,  these  results  play 
an  important  role  in  the  area  of  design.  Equations  (7-5),  (7-6)  and 
(7-7)  will  serve  as  the  basis  of  our  design  procedure.  The  design  stat- 
egy  using  a  first  order  approximation  is  as  follows: 

We  are  given  rR  vs.  6^  for  fixed  A  and  polarization  or  rR  vs.  a 
for  fixed  oQ  and  polarization.  Write  rR/a(e0)  as  a  function  of 

2n0  rR 

f,  =  — —  cos  0n.  Then  take  the  inverse  Fourier  transform  of  —r, — v  To 
Z  x  o  a  ( e  Q ) 

obtain  n'(z)  rect  (-^)  and  consequently  n'(z),  0  <  z  <  D.  Finally, 

n(z),  0  <  z  <  D  can  be  obtained  from  n'(z)  by  integration.  For  a  design 

using  a  second  order  approximation,  the  above  procedure  is  modified. 

We  are  given  rR  and  tR  vs.  6g  for  fixed  x  and  polarization  or  rR 

and  tR  vs.  A  for  fixed  0g  and  polarization.  Make  an  initial  estimate 

of  n(z)  (and  consequently  the  A.'s  and  f  .'s).  Use  our  analysis  technique 

J  J  ?nn 

to  compute  the  r  's  as  a  function  of  f  =  -  cos  oQ.  Write 
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to  obtain  a  new  estimate  of  n'(z)  rect  (-^) .  n(z),  0  <  :  <  D  and  conse¬ 

quently  the  A,’s  and  f.'s  can  then  be  obtained  by  integration.  Repeat 
J  J 

the  above  process  until  n(z)  does  not  change  significantly  from  one  iter¬ 
ation  to  the  next.  There  are  a  number  of  important  points  which  were 
not  explicitly  mentioned  in  the  above  discussion.  Let  us  first  reconsider 
the  design  technique  which  was  based  on  the  first  order  approximation. 

Note  that 

(1)  It  is  valid  only  in  situations  where  the  maximum  value  of 

2 

| rR |  is  small  . 

2  2 

(2)  Energy  conservation  requires  that  |rD|  +  |tD|  =  1.  Thus 

2  2  K  i\ 

|r^|  and  |t^j  cannot  be  specified  independently. 

(3)  Inversion  of  the  Fourier  transform  requires  knowledge  of 

rR  2n 

— n-Ti  expressed  as  a  function  of  f  =  — —  cos  on,  for  all 

u  i  e  q  ;  z  a  u  — 

fz,  -  <•  f  <  x  and  Oq,  however,  lie  in  the  intervals 

4  pm  <  A  <•  7  pm  and  0°  <  o  42°  (42°  is  approximately 
the  critical  angle  of  gelatin  when  imbedded  in  air).  Thus, 

7~-  t  is  at  most  known  for  f  in  the  interval  cos  4;- 

alon)  — — -  z  .7  ,,m 

3.2  . -m~ 1  to  cos  f)  5  7.3  .in'  .  Values  must  be  assigned  ti 

.4  .  .m 

rR 

-77—r  for  all  f  outside  of  the  known  interval.  Although  this 
1  0; 
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assignment  will  not  effect  the  rR  vs.  e q  (or  x)  character¬ 
istic  over  that  interval,  it  will  affect  the  n(z)  obtained 


from  the  design  procedure.  Note  that  y 
at  -f  if  n'(z)  is  to  be  real  valued. 


at  f  must  equal 


(4)  It  is  likely  that  the  inverse  Fourier  transform  of 


R 

Oi(0q) 


wi  1 1 


not  vanish  outside  of  some  finite  interval  on  the  positive 

z  axis.  Our  theory,  however,  requires  it  to  be  of  the  form 

n‘(z)  rect  (^).  Therefore,  it  will  be  necessary  to  truncate 

this  inverse  transform  to  a  finite  size  and  to  shift  it  along 

the  z  axis  so  that  it  will  vanish  outside  of  some  interval 

[0,  D],  D  >  0.  This  operation  will  yield  a  r  vs.  0  (or  X) 

“  K  U 

characteristic  which  differs  from  the  specified  one.  The 
effect  of  the  shifting  operation  is  to  alter  the  rR  vs. 

(or  A)  characteristic  by  a  multiplicative  phase  factor  exp 

r  V  ^ 


-1 


2n 


0 


j-  cos  Sq  )  s^j,  where  Sq  is  the  distance  along  the  z  axis 
by  which  the  transform  was  shifted  (this  result  follows  im¬ 
mediately  from  the  theory  of  Fourier  transforms).  Thus,  the 
phase  characteristic  of  rD  vs.  o  (or  1 )  can  be  designed  only 
up  to  this  multiplicative  phase  factor.  The  effect  of  trun¬ 
cating  the  inverse  transform  to  a  finite  extent  D  (note  D  is 
the  element  thickness)  is  to  alter  the  specified 

rR 

«/o"T  vs‘  ^ z  characteristic  by  convolving  it  with  the  function 

sine  (Dfz).  This  conclusion  also  follows  from  Fourier  trans¬ 
form  theory.  The  larger  the  value  of  D,  the  less  effect  this 
truncation  will  have.  D,  however,  is  the  element  thickness. 

We  have  already  noted  that  this  is  limited  to  about  30  <<m  by 
current  fabrication  technology. 
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(5)  There  is  no  guarantee  that  the  solution  obtained  will  be  in 
the  class  of  achievable  refractive  index  profiles.  If  it 
does  not  lie  in  this  class,  then  it  isn't  a  valid  solution, 
fr. 

_ _  anH  n 1  f  7  ^  rorf  I 

1 D' 


R  z 

and  n'(z)  rect  (•=.)  are  approximate  Fourier  transform 


Ot  ( 0  Q ) 


pairs.  Therefore 


R 


a(6  q) 


,  expressed  as  a  function  of  f  ,  should 

-1 


(6) 


be  a  narrowband  signal  with  center  frequency  between  4  pm  and 
6  um" ^ .  This  follows  from  the  fact  that  n'(z)  is  a  quasi  - 
sinusoidal  function  of  nominal  frequency  4  pm~ ^  to  6  um  \ 

2 

In  many  cases  only  the  reflectance,  i.e.,  | r^ |  ,  is  specified. 
Therefore  only  the  modulas  of  the  Fourier  transform  of  n'(z) 

rect  (^)  is  known.  Phase  information  must  be  assigned  to 
rR  in  order  to  obtain  n(z).  This  information  should  be  chosen 
so  as  to  maximize  the  probability  of  achieving  a  realizable 
refractive  index  profile. 

For  the  design  technique  which  was  based  on  the  second  order  approxima¬ 
tion  we  note  that 


R1 


Thus 


(1)  It  is  valid  in  situations  where  the  maximum  value  of  !rn^  is 
fairly  large. 

p  n 

(2)  Energy  conservation  requires  that  |rn;  +  |tn|  =  1. 

2  2  R  1  R 

| rp |  and  |tR|  cannot  be  specified  independently. 

(3)  It  is  a  design  procedure  for  rR/tR  rather  than  for  rR  or  tR 
alone. 

(4)  Inversion  of  the  Fourier  transform  requires  knowledge  of 


R 


N  _ 

II  VTT 


j=l 


R  . 
J 


exp 


\ 


cos 


expressed  as  function 
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0 

of  f  =  - cos  0n,  for  all  f  ,  <  f  <  ■».  As  before, 

2  x0  Z 

this  information  will  only  be  known  over  some  finite  interval 

and  must  be  assigned  values  elsewhere.  Note  that 


at  -  fz  if  n'(z)  is  to  be  real  valued. 

(5)  It  is  likely  that  the  inverse  Fourier  transform  of 


will  not  vanish  outside  of  some  finite  interval  on  the  posi¬ 
tive  z  azis.  Our  theory,  however,  requires  it  to  be  of  the 
form  n’(z)  rect  (--).  Therefore,  it  will  be  necessary  to  trun¬ 
cate  this  inverse  transform  to  a  finite  size  and  to  shift  it 
along  the  z  axis  so  that  it  will  vanish  outside  of  some  inter¬ 
val  [0,  0],  D  >  0.  This  operation  will  yield  a  r  /tD  vs. 

0q  (or  \)  characteristic  which  differs  from  the  one  specified. 

The  effect  of  the  shifting  operation  is  to  alter  the  r  /tD 

R  R 

characteristic  by  a  multiplicative  phase  factor 
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0 

exp  [-i(-^—  cos  0 q ) s q ] ,  where  Sg  is  the  distance  along  the  z 

axis  by  which  the  transform  was  shifted  (this  result  follows 
immediately  from  the  theory  of  Fourier  transforms).  Thus,  the 
phase  characteristic  of  r  /tD  vs.  en  (or  x)  can  be  designed 
only  up  to  this  multiplicative  phase  factor.  The  effect  of 
truncating  the  inverse  transform  to  a  finite  extend  D  (note  D 
is  the  element  thickness)  does  not  have  a  simple  analytic  de¬ 
scription.  The  larger  the  value  of  D  of  course,  the  less 
effect  this  truncation  will  have. 

(6)  There  is  no  guarantee  that  the  solution  obtained  will  be  in 
the  class  of  achievable  refractive  index  profiles.  If  it 
does  not  lie  in  the  class  then  it  isn't  a  valid  solution. 

2 

(7)  In  many  cases  only  the  reflectance,  i.e.,  |rR(  ,  is  specified. 
Therefore  only  the  modulus  of  the  Fourier  transform  of 

n'(z)  rect  (-)  is  known.  Phase  information  must  be  assigned 
to  rR/tR  in  order  to  obtain  n(z).  This  information  should 
be  chosen  so  as  to  maximize  the  probability  of  achieving 
realizable  refractive  index  profiles. 

(8)  The  thickness  of  the  layers  (i.e.,  the  j ' s )  must  be  chosen. 
Experience  indicates  that  the  quality  of  the  approximation 

sum  of  all  gop(m),  <  N  -  go p ( 1 ) 

increases  as  the  maximum  values  of  rR .  c  j  =  1 . N 

increase.  On  the  other  hand,  the  quality  of  the  appro* i- 
ma  t  i  on 


gop(  1 ) 


(>3  )C1 


n* (z)  rect  (g) 


2n, 


cos 
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decreases  as  the  maximum  values 


Furthermore, 


the  maximum  values  of 


N  increase. 
. .  .  ,  N  are 


directly  proportional  to  the  a^.'s, 
be  made  in  choosing  these  A  • ' s . 


thus  some  compromise  must 


We  have  already  pointed  out  that  Eqs .  (7-1),  (7-2)  and  (7-3)  are 
analogous  to  key  formulas  in  the  theory  of  multilayer  dielectric  stacks. 

In  that  theory  the  a ^ ' s  are  roughly  equivalent  to  the  reflectivities  at 
the  boundaries  between  pairs  of  homogeneous  dielectric  layers,  and  the 
i+> j ' s  are  essentially  equivalent  to  the  optical  thicknesses  of  the  homo¬ 
geneous  dielectric  layers.  These  equivalent  a ^ ' s  are  functions  of  the 
refractive  indicies  of  the  layers,  the  incident  angle  and  the  polariza¬ 
tion,  but  not  the  wavelength.  The  equivalent  c j ’ s  are  functions  of  the 
incident  angle,  the  wavelength,  the  thicknesses  of  the  layers,  and  the 
refractive  indicies  of  the  layers.  In  most  multilayer  dielectric  design 
techniques  the  reflectivity  or  the  transmittance  is  specified  as  a  func¬ 
tion  of  wavelength  for  some  fixed  angle  of  incidence.  Note  that  phase 
information  is  usually  not  specified,  and  consequently  one  does  not  exercise 
any  control  over  the  resulting  phase  characteristics.  Furthermore  since 
the  angle  of  incidence  is  fixed,  the  a^'s  will  be  unknown  constants. 

The  multilayer  dielectric  stack  is  generally  assumed  to  be  tuned,  i.e., 
the  ‘ ' s  are  all  integer  multiples  of  one  another.  The  expression  for 
; rR |  can  then  be  written  as  a  rational  function  of  z  by  making  the 

substitution  z  =  exp  (-if).  The  desired  characteristic  for  r  6  is 

K 

also  written  as  a  rational  function  of  z.  The  two  functions  are  equated 
which  yields  algebraic  equations  for  unknown  a^'s.  These  equations  are 
then  solved.  From  the  a^'s,  the  refractive  indices  of  each  of  the 
layers  can  be  determined.  Finally,  the  thicknesses  of  the  layers  are 
chosen  so  that  the  system  is  tuned. 
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Below  we  list  some  of  the  major  differences  between  our  holographic 
design  technique  and  the  general  multilayer  dielectric  design  technique. 

(1)  In  the  holographic  case  the  a.'s  are  a  function  of  polariza- 

J 

tion,  incident  angle  and  wavelength.  Therefore  for  a  fixed 
polarization  and  incident  angle,  the  a^s  cannot  be  considered 
as  unknown  constants. 

(2)  The  holographic  design  approach  allows  for  the  systematic 
design  of  detuned  systems,  i.e.,  systems  where  the  $ . 1 s  are  not 

J 

integer  multiples  of  one  another. 

(3)  The  holographic  design  approach  is  based  on  approximate  methods 
which  require  a  limit  to  be  placed  on  the  maximum  value  of 


(4)  The  holographic  design  approach  is  based  on  performing  Fourier 
transforms,  whereas  the  multilayer  dielectric  design  method 

is  based  on  solving  algebraic  equations. 

(5)  Even  if  the  difference  listed  in  (1)  and  (2)  above  did  not 
exist,  the  number  of  a^'s  in  the  holographic  design  problem 
generally  is  much  larger  than  the  number  in  the  multilayer 
design  problem.  This  factor  would  greatly  increase  the  com¬ 
plexity  of  solving  the  required  algebraic  equations  for  the 


We  used  the  first  order  design  technique  to  determine  a  quasi- 
sirosoidal  refractive  index  profile  which  corresponds  to  a  linear  reflec¬ 
tance  vs.  angle  of  incidence  characteristic  (i.e.,  |rR|2  =  20  -  |e0). 
Figure  5-11  shows  the  computed  refractive  index  profile  truncated  to  nom¬ 
inally  a  15  pm  extent  and  Figure  5-12  shows  the  corresponding  plot  of  re¬ 
flectance  vs.  angle  of  incidence.  It  can  be  seen  from  Figure  5-12  that 
the  reflectance  curve  departs  from  the  specified  linear  characteristic 
most  noticeably  in  the  region  around  e0  =  0.  This  could  have  been  antic¬ 
ipated  from  Eq.  (7-5).  Note  that  f  =  2ng  cos  6g/A  varies 
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very  slowly  as  a  function  of  eQ  in  the  region  around  eQ  =  0.  Therefore, 
r^/a  will  vary  quite  rapidly  as  a  function  of  f?  in  this  same  region. 

This  implies  that  the  function  rg('fz)/0l^2)  has  a  significant  amount  of 
high  frequency  content.  It  then  follows  from  Eq.  (7-5)  that  n'(z) 
will  have  a  considerable  amount  of  "energy"  which  lies  outside  of  some 
small  region  around  z  =  0.  Restricting  the  grating  thickness  to  nomi¬ 
nally  15  ym,  as  was  done  in  our  design  calculations,  resulted  in  a  dis¬ 
torted  linear  characteristic.  This  distortion  would  not  have  been  present 
if  this  thickness  restriction  had  not  imposed. 

The  second  order  design  technique  was  used  to  achieve  sidelobe  sup¬ 
pression  in  a  wavelength  selective  filter.  Figure  7-1  is  a  plot  of  re¬ 
flectance  vs.  wavelength  for  a  grating  which  has  a  sinusoidal  refractive 
index  variation  (a.  =  0.02,  f.  =  5.28,  j  =  1,  ....  80).  Notice  the  large 

J  J 

sidelobe  levels.  We  sought  a  design  which  would  have  a  reflectance  (vs. 
wavelength)  identical  to  the  central  lobe  of  the  curve  in  Figure  7-1  and 
which  would  be  zero  elsewhere.  Using  the  2nd  order  design  technique 
(with  the  grating  thickness  restricted  to  nominally  15  ym),  we  arrived  at 
the  refractive  index  profile  shown  in  Figure  7-3  after  four  iterations. 

The  corresponding  reflectance  vs.  wavelength  characteristic  is  shown  in 
Figure  7-2.  For  the  initial  starting  point  of  this  second  order  design 
we  chose  the  a^'s  and  f^'s  to  be  the  same  as  those  of  the  sinusoidal  grat¬ 
ing  (a.  =  0.02,  f.  =  5.28,  j  =  1,  ...,  80).  The  thickness  of  the  layers 

J  %) 

(i.e.,  the  A^'s  k  =  1,  ...,  M)  were  chosen  so  that  maximum  value  of 

| r n  (f,)|^  ~  7%  K  =  1 ,  . ..,  M.  Notice  that  the  sidelobe  level  has  been 
kk  z 

reduced  considerably.  The  width  of  the  reflection,  however,  band  has 
increased.  We  attribute  this  reflection  band  broadening  to  the  thickness 
restriction  which  was  imposed  on  the  grating  design. 

We  wish  to  make  some  qualitative  comments  regarding  the  class  of 
achievable  spectral  and  angular  grating  characteristics.  A  strong  re- 
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Figure  7-1.  Typical  Sinusoidal  Reflection  Grating  Character 
Istlc  Exhibiting  High  Side  Lobe  Levels 
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Figure  7-2.  Quasi -Sinusoidal  Reflection  Grating  Characteristic 
Exhibiting  Reduced  Side  Lobe  Levels 
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Figure  7-3.  Quasi-Sinusoldal  Refractive  Index  Profile  of  a 
Reflection  Grating  Which  Exhibits  Reduced  Side 
Lobe  Levels 
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lationship  has  been  shown  to  exist  between  the  Fourier  transform  and  the 
design  problem.  This  relationship  leads  us  to  make  the  following  con¬ 
jectures: 

2n0  cos  e 

1.  The  rate  at  which  rD  changes  as  a  function  of  f  = - - - - 

K  Z  A 

will  determine  the  nominal  thickness  of  the  grating  required 
to  yield  this  rR(fz). 

2.  If  rR  changes  too  rapidly  as  a  function  of  fz  a  very  thick 
grating  will  be  required. 

3.  If  rR  changes  too  slowly  as  a  function  of  fz,  a  thin  grating 
will  be  required.  If  this  thin  grating  is  to  exhibit  signifi- 

p 

cant  diffraction  efficiency  (i.e.,  |rR(fz)|  is  to  be  large 
for  some  set  of  fz's),  then  its  peak  index  modulation  must  be 
high.  It  may  be  sufficiently  high  that  our  analysis  techniques 
are  no  longer  applicable  or  that  the  grating  cannot  be  produced 
in  currently  available  materials. 
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FABRICATION 

In  the  previous  chapters  we  developed  both  analysis  and  design 
techniques  for  holographic  gratings.  Recall  that  the  class  of  gratings 
considered  were  flat  dielectrics  having  refractive  index  variations 
only  along  the  z  direction  (see  Figure  2-1).  These  refractive  index 
variations  were  quasi-sinusoidal .  This  chapter  explores  techniques 
for  holographically  producing  these  refractive  index  variations  in  a 
recording  material.  We  assume  an  ideal  recording  material  which  is  non¬ 
absorbing  and  which  has  a  linear  response.  The  nonabsorbing  nature  of 
the  material  implies  that  the  optical  field  incident  on  the  material 
during  recording  is  not  attenuated  as  propagates  through  it.  Further¬ 
more,  if  I(z)  designates  the  intensity  of  the  optical  field  within  the 
material  during  recording,  then  from  our  assumption  of  linearity,  the 
resulting  refractive  index  profile,  n(z),  will  be  given  by 

n(z)  =  cQ  +  c1  I (z )  (3-1 ) 

where  Cq  and  are  constants.  It  is  usually  true  that  c-|I(z)  <<  c^, 
and  we  will  assume  this  to  be  the  case.  The  practical  effects  of  de¬ 
partures  from  an  ideal  recording  material  will  be  discussed  later. 

Now  let  us  write  I(z)  as  some  bias,  Ig,  plus  a  spatially  vary¬ 
ing  term  M(z). 

I ( z )  =  I0  +  A I (z )  (8-2) 

Then  Eq.  (8-1)  becomes 

n(z)  =  ( c  0  +  c1  IQ)  +  c^I(z)  (8-3) 
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Similarly  we  can  write  n (z )  as  some  bias,  ng,  plus  a  spatially  varying 
term  n(z). 


n(z)  *=  nQ  +  An(z) 

Comparing  Eqs.  (8-3)  and  (8-4)  yields 
n0  =  C0  +  C1 !0  =  c0 


(8-4) 


(8- 5a) 


Al  ( z )  =  Mtii  (8-5b) 

C1 

But  as  discussed  earlier,  n(z)  is  quasi-sinusoidal .  Therefore  An(z) 
can  be  written  as 

An(z)  =  A(z)  cos  [2tt f ( z ) z]  (8-6) 


where 


4  pm  ^  <  f(z)  <  6  pm  ^ 

and  where  A(z)  and  f(z)  are  slowly  varying  functions  of  z  (i.e.,  A(z) 
and  f(z)  can  be  assumed  constant  over  intervals  of  length  1/4  ^m"1). 
Combining  Eqs.  (8-5b)  and  (8-6)  now  yields 

aKz)  =  hill  cos  [2irf(z)z]  (8-7) 

C1 

Our  goal  then  is  to  holographically  produce  an  optical  field  having  a 
spatially  varying  portion,  Al(z),  given  by  Eq.  (8-7)  above.  Note  that 
the  A I ( z )  of  Eq.  (8-7)  is  only  defined  inside  the  recording  material 
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2FR!E 


which  is  assumed  to  be  of  thickness  D.  We  define  the  extension,  Al(z), 
of  Al  (2 )  by 


Al(z)  = 


Al (z ) ,  0  <  z  <  D 


lg(z),  otherwise 


(8-8) 


where  g(z)  is  some  real  valued  function  which  is  arbitrary  and  which  will  be 
chosen  later. 

Let  Fa  1(B)  denote  the  Fourier  transform  of  aT(z),  that  is 


■/* 


FaI(B)  =  /  A I  (2)  exp  {-i2ir  Bz)  dz 


(8-9) 


Then 


c 

■/ 


A I  ( z )  =  /  Fa  I  ( 3 )  exp  ( i  2tt8z  )  dz 


(8-10) 


Now  Fa  I ( B )  =  FaI(-B)*  since  a  I ( z )  is  real  valued.  Therefore 


a  I ( z )  ■  /  2. Real  Part  of  [FaI(B)  exp  (i2nBz)]  dz 

^0 


(8-11) 


If  we  write  FaI(B)  as 


FAI(B)  =  1 FaI(B) (  exp  tia(B)] 


(8-12) 
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then  Eq.  (8-10)  becomes 

oo 

2 1  FaT(B)  j  cos  (2ttBz  +o(B))dB  (8-13) 

0 


Let  [B  .  ,  B  ]  designate  the  frequency  content  range  of  M(z),  i.e., 
_ mi n  max  J  * 

I  FA  I  ( B )  [  will  be  nominally  zero  for  I  B I  >  Bmflx  or  I  B  [  <  Bml-n.  Then  Eq. 
(8-13)  can  be  written  as 


Al(z) 


2 


FA  I  ( B )  |  cos  (  2ttBz  +  a  (B) )  dB 


(8-14) 


We  choose  the  function  g(z)  so  as  to  minimize  B  -  B  .  . 

JJJo  ^  JJ)  1 1 1 

The  integral  in  Eq.  (8-14)  can  be  approximated  by  the  discrete  sum 
given  below 


(8-15) 


where  <5  and  N  are  constants  such  that  B  +  (N  -  l)6  =  Bmav.  Therefore 

nil  n  max 


N-l 

aT(z)  “  £  25lprr(Bmin  +  j6)l 

j=0 

The  "goodness"  of  the  approximation 
now  the  setup  shown  in  Figure  8-1. 
in  a  holographic  material  using  two 


tos  (2„<Bm1n  ♦  j«)z  *  c(Bm1n  *  os)) 

(8-16) 

increases  as  5  decreases.  Consider 
A  reflection  hologram  is  recorded 
plane  waves.  Both  beams  make  an 
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angle  of  e  .  to  the  normal  of  the  recording  material  where 
3 


9j  =  C0S 


-1 


A(B  •  +  j<5 ) 

A'  min  J  ' 


2n, 


(8-17) 


The  beamsplitter  is  adjusted  so  that  the  beam  intensities,  1^  and  1^, 
are  equal,  and  the  attenuator  is  adjusted  so  that 


hr 1 “5iFiI<w*>i 


(8-18) 


The  phase  shifter  is  set  so  that  the  phase  of  beam  one,  relative  to 

J 

beam  two  at  z  -  0  is  given  by 


‘j  *  °<Bmin  * 


(8-19) 


Then  the  resulting  intensity  pattern,  I.(z),  due  to  the  interference 

3 

between  the  two  beams  will  be  given  by 


'j'2)  ■  ln  *  hi  * 2 


VvS ...  (^P>  • .) 


=  26  1  FAl(Bmin  +  j6)|  +  26|FAl(Bmin  +  j6  | 


cos  (2ir(Bmin  +  J6)z  +  a(Bmin  +  j6)) 


(8-20) 
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Now  by  making  a  series  of  N  different  exposures,  as  described  above  the 
resulting  exposure  will  be  equivalent  to  a  single  exposure  having  an 
intensity  pattern  given  by 


N-l  N-l 

Y.  ‘j(2)  ■!  26ifIT(iu  ♦*>! 

j=0  j=o 


N-l 

Z  26|FZI(Bmin  +  j6)l  cos  (2*(Bmin  +  ^  z  +  o(Bmin  + 
j=0 

(8-21) 

N-l 

Note  that  the  term  ^  26|F^T(Bmin  +  in  Eq’  (B_21 )  Just  a  bias, 
j=0 

Therefore, 


N-l 

Ij(z)  =  bias  +  aT(z)  (8-22) 

j=0 

Thus,  we  have  been  successful  in  producing  holographically  the  desired 
spatially  varying  intensity  pattern. 

We  note  that  a  refractive  index  mismatch  exists  between  the  recording 
material  and  air.  This  mismatch  will  result  in  multiple  reflections  at 
the  air/recording  material  interfaces  during  fabrication.  These  multiple 
reflections  will  result  in  a  distorted  intensity  pattern  which  can  be 
significantly  different  from  the  one  desired.  As  a  consequence,  these 
multiple  reflections  must  be  elminated.  This  is  accomplished  by  immersing 
the  recording  material  in  a  refractive  index  matching  medium. 
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The  above  scheme  for  producing  holographically  the  desired  spatially 
varying  intensity  pattern  has  a  number  of  serious  limitations.  First, 
it  is  extremely  difficult  to  determine  whether  the  phase  shifter  is  set 
appropriately  (i.e.,  whether  =  o(Bm^n  +  j<$)).  Techniques  for  deter¬ 
mining  this  would  most  likely  utilize  reflections  which  occur  at  the  air/ 
recording  material  interfaces.  As  we  have  just  seen  these  reflections 
will  be  non-existant,  since  the  recording  material  must  be  immersed  in  a 
refractive  index  matching  medium.  Secondly,  the  above  scheme  is  quite 
tedious.  It  requires  one  to  make  a  large  number  of  sequential  exposures, 
with  a  critical  alignment  being  performed  before  each  new  exposure. 

An  alternate  approach  for  producing  holographically  aT(z)  is  to  use 
the  setup  shown  in  Figure  8-1  with  e.  fixed  j  =  1,  ....  N,  but  vary  X 
from  exposure  to  exposure.  Thus  with 


ej  =  9  fi  xed 


j  =  1, 


Aj  should  be  given  by 


xj 


;"o  cos  6  fi x€d 
Bmin  + 


(8-23) 


This  is  not  a  suitable  solution  since  the  required  a.'s  as  given  by 

J 

Eq.  (8-23)  will  generally  include  some  wavelengths  for  which  the  recording 

material  is  insensitive  (currently  available  phase  holographic  recording 

materials  are  not  panchromatic).  Furthermore,  a  dye  laser  plus  a  variety 

of  dyes  would  be  necessary  to  generate  all  the  a.'s. 

J 

Consider  now  Figure  8-2.  It  provides  still  a  third  method  for 
producing  holographically  &I(z).  The  set-up  is  equivalent  to  the  one 
shown  In  Figure  8-1  except  beam  number  two  is  produced  from  the  reflec¬ 
tion  off  of  a  mirror  placed  directly  in  back  of  the  recording  material. 
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Recording  Material 


Gap  Filled  With 
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S  polarization 


y 


Figure  8-2.  Single  Beam,  Multiple  Exposure  Interferometric 
Technique  for  Fabricating  Quasl-SInusoldal  Holo¬ 
graphic,  Reflection  Gratings 
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The  phase  of  beam  one  relative  to  that  of  beam  two  Is  determined  by 
the  spacing,  s.,  between  the  recording  material  and  the  mirror.  This 

J 

spacing  will  be  different  for  each  exposure.  This  method  suffers  from 
the  same  problems  as  before,  namely  the  Inability  to  determine  the  phase 
difference  between  the  two  beams  and  the  necessity  of  making  many  ex¬ 
posures. 

Below  we  describe  an  approach  which  eliminates  these  two  problems 
in  a  large  number  of  situations.  First,  choose  g(z)  in  Eq.  (8-8)  so 
that  Al (z )  is  an  even  function  of  z,  l.e., 

ZT(z)  =  aT(-z)  (8-24) 

Then  by  Eq.  (8-9) 

00 

FaT(B)*  =f  aT(z)*  exp  (i2irBz)  dz 


Al(-s)*  exp  (12wB(-c))d(-s) 

(we  made  the  substitution  z  =  -s) 
aT(c)  exp  (-12ttBs)  d£ 

(from  Eq.  (8-24)  and  the  fact  that  ZT(z)  is  a  real 
function) 


=  FAI(B) 


(8-25) 
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Equation  (8-25)  Indicates  that  FaT(B)  Is  a  real  function.  Consequently 
o(Bmin  +  j6 )  in  Eq.  (8-20)  will  be  either  0  or  *.  Equation  (8-21)  can 
now  be  rewritten  as 

N-l  N-l 

Z  'j(2)  •  Z  25!fST<B,nin  *  J«>l  + 

j=0  j-0 

N-l 

Y  2uj6  |  faT( Bmi n  ♦  j«)|  cos  (2T>(Bmin  +  J6 )z)  (8-26) 

where 

A  I1  for  FCT(em1n  *  JO  ;  o 

“j 

-1  otherwise 

Thus,  if  we  use  the  scheme  Illustrated  in  Figure  8-2,  then  the  required 
relative  phase  shift  between  the  two  beams  will  always  be  either  0  or  i. 

A  relative  phase  shift  of  rr  corresponds  to  s^  equals  zero,  whereas  a  rel¬ 
ative  phase  shift  of  zero  corresponds  to  [1]. 

n0(cos  0j>Sj  *  ' 

Now  0 j  is  the  angle  the  incoming  beam  makes  with  the  normal  to  the  record¬ 
ing  material.  This  angle  Is  measured  Internal  to  the  recording  material 
and  hence  cannot  exceed  the  critical  angle,  ec<  The  most  likely  candi¬ 
date  for  a  recording  material  Is  dichromated  gelatin,  which  has  a  nominal 
refractive  index  of  1.5. 
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Therefore 

9c  “  s1"  -  42° 

and  so 

ej  <  42° 

Suppose  s.  j  «  1,  ....  N  is  chosen  such  that 

J 

Y1  nn(cos  30°) s.  =  n 
A0  J 

then 

Y2-  nn(cos  0°)  s .  =  1.15  « 

A0  J 

and 

Y^  nQ(cos  42°) Sj  *  0.86  * 

Therefore  if  s.  j  =  1,  ....  N  is  a  constant  having  a  value 

J 

=  l/4(x0/n0)/(cos  30°)  ^q/(”0  cos  eQ) 

then  for  arvyr  angle,  e^,  between  zero  degrees  and  the  critical  angle,  the 
relative  phase  between  beam  one  and  beam  two  will  lie  within  the  interval 
[0.86  n,  1.15  it].  Thus  the  relative  phase  difference  won't  differ  appre¬ 
ciably  from  ir.  It  should  also  be  noted  that  the  worst  case  situation 
has  analyzed.  In  most  fabrication  problems  the  range  of  e  .  will  be  smaller, 

J 

and  consequently,  the  above  Interval  will  also  be  smaller. 

We  have  shown  that  a  good  approximation  to  the  required  Intensity 
pattern  can  be  generated  using  the  setup  shown  in  Figure  8-2.  Two  sets 
of  exposures  are  made.  The  first  set  corresponds  to  a  <  phase  difference 
between  the  two  beams  (at  z  *  0)  and  is  accomplished  by  placing  the  mir- 
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ror  approximately  a  quarter  wavelength  away  from  the  recording  material . 

It  is  impossible  to  place  a  conventional  mirror  in  direct  contact  with 
the  recording  material  since  neither  the  mirror  nor  recording  material 
can  be  made  perfectly  flat.  We  circumvent  this  problem  by  using  liquid 
mercury  as  the  mirror.  We  propose  that  the  quarter  wave  separation  be¬ 
tween  the  recording  material  and  the  mercury  is  realized  by  coating  the 
recording  material  with  a  thin  layer  of  plastic.  This  is  to  be  done  by 
dip  coating.  The  recording  material  is  dipped  into  a  solvent  containing 
the  plastic,  and  then  is  slowly  withdrawn  at  a  constant  speed[2].  As  it  is 
withdrawn,  it  pulls  along  some  of  the  solvent  which  rapidly  evaporates 
leaving  behind  a  thin  film  of  plastic.  The  thickness  of  this  film  will 
be  a  function  of,  among  other  things,  the  plastic  concentration  in  the 
solvent,  the  ambient  humidity  and  temperature,  and  the  speed  at  which 
the  recording  material  is  withdrawn.  The  thin  film  of  plastic  can  be 
removed  by  redipping  the  recording  material  in  the  solvent  alone. 

We  have  suggested  a  technique  for  avoiding  the  difficult  problem  of 
setting  the  phase  difference  between  the  two  beams.  We  now  wish  to  present 
a  method  which  will  eliminate  the  need  to  make  many  exposures.  With  this 
method  two  exposures  will  be  sufficient.  Consider  the  set  up  shown  in 
Figure  8-3.  The  input  transparency  has  an  amplitude  transmittance  variation 
only  along  the  x  direction.  The  combination  of  spherical  lens  #1  and  cylin¬ 
drical  lens  #1  image  this  transparency,  in  the  form  of  a  slit  along  the  x 
direction,  onto  a  rotating  diffuser.  The  rotating  diffuser  makes  each  point 
on  this  slit  spatially  incoherent  with  any  other.  Therefore,  each  of  these 
points  will  act  as  independent  point  source  radiators.  Spherical  1ms  #2 
collimates  each  of  these  point  sources.  The  collimated  beams  cannot  inter¬ 
fere  with  one  another  since  they  are  derived  from  independent  point  sources. 

A  point  on  the  slit,  having  x  coordinated  xp,  will  result  in  a  col¬ 
limated  beam  making  an  angle  of 
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,  X 

tan  ~jT^~  (8-37 ) 

s2 

measured  from  the  normal  to  the  recording  material.  Furthermore,  this 
point  source  will  have  an  intensity  proportional  to  the  square  of  the 
amplitude  transmittance  of  the  input  transparency  at  x  =  Xp.  Thus  by 
choosing  two  appropriate  input  transparencies  and  making  two  exposures, 
one  with  the  first  transparency  in  the  system  and  a  plastic  film  coating 
on  the  recording  material  and  the  other  with  the  second  transparency 
in  the  system  and  no  coating  on  the  recording  material,  the  intensity 
profile  described  by  Eq.  (8-26)  can  be  generated.  The  first  transparency 
corresponds  to  the  terms  with  u-  =  1  in  that  equation,  while  the  second 

J 

transparency  corresponds  to  terms  with  u.  =  -1 . 

3 

Recall  that  we  originally  approximated  the  integral  of  Eq.  (8-14) 
by  the  discrete  sum  of  Eq.  (8-15).  The  discrete  sum,  of  course,  becomes 
the  integral  as  6  ->  0  and  N  -*■  «.  The  number  N  equaled  the  number  of 
exposures  one  had  to  make.  Consequently,  practical  considerations  put  an 
upper  bound  on  N.  In  the  new  technique  described  above,  however,  there 
is  no  practical  limitation  on  N,  since  each  point  on  the  input  transparency 
represents  one  exposure.  Thus,  the  integral  form  of  Eq.  (8-14)  can  be 
real ized. 

Till  now,  we  have  not  discussed  the  limitations  of  currently  avail¬ 
able  holographic  recording  materials.  This  topic  has  been  extensively 
studied  elsewhere,  and  it  is  not  our  aim  to  explore  this  area  here.  We 
will,  however,  mention  some  of  the  salient  points  which  are  directly 
applicable  to  the  problem  addressed  by  this  report. 

The  holographic  gratings  studied  in  this  report  are  all  volume  phase 
holograms.  The  important  characteristics  of  any  volume  phase  recording 
material  are  its  resolving  power,  its  spectal  sensitivity,  its  maximum 
index  modulation;  its  linearity,  its  dynamic  range,  its  absorbtion,  its 
maximum  thickness,  its  stability  and  its  signal  to  noise  ratio.  A 
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volume  phase  recording  material  suitable  for  reflection  holograms  must 
be  capable  of  resolving  over  5000  line  pairs/mm.  There  are  only  three 
materials  which  have  this  resolving  power.  They  are  dichromated  gelatin 
(DG),  bleached  high  resolution  silver  halide  photographic  plates  (BHRP), 
and  some  photorefracti ve  crystals  such  as  lithium  niobate.  Photorefrac- 
tive  crystals,  however,  are  not  stable,  and  recorded  holograms  degrade 
with  time.  Photorefracti ve  materials  also  have  maximum  refractive  index 
modulations  four  orders  of  magnitude  less  than  either  dichromated  gela¬ 
tin  (~  .05)  or  bleached  photographic  plates  and  are  a  factor  of  six 
orders  of  magnitude  slower  than  bleached  photographic  plates  (three  orders 
of  magnitude  slower  than  dichromated  gelatin).  They  can,  however,  be 
orders  of  magnitude  thicker  than  the  other  two  phase  materials.  Dichro¬ 
mated  gelatin  has  very  high  signal  to  noise  ratio  and  is  far  superior 
in  this  characteristic  to  bleached  photographic  plates.  It’s  less  sen¬ 
sitive  by  some  three  orders  of  magnitude  than  BHRP  and  is  totally  insen- 

o 

sitive  to  wavelengths  longer  than  about  5500  A.  BHRP,  on  the  otherhand, 
can  have  a  panchromatic  response.  The  thickness  of  DG  and  BHRP  are  approx¬ 
imately  the  same  (nominally  15  urn).  Both  materials  exhibit  considerable 
absorbtion.  Therefore,  the  incident  beams  attenuate  as  they  propagate 
into  the  materials,  making  fabrication  of  specified  refractive  index 
profiles  a  difficult  problem.  The  attenuation  appears  to  be  greater  in 
BHRP  than  in  DG.  As  far  as  stability  is  concerned  dichromated  gelatin 
is  quite  stable,  but  BHRP  tends  to  degrade  with  exposure  to  light  (espe¬ 
cially  UV).  No  information  seems  to  be  available  on  the  linearity  or  dy¬ 
namic  range  of  either  dichromated  gelatin  or  bleached  high  resolution  plates. 
These  two  characteristics  will  significantly  affect  the  ability  to  fab¬ 
ricate  specified  quasi-sinusoidal  refractive  index  profiles  in  any  phase 
recording  material.  We  point  out  that  excellent  sinusoidal  reflection 
gratings,  having  almost  100%  diffraction  efficiency,  have  been  fabricated 
in  dichromated  gelatin  [  3].  Thus,  there  is  reason  to  believe  that  it  may 
be  possible  to  fabriate  specified  quasi-sinusoidal  reflection  gratings  in 
0G. 
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Notes  to  Chapter  8 

[1]  A  perfect  mirror  changes  the  phase  of  the  electric  field  (s  polariza¬ 
tion)  upon  reflection  by  u. 

[2]  We  have  had  some  success  with  Foral  105  (Organics  Department, 

Hercules  Incorporated,  Wilmington,  Delaware,  U.S.A.). 

[3]  Holographic  Combiners  for  Head-Up  Displays.  Final  Technical  Report 
No.  AFAL-TR-77-1 1 0,  Radar  and  Optics  Division,  Environmental  Re¬ 
search  Institute  of  Michigan,  Ann  Arbor,  Michigan,  1977. 
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9 

CONCLUSIONS  AND  AREAS  FOR  FUTURE  RESEARCH 

This  report  has  been  a  study  of  holographically  produced  phase  re¬ 
flection  gratings.  The  objectives  have  been  to  obtain  a  fundamental 
understanding  at  optical  propagation  through  such  gratings,  to  develop 
techniques  for  designing  gratings  having  specified  characteristics ,  and 
to  indicate  approaches  for  fabricating  these  gratings.  The  present  re¬ 
search  was  motivated  by  the  possibility  of  holographically  producing 
reflection  grating  structures  that  could  perform  spectral  and  angular 
filtering  operations  which  currently  available  devices  cannot.  In  this 
chapter,  the  report  material  is  briefly  summarized,  and  areas  for  future 
research  are  proposed. 

As  a  first  step  toward  achieving  the  above-mentioned  goals,  we  demon 
strated  the  applicability  of  frequency  domain  concepts  to  holographic 
reflection  gratings.  The  use  of  such  gratings  as  position-invariant 
frequency  domain  filters  was  also  discussed. 

Next,  computationally  efficient  techniques  were  developed  for  analyz 
ing  optical  propagation  through  holographic  reflection  gratings.  Pre¬ 
vious  theories  were  too  restrictive.  They  generally  required  periodic 
refractive  index  variations  and  strong  depletion  of  the  undiffracted 
field.  Furthermore,  they  were  not  computationally  efficient.  We  used 
our  theory  to  determine  some  of  the  basic  properties  of  holographic 
reflection  gratings  such  as  the  effects  of  polarization  and  the  inter¬ 
relationship  between  the  angle  of  incidence  and  wavelength.  Simple 
analytic  expressions  for  the  reflection  and  transmission  coefficients  of 
sinusoidal  reflection  gratings  were  also  derived. 

In  chapter  7  the  design  problem  (or  inverse  scattering  problem)  was 
explored.  We  devised  approximate  techniques  for  determining  the  grating 
refractive  index  profile  required  to  yield  specified  spectral  and  angular 
characteristics.  The  relationship  between  these  techniques  and  design 
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techniques  developed  for  multilayer  dielectric  interference  filters 
was  discussed. 

The  following  chapter  presented  a  nonconventional  interferometric 
approach  for  fabricating  holographic  reflection  gratings  having  quasi - 
sinusoidal  refractive  index  profiles.  Recording  material  considerations 
were  also  mentioned. 

There  are  a  number  of  topics  which  this  report  did  not  address 
which  warrant  further  investigation.  First,  it  would  be  instructive  to 
compare,  in  detail,  Kogelnik's  coupled  wave  theory  predictions  and  our 
analytic  results  for  reflection  gratings  having  sinusoidal  refractive 
index  profiles.  Secondly,  our  design  technique  should  be  thoroughly 
investigated  in  order  to  determine  its  ultimate  limitations.  The  tech¬ 
nique  should  then  be  used  to  specify  the  types  of  angular  and  spectral 
characteristics  which  are  achievable  by  holographic  gratings.  Special 
purpose  holographic  gratings  haveing  quasi-sinusoidal  refractive  index 
profiles  should  also  be  fabricated  and  their  performances  evaluated. 
Finally,  the  connection  between  corrugated,  thin  film  waveguide  gratings 
and  holographic  reflection  gratings  might  be  explored.  Such  waveguide 
gratings  already  have  a  number  of  significant  applications.  The  theory 
developed  in  this  report  may  be  useful  in  analyzing  and  designing  these 
waveguide  gratings. 
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APPENDIX 

Papers  and  technical  presentations  prepared  during  the  second 
year  effort  of  this  contract. 

A  paper  by  K.A.  Winick  entitled  "Holographic  Interference  Filters" 
was  presented  at  the  Annual  Optical  Society  of  America  Meeting,  October 
13-17,  1980,  in  Chicago.  This  paper  will  be  submitted  for  publication 
in  1981. 

This  final  technical  report  will  constitute  the  bulk  of  K.A.  Winick's 
doctoral  dissertation  at  the  University  of  Michigan. 
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